Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.     

多源地质数据综合三维建模及地质学分析——以青海木里三露天勘探区为例

常规的二维地质图件难以直观地表达煤炭资源的三维空间展布。基于青海木里三露天勘探区的钻孔数据、二维平面地形地质图、煤层底板等高线和地质勘探剖面图,对该地区地形、煤层底板和断层的三维空间形态及温度场进行了重构。三维重建结果清晰明确地展示了三露天勘探区煤层的形态特征、各个沉积旋回的聚煤范围和强度、构造轮廓特征以及地下地温特征等。结果表明,基于多源地质数据的三维可视化技术可直观形象地反映煤层等地质体的形态和属性特征,可有效地分析复杂断层在地质历史中产生的先后顺序及古构造应力场特征。使用三维建模技术重构地质体对沉积构造及地质体属性特征分析是一种行之有效的手段。      

Transferable wavelet method for grain‐size distribution from images of sediment surfaces and thin sections,and other natural granular patterns

In images of sedimentary or granular material, or simulations of binary (two‐phase) granular media, in which the individual grains are resolved, the complete size distribution of apparent grain axes is well‐approximated by the global power spectral density function derived using a Morlet wavelet. This approach overcomes many limitations of previous automated methods for estimating the grain‐size distribution from images, all of which rely on either: identification and segmentation of individual grains; calibration and/or relatively large sample sizes. The new method presented here is tested using: (i) various types of simulations of two‐phase media with a size distribution, with and without preferred orientation; (ii) 300 sample images drawn from 46 populations of sands and gravels from around the world, displaying a wide variability in origin (biogenic and mineralogical), size, surface texture and shape; (iii) petrographic thin section samples from nine populations of sedimentary rock; (iv) high‐resolution scans of marine sediment cores; and (v) non‐sedimentary natural granular patterns including sea ice and patterned ground. The grain‐size distribution obtained is equivalent to the distribution of apparent intermediate grain diameters, grid by number style. For images containing sufficient well‐resolved grains, root mean square errors are within tens of percent for percentiles across the entire grain‐size distribution. As such, this method is the first of its type which is completely transferable, unmodified, without calibration, for both consolidated and unconsolidated sediment, isotropic and anisotropic two‐phase media, and even non‐sedimentary granular patterns. The success of the wavelet approach is due, in part, to it quantifying both spectral and spatial information from the sediment image simultaneously, something which no previously developed technique is able to do.     

Delineation of fluvial sediment architecture of subalpine riverine systems using noninvasive hydrogeophysical methods

River management and restoration measures are of increasing importance for integrated water resources management (IWRM) as well as for ecosystem services. However, often river management mainly considers engineering and construction aspects only and the hydrogeological settings as the properties and functions of ancient fluvial systems are neglected which often do not lead to the desired outcome. Knowledge of the distribution of sediment units could contribute to a more efficient restoration. In this study, we present two noninvasive approaches for delineation of fluvial sediment architecture that can form a basis for the restoration, particularly in areas where site disturbance is not permitted. We investigate the floodplain of a heavily modified low-mountain river in Switzerland using different hydrogeophysical methods. In the first approach, we use data from electromagnetic induction (EMI) with four different integral depths (0.75–6 m) and gamma-spectrometry as well as the elevation data as input for a K-means cluster algorithm. The generated cluster map of the surface combines the main characteristics from multilayered input data and delineates areas of varying soil properties. The resulting map provides an indication of areas with different sedimentary units. In the second approach, we develop a new iterative method for the generation of a geological structure model (GSM) by means of various EMI forward models. We vary the geological input parameters based on the measured data until the predicted EMI maps match the measured EMI values. Subsequently, we use the best matched input data for the GSM generation. The derived GSM provides a 3D delineation of possible ancient stream courses. A comparison with an independent ground penetrating radar (GPR) profile confirmed the delineations on the cluster map as well as the vertical changes of the GSM qualitatively. Thus, each of the approaches had the capacity for detecting sedimentary units with distinct hydraulic properties as an indication of former stream courses. The developed methodology presents a promising tool for the characterization of test sites with no additional subsurface information.     

