Fracture spacing during hydro‐fracturing of cap‐rocks

Layered low permeability rock units, like shales, represent seals or ‘cap‐rocks’ in a variety of geological settings. A continuous increase in the fluid pressure gradients across a virtually impermeable rock layer will ultimately lead to hydro‐fracturing. Depending on the boundary conditions, such fracturing may lead to the formation of a set of sub‐parallel cracks oriented more or less perpendicular to the cap‐rock layer. In this article, we propose a new numerical model that describes interactions between multiple cross‐cutting fractures in an elastic low permeability rock layer. The width of each fracture and the spacing between them are modeled as a force balance between the fluid pressure and the elastic forces in the cap‐rock and between each fracture. The model indicates that the system of fractures evolves toward a spatially periodic steady‐state distribution with a fixed fracture spacing and aperture. The results are similar for incompressible and compressible fluids. The steady‐state conditions depend on only two dimensionless parameters, and the fracture spacing is only weakly dependent on the cap‐rock thickness. This is in contrast to fracturing produced by simple extension of an elastic rock layer beyond the fracture strength, in which case fracture spacing is proportional to layer thickness.     

Mechanisms of strain-rate- and reaction-dependent permeability in the mid-crust of the Southern Alps: insight from three-dimensional mechanical models

We used three-dimensional mechanical modelling to explore the interplay of rheology, modes of permeability creation and fluid flow in the mid-crust of an oblique orogen. We used the central Southern Alps and the Otago Schist of New Zealand to constrain our models. We also compared our models with the magnetotelluric survey along the Rangitata–Whataroa rivers that imaged a U-shaped zone of high conductivity, interpreted as interconnected fluids, beneath the central Southern Alps. Modelling was carried out using the numerical code FLAC^3D. We used a number of simple assumptions: an initially homogeneous starting material, deformation boundary conditions based on the tectonics of the South Island, the capability of the modelled material to develop an anisotropic permeability structure, strain rate and reaction-induced permeability increases, initial saturation and lithostatic pore pressures as a basis for our models. The initial isotropic permeability was 10⁻¹⁸ m². We modelled two possible mechanisms of permeability increase: (i) strain-rate dependent and (ii) reaction dependent. For a strong mid-crust, the models showed enhanced permeability and hence fluid interconnectivity in a symmetric region beneath the model main divide, both ends were turned up towards the brittle–ductile transition and fluid flow was the greatest in the across strike direction. For a weak mid-crust, the region of enhanced permeability was asymmetric and turned up towards the brittle–ductile regime close to the Alpine Fault. The strong mid-crust model reproduces the features that are common to all interpretations of the MT soundings and is our preferred model for the central Southern Alps.     

1995-2010年中国省际人口迁移态势与空间格局演变——基于社会网络分析的视角

采用第五次人口普查、2005年人口抽样调查和第六次人口普查数据,基于社会网络分析的视角,研究了省级尺度下中国人口迁移态势与空间格局演变特征。结果显示:(1)1995-2010年间,中国人口迁移规模迅速增大,省际间人口迁移规模差异明显;上海、浙江人口迁入率显著上升,河南、安徽人口迁出率显著上升。(2)北京、上海、安徽、四川等13省是我国人口迁移网络的核心节点,人口迁移网络表现出中心性、收敛性和地区非均衡性特点。(3)京津地区、长三角地区、广东是我国主要人口辐合流场,安徽、四川、河南、湖南等中部省份是我国主要人口辐散流场。研究表明,不同地区的人口迁移有显著的空间异质性,各地应根据自身人口迁移趋势及在全国人口迁移网络中角色,科学制定区域人口政策与经济社会发展规划。     

Evolution of deformation and fault‐related fluid flow within an ancient methane seep system (Eocene,Varna, Bulgaria)

Faults are often important in fuelling methane seep systems; however, little is known on how different components in fault zones control subsurface fluid circulation paths and how they evolve through time. This study provides insight into fault‐related fluid flow systems that operated in the shallow subsurface of an ancient methane seep system. The Pobiti Kamani area (NE Bulgaria) encloses a well‐exposed, fault‐related seep system in unconsolidated Lower Eocene sandy deposits of the Dikilitash Formation. The Beloslav quarry and Beloslav N faults displace the Dikilitash Formation and are typified by broad, up to 80 m wide, preferentially lithified hanging wall damage zones, crosscut by deformation bands and deformation band zones, smaller slip planes and fault‐related joints. The formation of a shallow plumbing system and chimney‐like concretions in the Dikilitash Formation was followed by at least two phases of fault‐related methane fluid migration. Widespread fluid circulation through the Dikilitash sands caused massive cementation of the entire damage zones in the fault hanging walls. During this phase, paths of ascending methane fluids were locally obstructed by decimetre‐thick, continuous deformation band zones that developed in the partly lithified sands upon the onset of deformation. Once the entire damage zone was pervasively cemented, deformation proceeded through the formation of slip planes and joints. This created a new network of more localized conduits in close vicinity to the main fault plane and around through‐going slip planes. ¹³C‐depleted crustiform calcite cements in several joints record the last phase of focused methane fluid ascent. Their formation predated Neogene uplift and later meteoric water infiltration along the joint network. This illustrates how fault‐related fluid pathways evolved, over time, from ‘plumes’ in unconsolidated sediments above damage zones, leading to chimney fields, over widespread fluid paths, deflected by early deformation structures, to localized paths along fracture networks near the main fault.     

