A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

Abstract

Non-Darcian flow in a finite fractured confined aquifer is studied. A stream bounds the aquifer at one side and an impervious stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and porous blocks which mainly store water. Unsteady flow in the aquifer due to a sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by an interaction term that controls the quasi-steady-state fracture—block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical solutions give the same results when the Reynolds number is less than 0.1. The effect of nonlinearity on the flow appears when the Reynolds number is greater than 0.1. The higher the rate of flow from the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters on the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.     

南极洲东部普里兹湾海域重磁场特征及地壳结构

普里兹湾位于南极洲东部大陆边缘,其深部地壳结构特征对认识白垩纪冈瓦纳古陆裂解和新生代大陆边缘形成具有重要意义.本文利用重磁、多道反射地震、声纳浮标折射地震和ODP钻井数据对普里兹湾海域的深部地壳结构进行了研究.研究结果显示,普里兹凹陷表现为典型的盆地负重力异常特征,其沉积基底较深,而在四夫人浅滩为高幅重力正异常,其沉积基底普遍抬升.在大陆架中部存在SW-NE向条带状基底的抬升,且呈朝NE向逐渐变深的趋势.在中大陆架外侧,均衡残余重力异常呈V字形负异常条带状分布,其两翼分别与四夫人浅滩和弗拉姆浅滩外的大陆坡相连.该异常带在大陆架中部向陆的偏移可能是由于古大陆架边缘的地形影响,推测其与普里兹冲积扇同属于洋陆过渡带向陆的部分,在重力模拟剖面表现为地壳向海逐渐减薄.普里兹冲积扇的地壳厚度较薄,平均为6 km,最薄处可达4.6 km,并且根据洋陆过渡带向海端的位置,推测可能属于接近洋壳厚度的过渡壳.重力异常分区的走向与兰伯特地堑在普里兹湾的构造走向基本一致,可能主要反映了二叠纪-三叠纪超级地幔柱对普里兹湾的裂谷作用的影响.该区域的自由空间重力异常和均衡残余异常均表现为超过100×10^-5m/s²的高幅正异常特征,可能由位于大陆架边缘的巨厚沉积体负载在高强度岩石圈之上的区域挠曲均衡作用所导致,可能与该区域第二期裂谷期之后的沉积间断以及快速进积加厚的演化过程有关.普里兹湾磁力异常的走向与重力异常明显不同,大致可分为东北高幅正异常区和西南低幅异常区.重磁异常在走向上的差异反映高磁异常主要来源于岩浆作用形成的铁镁质火成岩的影响,并且岩浆作用的时代不同于基底隆升的时代,而可能形成于前寒武纪或者南极洲和印度板块裂谷期间(白垩纪).     

INTERPLANETARY DUST AND CARBONACEOUS METEORITES: CONSTRAINTS ON POROSITY, MINERALOGY AND CHEMISTRY OF METEORS FROM RUBBLE-PILE PLANETESIMALS

This paper discusses measured textures, porosity, chemical compositions and the minerals of interplanetary dust and meteorites from primitive planetesimals and how this information can be used to explain some of the observed physical and chemical properties of meteor phenomena.     

确定恒星表面有效温度的曲面拟合方法

恒星表面有效温度是恒星的一个重要物理参量，是恒星光谱差异的决定因素。本文提出了一种确定恒星表面有效温度的曲面拟合方法，所使用的拟合曲面模型是多项式的指数函数。首先对历史光谱数据进行PCA处理，再根据PCA特征数据与其表面温度的对应关系计算拟合曲面。通过实验，我们发现使用2维PCA数据和指数为3次多项式，根为10的指数函数模型所得到的拟合曲面，不仅有效好的拟合精度而且有很好的鲁棒性。本文的研究结果对恒星表面有效温度的自动测量具有重要的意义。     

A coupled transport and reaction model for long column experiments simulating bank filtration

Within the scope of the interdisciplinary Natural and Artificial Systems for Recharge and Infiltration research project dealing with riverbank filtration processes at the Berlin water works, a semi‐technical column experiment has been ongoing since January 2003. Here, a 30 m long soil column is infiltrated with surface water sampled from Lake Tegel (Berlin, Germany) under saturated flow conditions. Changes in pore water hydrochemistry sampled on 21 non‐equidistant distributed points are verified by coupled transport and reaction modelling. The objective of reactive transport modelling was to identify the main biogeochemical processes within the soil column during the flushing experiment as a conceptual model for riverbank filtration. Modelling was done with a combination of MATLAB and PHREEQC. The main processes identified are: (1) biogeochemical degradation due to interaction of natural surface water with the soil matrix; (2) continuous dissolution of refractory air bubbles from the soil column matrix. Copyright © 2006 John Wiley & Sons, Ltd.     

济阳坳陷陡坡带深层砂砾岩体次生孔隙成因机制探讨

从岩芯、薄片、扫描电镜以及FMI成像测井等资料分析发现,济阳坳陷陡坡带深层砂砾岩体储层次生孔隙发育,大量发育的次生孔隙不仅增大了储层的有效储集空间,更重要的是改变了储层的孔隙结构,使储层孔喉增大,大大提高了储层的渗透性能。研究认为长石、岩屑、碳酸盐和硫酸盐等酸性组分的溶蚀以及构造运动是深层次生孔隙形成的主导因素。垂向上存在两个次生孔隙发育带,对应深度范围从上至下依次为3450～3850 m和4075～4550 m。其主要形成机理有：①有机质热演化过程中产生的酸性流体的溶蚀作用;②去石膏化和去白云化以及硫酸盐热化学作用;③异常高压加速有机质热演化和对已形成孔隙的保护作用;④构造作用、脱水收缩作用等。     