德国北部盆地晚侏罗世沉积环境及古气候分析*

德国北部盆地上侏罗统广泛发育,但野外露头地层普遍出露不全。Hildesheimer Wald地区Wendhausen 6井和Süntel山地区Eulenflucht 1井完整钻遇了上侏罗统牛津阶和启莫里阶地层,为分析该区晚侏罗世沉积演化过程及其所反映的古环境变化规律提供了丰富的资料。通过岩心描述和岩石薄片镜下观察,根据不同层段的颗粒成分、生物组合特征、沉积结构和构造等特征,在2口井的岩心中共划分出14个岩石类型,分别形成于碳酸盐岩斜坡和三角洲环境。建立了该区牛津阶和启莫里阶垂向沉积演化序列,垂向上由Heersum组到Süntel组,沉积环境逐步由外陆棚、内陆棚、临滨过渡到了开阔台地、潮坪环境,表现出了相对海平面降低的进积过程。同时对不同沉积相中保存较好的以低镁方解石为主要成分的牡蛎壳进行原位Mg/Ca值(古温度指标)测试,得出该区牛津期至启莫里期总体表现出了古气候变暖的趋势,且共有3次气候变暖过程。这一古气候变化与由沉积相分析得出的古气候变化一致,且同苏格兰、俄罗斯台地古温度变化趋势有很好的对应关系,表明古气候是控制该区沉积演化的一个重要因素,且牡蛎壳Mg/Ca值可以做为一个古气候指示指标应用于其他地区的古气候分析中。     

Flood inundation mapping and hazard assessment of Baitarani River basin using hydrologic and hydraulic model

Frequent flood is a concern for most of the coastal regions of India. The importance of flood maps in governing strategies for flood risk management is of prime importance. Flood inundation maps are considered dependable output generated from simulation results from hydraulic models in evaluating flood risks. In the present work, a continuous hydrologic-hydraulic model has been implemented for mapping the flood, caused by the Baitarani River of Odisha, India. A rainfall time-series data were fed into the hydrologic model and the runoff generated from the model was given as an input into the hydraulic model. The study was performed using the HEC-HMS model and the FLO-2D model to map the extent of flooding in the area. Shuttle Radar Topographic Mission (SRTM) 90 m Digital Elevation Model (DEM) data, Land use/Land cover map (LULC), soil texture data of the basin area were used to compute the topographic and hydraulic parameters. Flood inundation was simulated using the FLO-2D model and based on the flow depth, hazard zones were specified using the MAPPER tool of the hydraulic model. Bhadrak District was found to be the most hazard-prone district affected by the flood of the Baitarani River. The result of the study exhibited the hydraulic model as a utile tool for generating inundation maps. An approach for assessing the risk of flooding and proper management could help in mitigating the flood. The automated procedure for mapping and the details of the study can be used for planning flood disaster preparedness in the worst affected area.

     

High resolution facies architecture and digital outcrop modeling of the Sandakan formation sandstone reservoir,Borneo:Implications for reservoir characterization and flow simulation