Relationship between deformation bands and petroleum migration in an exhumed reservoir rock,Los Angeles Basin,California, USA

An oil‐bearing sandstone unit within the Monterey Formation is exposed in the Los Angeles Basin along the Newport‐Inglewood fault zone in southern California. The unit preserves structures, some original fluids, and cements that record the local history of deformation, fluid flow, and cementation. The structures include two types of deformation bands, which are cut by later bitumen veins and sandstone dikes. The bands formed by dilation and by shear. Both types strike on average parallel to the Newport‐Inglewood fault zone (317°–332°) and show variable dip angles and directions. Generally the older deformation bands are shallow, and the younger bands are steep. The earlier set includes a type of deformation band not previously described in other field examples. These are thin, planar zones of oil 1–2 mm thick sandwiched between parallel, carbonate‐cemented, positively weathering ribs. All other deformation bands appear to be oil‐free. The undeformed sandstone matrix also contains some hydrocarbons. The oil‐cored bands formed largely in opening mode, similar to dilation bands. The oil‐cored bands differ from previously described dilation bands in the degree of carbonate cementation (up to 36% by volume) and in that some exhibit evidence for plane‐parallel shear during formation. Given the mostly oil‐free bands and oil‐rich matrix, deformation bands must have formed largely before the bulk of petroleum migration and acted as semi‐permeable baffles. Oil‐cored bands provide field evidence for early migration of oil into a potential reservoir rock. We infer a hydrofracture mechanism, probably from petroleum leaking out of a stratigraphically lower overpressured reservoir. The deformation bands described here provide a potential field example of a mechanism inferred for petroleum migration in modern systems such as in the Gulf of Mexico.     

Computing geochemical mass transfer and water/rock ratios in submarine hydrothermal systems: implications for estimating the vigour of convection

A method of calculating chemical water/rock ratios is presented that enables the estimation of fluid velocities in open, flow‐through hydrologic systems. The approach is based on relating the gain/loss of a chemical species per kilogram of solid phase to the loss/gain of that species in the fluid phase, integrated across a specified length of the flowpath. After examining the underlying approximations of the approach using a one‐dimensional model of seawater moving through a basalt under nonisothermal conditions, the method is applied to representative zones within a two‐dimensional hydrothermal convective system. The method requires that regions within the flow system can be identified in which the direction of flow is steady for an extended period of time. Estimates of fluid velocity are spatial and temporal averages for the length of the flowpath used in the calculation. The location within the flow system and the nature of the alteration reactions determine which species can provide reliable values of the chemical water/rock ratio and useful estimates of fluid velocities. Over the length of the flowpath considered, the calculation of water/rock ratios works best when a species is controlled by a single reaction. Accurate estimates are favoured if the concentration profile of a species along the flowpath increases or decreases monotonically. If the length of the flowpath extends over more than one reaction zone, then erroneous estimates of the water/rock ratio and fluid velocity are more likely. Model calculations suggest that the quartz/silica system should provide reliable estimates for fluid velocity under a wide range of temperature and flow conditions, in particular in those regions of a system at or near quartz equilibrium, so that the aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.     

初论人性时空场

在缜密论述了人地关系中人性时空场基本概念和特征的基础上,从人性时空场意识流与人地关系的相互作用的角度认知了人性时空场的科学哲学内涵。最后分析了人性时空场意识流的随机变化特点之表达。提出:随着行为地理学研究深入,人性时空场意识流将成为它新的研究对象;与此同时,本文首次导出潜意识地理学(Potential Sense Geography)的学科概念,认为潜意识地理学会弥补行为地理学的不足,它将发挥人文关怀的当代学科优势,继而为人性的完善和人的可持续发展及人地关系的根本和谐作出贡献。     

Experimental study of the effect of variation in in‐situ stress on capillary residual trapping during CO2 geo‐sequestration in sandstone reservoirs

During a geo‐sequestration process, CO₂ injection causes an increase in reservoir pore pressure, which in turn decreases the reservoir net effective stress. Changes in effective stress can change all the reservoir and cap‐rock properties including residual saturations. This article presents the results of an experimental work carried out to understand the potential change in the volumes of residually trapped CO₂, while the porous medium tested underwent change in the net effective stress under in‐situ reservoir conditions of pore pressure and temperature. The experimental results obtained show that an initial 1725 psi (11.9 MPa) decrease in the net effective pressure caused 1.4% reduction in the volumes of residually trapped CO₂, while another 1500 psi (10.3 MPa) reduction caused a further 3.2% drop in the residual saturation of CO₂.