Effective precipitation in southern Spain ( 266 to 46 ka) based on a speleothem stable carbon isotope record

We present the longest-duration directly dated terrestrial palaeoclimate record from the western Mediterranean region: a flowstone speleothem from Gitana Cave, southeast Spain. The main phase of growth was 274 to 58 ka, dated by multi-collector inductively coupled plasma mass spectrometry (MC-ICPMS) U-series methods. Effective precipitation, which we consider primarily responsible for flowstone calcite δ¹³C variations, measured at 300 μm resolution, was higher during interglacials associated with marine oxygen isotope stages (MIS) 7 and 5, and lower during glacial MIS 6. There is a close correspondence between speleothem δ¹³C and sea surface temperature (SST) estimates from adjacent Atlantic Ocean cores during MIS 6, which implies that oceanic conditions are critical in controlling the western Mediterranean terrestrial moisture balance during glacial periods. Other features of our record, such as the sequence of termination II warming/moistening between approximately 133 and 127 ka, including a “pause” around 130–128 ka, and the lagged termination of MIS 5 warm intervals (5e, 5c and 5a) are similar to other terrestrial records within the Mediterranean basin, indicating climate synchroneity along the northern Mediterranean coast. The Gitana cave region also may have been a refugium for temperate species during short-lived cold/arid periods during MIS 5.     

Shallow-burial dolomite cement: a major component of many ancient sucrosic dolomites

Dolomite cement is a significant and widespread component of Phanerozoic sucrosic dolomites. Cements in dolomites that were never deeply buried are limpid, have planar faces (non‐saddle forms), often distinct zonation in cathodoluminescence and form syntaxial overgrowths on crystals facing pores. Five samples of sucrosic dolomites, interpreted as having had mostly lime‐mudstone or wackestone precursors in four carbonate aquifers, provide insights into the abundance of planar cements in sucrosic dolomites. Such cement comprises 11% to 45% (32% mean) of peritidal to sub‐tidal dolomites on an outcrop in the Edwards aquifer (Early Cretaceous) of central Texas; 19% to 33% (25% mean) of ramp dolomites in the Hawthorn Group (Oligo‐Miocene) and 50% to 70% in shelf dolomites of the Avon Park Formation (Eocene) in the Upper Floridan aquifer of sub‐surface peninsular Florida; 18% to 45% (32+% mean) of sub‐tidal shelf dolomites in quarry sections of the Burlington‐Keokuk Formation (Early Mississippian) in south‐eastern Iowa; and 18% to 76% (50% mean) in shallow cores and outcrops of outer‐shelf dolomites from the Gambier Limestone (Oligo‐Miocene) of South Australia. Backstripping the cement phases revealed by cathodoluminescence colour photomicrographs documents the effects of cements on textural coarsening, pore‐space reduction, induration and general ‘maturation’ of these dolomites. Most pre‐Holocene dolomites are multiphase crystalline rocks composed of: (i) seed crystals or ‘cores’; (ii) crystal cortices that concentrically enlarged the cores; and (iii) free‐space, syntaxial precipitates of limpid cement around the crystals. Remaining CaCO₃ grains and micrite can be replaced by dolomite, but typically they are dissolved between stages (ii) and (iii), creating systems of intercrystal and mouldic pores typical of sucrosic dolomites. Networks of cement overgrowths, aided by water‐filled pore systems under hydrostatic to lithostatic pressure, are judged to slow or prevent compaction in sucrosic dolomites. It can be argued that cortex growth involves both replacement of CaCO₃ particles and microcementation of their interparticle pores. This interpretation, and the abundance of cements in so many dolomites, would obviate the controversy over the volumetrics of ‘replacement dolomitization’. Limpid, planar and syntaxial dolomite cements of early diagenetic origin are interpreted to have precipitated from clear pore waters, at low temperatures (<30 to 35 °C) and shallow burial depths (<100 m), in water‐saturated networks of dolomite ‘silt’ and ‘sand’. Cements in many dolomites in island and continental–aquifer systems appear to result from event‐driven processes related to sea‐level highstands. Cementation events can follow ‘replacement dolomitization’ events by time intervals ranging from geologically ‘instantaneous’ to tens of million years.     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.     

克拉苏构造带盐下超深层储层裂缝组合模式及分布规律

克拉苏构造带盐下超深层储层属于裂缝性低孔砂岩储层,构造裂缝普遍发育,分布复杂且对单井产量影响较大。目前针对该地区储层裂缝的研究已投入大量的工作,主要体现在裂缝的描述、表征和预测等方面,但由于该地区地震资料分辨率低、岩心资料少、成像测井多解性强,裂缝的研究受到一定程度的限制,生产中仅应用了裂缝密度、开度等统计数据。笔者综合野外露头、岩心和测井等资料,从裂缝组合的角度出发,系统开展了该地区储层裂缝及其组合特征研究。发现构造裂缝存在主次之分,主裂缝规模相对较大,一般垂直于地层,多为张性或张剪性,次级裂缝一般是主裂缝的伴生缝。裂缝之间存在多种接触关系,主要表现为"T"、"λ"、"y"、"τ"字等形状。依据裂缝之间的排列接触关系,建立了该地区储层裂缝四种组合模式,即"平行状"、"扫帚状"、"调节状"、"雁列状"。结合各裂缝组合成因的差异和断背斜形成过程中的构造演化过程,划分了裂缝组合平面三类成因区,即调节成因区、褶皱成因区、断裂成因区,与之对应的是高产井部署的一类区、二类区、三类区。