Advances in photogrammetry have eased the acquisition of high-resolution digital information from outcrops, enabling faster, non-destructive data capturing and improved reservoir modeling. Geocellular models for flow dynamics with in the virtual outcrop in siliciclastic deposits at different sets of sandstone facies architecture remain, however, a challenge. Digital maps of bedding, lithological contrast, spatial-temporal variations of bedding and permeability characteristics make it more easy to understand flow tortuosity in a particular architecture. An ability to precisely model these properties can improve reservoir characterization and flow modeling at different scales. Here we demonstrate the construction of realistic 2 D sandstone facies based models for a pragmatic simulation of flow dynamics using a combination of digital point clouds dataset acquired from LiDAR and field investigation of the Sandakan Formation, Sabah, Borneo.Additionally, we present methods for enhancing the accuracy of outcrop digital datasets for producing high resolution flow simulation. A well-exposed outcrop from the Sandakan Formation, Sabah, northwest Borneo having a lateral extent of 750 m was chosen in order to implement our research approach. Sandstone facies and its connectivity are well constrained by outcrop observations, data from air-permeability measurements, bilinear interpolation of permeability, grid construction and water vector analysis for flow dynamics.These proportions were then enumerated in terms of static digital outcrop model(DOM) and facies model based on sandstone facies bedding characteristics. Flow simulation of water vector analysis through each of the four sandstone facies types show persistent spatial correlation of permeability that align with either cross-bedded orientation or straight with more dispersion high quality sandstone(porosity 21.25%-41.2%and permeability 1265.20-5986.25 mD) and moderate quality sandstone(porosity 10.44%-28.75% and permeability 21.44-1023.33 mD). Whereas, in more heterolithic sandstone(wavy-to flaser-bedded and bioturbated sandstone), lateral variations in permeability show spatially non-correlated patterns over centimeters to tens of meters with mostly of low quality sandstone(porosity 3.4%-12.31% and permeability < 1 mD to 3.21 mD). These variations reflect the lateral juxtaposition in flow dynamics. It has also been resulted that the vertical connectivity and heterogeneities in terms of flow are mostly pragmatic due to the interconnected sandstone rather than the quality of sandstone.     

西沙海槽盆地强限制性中央峡谷水道地震相与内部结构的分段特征

西沙海槽盆地中央峡谷深水水道规模大,下切深,弯曲度低,表现出强限制性与东西分段性特征。充分利用研究区二维地震资料,对其地震相特征进行识别与分段性研究,识别出中央峡谷水道体系五种典型地震相类型,分别为底部滞留沉积、浊积砂体沉积、块体流沉积、泥质沉积、滑塌沉积。水道内部充填序列与限制性强度区相关性密切。通过地震相图、沉积相图和古地貌恢复图分析,发现西沙海槽中央峡谷整体呈东西向展布,水道西宽东窄,周边地势西高东低。基于分段特征与构造断裂带的吻合性,以西沙次凹处作为拐点将峡谷分为东西两段,其中西段NWW向展布,限制性弱;东段NEE向展布,限制性强。根据展布、形态、充填等诸多方面因素,总结出了13种分段特征,其限制性与分段性主要受到构造活动、地貌形态与物源供应等的控制。     

Structural elements of the sedimentary cover in the Deryugin Basin of the Sea of Okhotsk

A map of the tectonic zoning of the sedimentary cover in the Deryugin Basin area was constructed based on the structural data. For this purpose, a structural-tectonic map of the sedimentary cover was compiled using CDP seismic reflection and gravimetry data. Depending on the size, spatial position, and orientation, the structural elements of different orders were grouped into tectonic systems. The latters were subdivided into structural zones consisting of higher order structural elements in the form of relative rises and troughs. The corresponding maps and schemes are also presented in this paper.     

沁水盆地晚古生代煤系层序地层及岩相古地理研究

运用层序地层学及沉积学的方法和原理,根据野外露头、钻井及区域地质资料分析,认为山西沁水盆地煤系主要存在碳酸盐岩台地、碎屑岩浅海和三角洲等沉积相类型.进行层序界面的识别、沉积相的分析和准层序的研究对比,将沁水盆地煤系划分为两个层序,每个层序又分为海侵体系域和高位体系域（缺低位体系域）.在此基础上,以体系域为单位编制了岩相古地理图,并分析了各个体系域的岩相古地理特征.     

塔里木盆地走滑带碳酸盐岩断裂相特征及其与油气关系

通过露头与井下资料的综合分析,塔里木盆地奥陶系碳酸盐岩走滑断裂带断裂相具有多样性,根据内部构造发育程度可以分为断层核发育、断层核欠发育两类。露头走滑带断层核部以裂缝带、透镜体、滑动面等断裂相发育为特征,断裂边缘的破碎带发育裂缝带、变形带。裂缝带主要分布在断层核附近50m的破碎带内,裂缝多开启,渗流性好。断裂核部透镜体发育,在破碎带也有分布,破碎角砾组合的透镜体多致密。滑动面具有平直截切型、渐变条带型等两种类型,多为开启的半充填活动面。变形带多为方解石与碎裂岩充填,破碎带局部部位裂缝与溶蚀作用较发育。利用地震剖面、构造图、相干图等资料可以判识塔里木盆地内部奥陶系碳酸盐岩走滑断裂相的特征及其发育程度,沿走滑断裂带走向上断裂相具有分段性与差异性,根据渗流性可以定性区分高渗透相、致密相区。沿断裂带高渗透相区是碳酸盐岩缝洞体储层发育的有利部位。断裂相的横向变化造成油气分布的区段性,形成高渗透相输导模式、致密相遮挡模式等两类成藏模式。走滑断裂带碳酸盐岩断裂相的特征及其控藏作用对油气勘探开发储层建模具有重要意义。     

四川盆地下寒武统龙王庙组浅水碳酸盐岩台地沉积相特征及模式*

准确建立四川盆地龙王庙组沉积模式和精细刻画颗粒滩时空分布规律,是提高天然气勘探成功率的关键。依据野外露头、钻井、测井和地震等资料,综合分析四川盆地及周缘下寒武统龙王庙组沉积相类型及特征;以四级层序为单元,精细刻画等时地层格架内沉积相时空展布特征,建立沉积模式,揭示颗粒滩发育规律,指出下一步天然气勘探方向。结果表明: (1)四川盆地及周缘龙王庙组发育典型浅水碳酸盐岩镶边台地沉积体系,自盆地西缘向东依次发育混积潮坪、局限台地(颗粒滩、云坪和咸化潟湖等微相类型)、开阔台地(颗粒滩、灰坪等微相类型)、台地边缘和斜坡—盆地等沉积相类型;(2)建立了“一凹两隆,双滩成带”的沉积模式,颗粒滩的发育主要受控于古地貌和相对海平面升降,具体表现为垂向上主要发育于四级层序海退时期,平面上沿台地边缘带和围绕咸化潟湖的相对隆起区呈带状分布。由于颗粒滩是优质储集层发育的有利相带,因此该结果可为四川盆地寒武系龙王庙组下一步天然气勘探提供沉积地质基础。     

单因素分析多因素综合作图法——定量岩相古地理重建

笔者倡导并持续采用单因素分析多因素综合作图法这一定量的岩相古地理学方法论。单因素是能独立地反映某地区、某地质时期、某沉积层段沉积环境某些特征的因素,它的有无或含量的多少均可独立地反映该地区、该层段沉积环境的某些特征,如沉积环境水体的深浅、能量大小、性质等。某沉积层段的厚度、岩石类型、结构组分、矿物成分、化学成分、化石及其生态组合等，均可作为单因素。单因素分析多因素综合作图可分三个步骤：第一，是对各剖面尤其是各基干剖面进行认真的地层学和定量岩石学研究，取得各种第一手的定性和定量资料，尤其是定量资料，了解各剖面各沉积层段的沉积环境特征。第二，在已取得的各剖面的定量资料中，按要求的作图单位，选择出那些能独立地反映其沉积环境特征的因素，即单因素；并把全区各剖面各作图单位的各种单因素的百分含量都统计出来，作出各种相应的单因素图，主要是等值线图。这些单因素图可以从不同的侧面定量地反映该地区该沉积层段的沉积环境。这就是单因素分析。第三，把这些定量的单因素图叠加起来，并结合该地区该沉积层段的其他定量和定性资料，去粗取精，去伪存真，全面分析，综合判断，即可编制出该地区该沉积层段的定量岩相古地理图。这就是多因素综合作图。这一方法论的核心是定量，即以各剖面的定量单因素资料为基础，从各单因素图的分析入手，再通过各单因素图的叠加和综合分析判断，最后作出定量的岩相古地理图。在这种岩相古地理图中，各古地理单元的确定都有确切的定量资料和单因素图为依据。定量是笔者等用单因素分析多因素综合作图法编制出的岩相古地理图的最大特色。这就使岩相古地理图和岩相古地理学发展到了定量的阶段。这在古地理学中是个重大的进展。     

Sedimentary facies and model of shallow water carbonates platform of the Lower Cambrian Longwangmiao Formation in Sichuan Basin

Building the sedimentary model accurately and delineating the spatial-temporal distribution of grain bank are the keys to increase the rate of success in natural gas exploration. Thus,the types and characteristics of sedimentary facies of the Longwangmiao Formation in the Sichuan Basin were studied based on field outcrops,drilling,logging and seismic data. Taking a fourth-order sequence systems tract as a unit,the spatial-temporal distribution characteristics of sedimentary facies within the chronostratigraphic framework were precisely described,and the depositional model of the Longwangmiao Formation was established to reveal the development mechanisms of grain bank and point out the future natural gas exploration. The results showed that: (1)The typical shallow carbonate platform depositional system was developed in the Longwangmiao Formation of the Sichuan Basin and its surrounding areas. Mixed tidal flat,restricted platform(microfacies types including grain bank,dolomitic tidal flat,lagoon,and gypsum lagoon),open platform(microfacies types including grain bank,lime tidal flat),platform margin,and slope-basin were developed from the edge to the interior of the basin;(2)A sedimentary model was established,which was characterized by “one depression,two uplifts,and two grain banks in the form of belts”. The development of grain banks was mainly controlled by palaeo-topography and relative sea-level fluctuation. It was mainly developed in the highstand systems tract of a fourth-order sequence vertically,and distributed as a belt along the platform margin and around the lagoon area horizontally. Our paper pointed out the grain shoal deposits are the material basis for the development of high-quality reservoirs,the results thus provide a sedimentary geological basis for future natural gas exploration of the Cambrian Longwangmiao Formation in the Sichuan Basin.     

Sedimentary Microfacies and Porosity Modeling of Deep-Water Sandy Debris Flows by Combining Sedimentary Patterns with Seismic Data: An Example from Unit I of Gas Field A, South China Sea

Sandy debris flow deposits are present in Unit I during Miocene of Gas Field A in the Baiyun Depression of the South China Sea. The paucity of well data and the great variability of the sedimentary microfacies make it difficult to identify and predict the distribution patterns of the main gas reservoir, and have seriously hindered further exploration and development of the gas field. Therefore, making full use of the available seismic data is extremely important for predicting the spatial distribution of sedimentary microfacies when constructing three-dimensional reservoir models. A suitable reservoir modeling strategy or workflow controlled by sedimentary microfacies and seismic data has been developed. Five types of seismic attributes were selected to correlate with the sand percentage, and the root mean square (RMS) amplitude performed the best. The relation between the RMS amplitude and the sand percentage was used to construct a reservoir sand distribution map. Three types of main sedimentary microfacies were identified: debris channels, fan lobes, and natural levees. Using constraints from the sedimentary microfacies boundaries, a sedimentary microfacies model was constructed using the sequential indicator and assigned value simulation methods. Finally, reservoir models of physical properties for sandy debris flow deposits controlled by sedimentary microfacies and seismic inversion data were established. Property cutoff values were adopted because the sedimentary microfacies and the reservoir properties from well-logging interpretation are intrinsically different. Selection of appropriate reservoir property cutoffs is a key step in reservoir modeling when using simulation methods based on sedimentary microfacies control. When the abnormal data are truncated and the reservoir properties probability distribution fits a normal distribution, microfacies-controlled reservoir property models are more reliable than those obtained from the sequence Gauss simulation method. The cutoffs for effective porosity of the debris channel, fan lobe, and natural levee facies were 0.2, 0.09, and 0.12, respectively; the corresponding average effective porosities were 0.24, 0.13, and 0.15. The proposed modeling method makes full use of seismic attributes and seismic inversion data, and also makes the property data of single-well depositional microfacies more conformable to a normal distribution with geological significance. Thus, the method allows use of more reliable input data when we construct a model of a sandy debris flow.     

基于栅格叠合的沉积物底质图生成方法

文章提出了一种基于GIS空间插值和栅格图叠合技术,分析生成沉积物底质图的方法,并以曹妃甸海区为例,以170个海底表层沉积物样的粒度分析结果为依据,应用本文设计的方法生成这一海区的底质分布图,并用预留的20个采样点进行了精度分析和检验,在此基础上比较了不同的空间插值方法对生成的底质图结果的影响.结果显示出这一方法成图速度快、自动化程度高,但也存在对原始数据精度要求高、插值方法带来误差等问题.在借鉴传统底质图生成方法的基础上,运用三维TIN模型对成图结果进行改进,最终形成一种结合地貌要素的栅格叠合的沉积物底质图生成方法.     

Physical properties of Mesozoic sedimentary rocks from the Perth Basin,Western Australia

The Perth Basin (PB) hosts important aquifers within the Yarragadee Formation and adjacent geological formations with potential for economic exploitation by both geothermal energy and carbon capture and sequestration. Published studies on the reservoir quality of the sedimentary units of the PB are very few. This study reports some petrophysical and lithological characteristics of the sedimentary units of interest for geothermal and geosequestration scenarios and help interpolation toward non-sampled intervals. A new fluvial-dominated lithofacies scheme was developed for the Mesozoic stratigraphy from four wells drilled in the central PB (Pinjarra-1, Cockburn-1, Gingin-1 and Gingin-2) based on grainsize, sorting, sedimentary structures and colour that relate to the environment of deposition. Systematic laboratory measurements of permeability, porosity, and thermal conductivity were conducted on core samples to investigate a variety of lithofacies and depths from these wells. Empirical correlations are established among the different physical properties, indicating encouraging relationships for full PB basin interpolation such as between porosity and permeability, when the samples are grouped into ‘hydraulic units’ defined by a ‘flow zone indicator’ parameter. The common principal controls on the PB thermal conductivity are the pore space arrangement and mineralogical content, which are strongly lithofacies-specific. Therefore, the lithofacies type could be a good first-order discriminator for describing spatial variations of thermal conductivity and then estimate their flow zone indicator.     

浙东沿海地区陆海不同介质采样区地球化学成图方法

针对浙东沿海地区陆海不同介质采样区地球化学图中出现的自然因素,以及非点源污染和点源污染等因素引起的随机变异,提出选择一种适宜的数据处理模型和色区空间分配为内容的陆海不同介质地球化学成图方法,以改善图面效果。     

Geostatistical Mapping with Continuous Moving Neighborhood

An issue that often arises in such GIS applications as digital elevation modeling (DEM) is how to create a continuous surface using a limited number of point observations. In hydrological applications, such as estimating drainage areas, direction of water flow is easier to detect from a smooth DEM than from a grid created using standard interpolation programs. Another reason for continuous mapping is esthetic; like a picture, a map should be visually appealing, and for some GIS users this is more important than map accuracy. There are many methods for local smoothing. Spline algorithms are usually used to create a continuous map, because they minimize curvature of the surface. Geostatistical models are commonly used approaches to spatial prediction and mapping in many scientific disciplines, but classical kriging models produce noncontinuous surfaces when local neighborhood is used. This motivated us to develop a continuous version of kriging. We propose a modification of kriging that produces continuous prediction and prediction standard error surfaces. The idea is to modify kriging systems so that data outside a specified distance from the prediction location have zero weights. We discuss simple kriging and conditional geostatistical simulation, models that essentially use information about mean value or trend surface. We also discuss how to modify ordinary and universal kriging models to produce continuous predictions, and limitations using the proposed models.