3D geological modeling for mineral resource assessment of the Tongshan Cu deposit,Heilongjiang Province,China

Three-dimensional geological modeling （3DGM） assists geologists to quantitatively study in three-dimensional （3D） space structures that define temporal and spatial relationships between geological objects. The 3D property model can also be used to infer or deduce causes of geological objects. 3DGM technology provides technical support for extraction of diverse geoscience information, 3D modeling, and quantitative calculation of mineral resources. Based on metallogenic concepts and an ore deposit model, 3DGM technology is applied to analyze geological characteristics of the Tongshan Cu deposit in order to define a metallogenic model and develop a virtual borehole technology; a BP neural network and a 3D interpolation technique were combined to integrate multiple geoscience information in a 3D environment. The results indicate： （1） on basis of the concept of magmatic-hydrothermal Cu polymetallic mineraliza- tion and a porphyry Cu deposit model, a spatial relational database of multiple geoscience information for mineralization in the study area （geology, geophysics, geochemistry, borehole, and cross-section data） was established, and 3D metallogenic geological objects including mineralization stratum, granodiorite, alteration rock, and magnetic anomaly were constructed; （2） on basis of the 3D ore deposit model, 23,800 effective surveys from 94 boreholes and 21 sections were applied to establish 3D orebody models with a kriging interpolation method; （3） combined 23,800 surveys involving 21 sections, using VC＋＋ and OpenGL platform, virtual borehole and virtual section with BP network, and an improved inverse distance interpolation （IDW） method were used to predict and delineate mineralization potential targets （Cu-grade of cell not less than 0.1%）; （4） comparison of 3D ore bodies, metallogenic geological objects of mineralization, and potential targets of mineralization models in the study area, delineated the 3D spatial and temporal relationship and causal processes among the ore bodies, alteration rock, metallo- genic stratum, intrusive rock, and the Tongshan Fault. This study provides important technical support and a scientific basis for assessment of the Tongshan Cu deposit and surrounding exploration and mineral resources.

3-D Geological Modeling-Concept,Methods and Key Techniques

3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.

Geological characteristics and models of fault-foldfracture body in deep tight sandstone of the second member of Upper Triassic Xujiahe Formation in Xinchang structural belt of Sichuan Basin,SW China

In the second member of the Upper Triassic Xujiahe Formation(T_(3)x_(2))in the Xinchang area,western Sichuan Basin,only a low percent of reserves has been recovered,and the geological model of gas reservoir sweet spot remains unclear.Based on a large number of core,field outcrop,test and logging-seismic data,the T_(3)x_(2) gas reservoir in the Xinchang area is examined.The concept of fault-fold-fracture body(FFFB)is proposed,and its types are recognized.The main factors controlling fracture development are identified,and the geological models of FFFB are established.FFFB refers to faults,folds and associated fractures reservoirs.According to the characteristics and genesis,FFFBs can be divided into three types:fault-fracture body,fold-fracture body,and fault-fold body.In the hanging wall of the fault,the closer to the fault,the more developed the effective fractures;the greater the fold amplitude and the closer to the fold hinge plane,the more developed the effective fractures.Two types of geological models of FFFB are established:fault-fold fracture,and matrix storage and permeability.The former can be divided into two subtypes:network fracture,and single structural fracture,and the later can be divided into three subtypes:bedding fracture,low permeability pore,and extremely low permeability pore.The process for evaluating favorable FFFB zones was formed to define favorable development targets and support the well deployment for purpose of high production.The study results provide a reference for the exploration and development of deep tight sandstone oil and gas reservoirs in China.

A New Method of Multi-Scale Geologic Modeling and Display

A new method of multi-scale modeling and display of geologic data is introduced to provide information with appropriate detail levels for different types of research. The multi-scale display mode employs a model extending existing 2D methods into 3D space. Geologic models with different scales are organized by segmenting data into orthogonal blocks. A flow diagram illustrates an octree method for upscaling between blocks with different scales. Upscaling data from the smallest unit cells takes into account their average size and the Burgers vector when there are mismatches. A geocellular model of the Chengdao Reservoir of the Shengli Oilfield, China is taken as an illustrative case, showing that the methods proposed can construct a multi-scale geologic model correctly and display data from the multi-scale model effectively in 3D.

Construction of carbonate reservoir knowledge base and its application in fracture-cavity reservoir geological modeling

To improve the efficiency and accuracy of carbonate reservoir research,a unified reservoir knowledge base linking geological knowledge management with reservoir research is proposed.The reservoir knowledge base serves high-quality analysis,evaluation,description and geological modeling of reservoirs.The knowledge framework is divided into three categories:technical service standard,technical research method and professional knowledge and cases related to geological objects.In order to build a knowledge base,first of all,it is necessary to form a knowledge classification system and knowledge description standards;secondly,to sort out theoretical understandings and various technical methods for different geologic objects and work out a technical service standard package according to the technical standard;thirdly,to collect typical outcrop and reservoir cases,constantly expand the content of the knowledge base through systematic extraction,sorting and saving,and construct professional knowledge about geological objects.Through the use of encyclopedia based collaborative editing architecture,knowledge construction and sharing can be realized.Geological objects and related attribute parameters can be automatically extracted by using natural language processing(NLP)technology,and outcrop data can be collected by using modern fine measurement technology,to enhance the efficiency of knowledge acquisition,extraction and sorting.In this paper,the geological modeling of fracture-cavity reservoir in the Tarim Basin is taken as an example to illustrate the construction of knowledge base of carbonate reservoir and its application in geological modeling of fracture-cavity carbonate reservoir.

Reservoir geological modeling and significance of Cambrian Xiaoerblak Formation in Keping outcrop area, Tarim Basin, NW China

Take the Cambrian Xiaoerblak Formation in the Keping(Kalpin) outcrop area as an example, a 28 km reservoir scale geological model was built based on description of 7 profiles, observation of more than 1000 thin sections, petrophysical analysis of 556 samples and many geochemical tests. The Xiaoerblak Formation, 158–178 m thick, is divided into three members and 5 submembers, and is composed of laminated microbialite dolomite(LMD), thrombolite dolomite(TD), foamy-stromatolite dolomite(FSD), oncolite dolomite(OD), grain dolomite(GD)/crystalline dolomite with grain ghost and micritic dolomite(MD)/argillaceous dolomite. The petrology features show that its sediment sequence is micro-organism layer – microbial mound/shoal – tidal flat in carbonate ramp background from bottom up. The reservoir has 5 types of pores, namely, framework pore, dissolved vug, intergranular and intragranular dissolved pore and intercrystalline dissolved pore, as main reservoir space. It is found that the development of pore has high lithofacies selectivity, FSD has the highest average porosity, TD, OD and GD come second. The reservoir is pore-vug reservoir with medium-high porosity and medium-low permeability. The dolomite of Xiaoerblak Formation was formed in para-syngenetic to early diagenetic stage through dolomitization caused by seawater. The reservoir development is jointly controlled by sedimentary facies, micro-organism type, high frequency sequence interface and early dolomitization. The classⅠand Ⅱ reservoirs, with an average thickness of 41.2 m and average reservoir-stratum ratio of about 25.6%, have significant potential. It is predicted that the microbial mounds and shoals in the middle ramp around the ancient uplift are the favorable zones for reservoir development.

Optimization of integrated geological-engineering design of volume fracturing with fan-shaped well pattern

According to the variable toe-to-heel well spacing, combined with the dislocation theory, discrete lattice method, and finite-element-method(FEM) based fluid-solid coupling, an integrated geological-engineering method of volume fracturing for fan-shaped well pattern is proposed considering the geomechanical modeling, induced stress calculation, hydraulic fracturing simulation, and post-frac productivity evaluation. Besides, we propose the differential fracturing design for the conventional productivity-area and the potential production area for fan-shaped horizontal wells. After the fracturing of the conventional production area for H1 fan-shaped well platform, the research shows that the maximum reduction of the horizontal principal stress difference in the potential productivity-area is 0.2 MPa, which cannot cause the stress reversal, but this reduction is still conducive to the lateral propagation of hydraulic fractures. According to the optimized fracturing design, in zone-Ⅰ of the potential production area, only Well 2 is fractured, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage;in zone-Ⅱ, Well 2 is fractured before Well 3, with a cluster spacing of 30 m and an injection rate of 12 m^(3)/min per stage. The swept area of the pore pressure drop in the potential production area is small, showing that the reservoir is not well developed. The hydraulic fracturing in the toe area can be improved by, for example, properly densifying the fractures and adjusting the fracture distribution, in order to enhance the swept volume and increase the reservoir utilization.

A 3D attention U-Net network and its application in geological model parameterization

To solve the problems of convolutional neural network–principal component analysis(CNN-PCA)in fine description and generalization of complex reservoir geological features,a 3D attention U-Net network was proposed not using a trained C3D video motion analysis model to extract the style of a 3D model,and applied to complement the details of geologic model lost in the dimension reduction of PCA method in this study.The 3D attention U-Net network was applied to a complex river channel sandstone reservoir to test its effects.The results show that compared with CNN-PCA method,the 3D attention U-Net network could better complement the details of geological model lost in the PCA dimension reduction,better reflect the fluid flow features in the original geologic model,and improve history matching results.

The Method and Practice of Constructing 3D Geological Model from Coalfield Exploration 2D Maps

3D geological modeling is an inevitable choice for coal exploration to adapt to the transformation of coal mining for green, fine, transparent and Intelligent mining. In the traditional Coalfield exploration geological reports, the spatial expression form for the coal seams and their surrounding rocks are 2D maps. These 2D maps are excellent data sources for constructing 3D geological models of coal field exploration areas. How to construct 3D models from these 2D maps has been studying in coal exploration industry for a long time, and still no breakthrough has been achieved so far. This paper discusses the principle, method and software design idea of constructing 3D geological model of an exploration area with 2D maps made by AutoCAD/MapGIS. At first, the paper analyzes 3D geological surface expression mode in 3D geological modeling software. It is pointed out that although contour method has unique advantages in coal field exploration, TIN (Triangular Irregular Network) is still the standard configuration of 3D modeling software for coal field. Then, the paper discusses the method of 2D line features obtaining elevation and upgrading 2D curve to 3D curve. Next, the method of semi-automatic partition is introduced to build the boundary ring of the surface patch, that is, the user clicks and selects the line feature to build the outer boundary ring of the surface patch. Then, Auto-process method for fault line inside of the outer boundary ring is discussed, it including construction of fault ring, determining fault ring being normal fault ring or reverse fault ring and an algorithm of dealing with normal fault ring. An algorithm of dealing with reverse fault ring is discussed detailly, the method of expanding reverse fault ring and dividing the duplicate area in reverse fault into two portions is introduced. The paper also discusses the method of extraction ridge line/valley line, the construction of fault plane, the construction of stratum and coal body. The above ideas and methods have been initially implemented in the “3D modeling platform for coal field exploration” software, and applied to the 3D modeling practice of data from several coal field exploration areas in Ningxia, Shanxi, Qinghai, etc.

Framework system and research flow of uncertainty in 3D geological structure models

Uncertainty in 3D geological structure models has become a bottleneck that restricts the development and application of 3D geological modeling.In order to solve this problem during periods of accuracy assessment,error detection and dynamic correction in 3D geological structure models,we have reviewed the current situation and development trends in 3D geological modeling.The main context of uncertainty in 3D geological structure models is discussed.Major research issues and a general framework system of uncertainty in 3D geological structure models are proposed.We have described in detail the integration of development practices of 3D geological modeling systems,as well as the implementation process for uncertainty evaluation in 3D geological structure models.This study has laid the basis to build theoretical and methodological systems for accuracy assessment and error correction in 3D geological models and can assist in improving 3D modeling techniques under complex geological conditions.

A 3D Geological Model Constrained by Gravity and Magnetic Inversion and its Exploration Implications for the World-class Zhuxi Tungsten Deposit, South China

The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallogeny of the deposit have included its timing,the ore-controlling structures and sedimentary host rocks and their implications for mineral exploration.However,the deep nappe structural style of Taqian-Fuchun metallogenic belt that hosts the W deposit,and the spatial shape and scale of deeply concealed intrusions and their sedimentary host rocks are still poorly defined,which seriously restricts the discovery of new deposits at depth and in surrounding areas of the W deposit.Modern 3 D geological modeling is an important tool for the exploration of concealed orebodies,especially in brownfield environments.There are obvious density contrast and weak magnetic contrast in the ore-controlling strata and granite at the periphery of the deposit,which lays a physical foundation for solving the 3 D spatial problems of the ore-controlling geological body in the deep part of the study area through gravity and magnetic modeling.Gravity data(1:50000)and aeromagnetic data(1:50000)from the latest geophysical surveys of 2016-2018 have been used,firstly,to carry out a potential field separation to obtain residual anomalies for gravity and magnetic interactive inversion.Then,on the basis of the analysis of the relationship between physical properties and lithology,under the constraints of surface geology and borehole data,human-computer interactive gravity and magnetic inversion for 18 cross-sections were completed.Finally,the 3 D geological model of the Zhuxi tungsten deposit and its periphery have been established through these 18 sections,and the spatial shape of the intrusions and strata with a depth of 5 km underground were obtained,initially realizing―transparency‖for ore-controlling bodies.According the analysis of the geophysical,geochemical,and geological characteristics of the Zhuxi tungsten deposit,we discern three principles for prospecting and prediction in the research area,and propose five new exploration targets in its periphery.

Geological Model of Member 3 of Shahejie Formation Reservoir in Liuzan Oilfield, Eastern Hebei Province

A new geological model of Member 3 of Shahejie Formation reservoir in the Liuzan Oilfield, eastern Hebei Province was constructed by using modem reservoir modeling technology as sequence stratigraphy and conditional simulation combined with traditional geological analysis. The model consists of a stratigraphic framework model, a structural model, a sedimentary model and a reservoir model. The study shows that the reservoir is a relatively integrated nose structure, whose strata can be divided into 4 sets of parasequence, 12 parasequences. The submerged branched channel of fan delta front is the favorable microfacies, which controls the geometric shape and physical properties of reservoir sandstone. Oil is distributed in premium reservoir sandstones at structural high positions. According to the new geological model, not only the geological contradictions appearing during oil field development are resolved, but also the oil-bearing area is enlarged by 2.7km^2 and geological reserves increased by 156.9 million tons. The production capacity of the Liuzan Oilfield is increased by 0.27 million tons per year.

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.

A method of reconstructing 3D model from 2D geological cross-section based on self-adaptive spatial sampling:A case study of Cretaceous McMurray reservoirs in a block of Canada

An orthogonal 2D training image is constructed from the geological analysis results of well logs and sedimentary facies;the 2 D probabilities in three directions are obtained through linear pooling method and then aggregated by the logarithmic linear pooling to determine the 3 D multi-point pattern probabilities at the unknown points,to realize the reconstruction of a 3 D model from 2D cross-section.To solve the problems of reducing pattern variability in the 2 D training image and increasing sampling uncertainty,an adaptive spatial sampling method is introduced,and an iterative simulation strategy is adopted,in which sample points from the region with higher reliability of the previous simulation results are extracted to be additional condition points in the following simulation to improve the pattern probability sampling stability.The comparison of lateral accretion layer conceptual models shows that the reconstructing algorithm using self-adaptive spatial sampling can improve the accuracy of pattern sampling and rationality of spatial structure characteristics,and accurately reflect the morphology and distribution pattern of the lateral accretion layer.Application of the method in reconstructing the meandering river reservoir of the Cretaceous McMurray Formation in Canada shows that the new method can accurately reproduce the shape,spatial distribution pattern and development features of complex lateral accretion layers in the meandering river reservoir under tide effect.The test by sparse wells shows that the simulation accuracy is above 85%,and the coincidence rate of interpretation and prediction results of newly drilled horizontal wells is up to 80%.

Theoretical Study on Stochastic Modeling of Combined Gravity-Magnetic-Electric-Seismic Inversion and Its Application

As gravity field, magnetic field, electric field and seismic wave field are all physical fields, their object function, reverse function and compound function are certainly infinite continuously differentiable functions which can be expanded into Taylor (Fourier) series within domain of definition and be further reduced into solving stochastic distribution function of series and statistic inference of optimal approximation. This is the basis of combined gravity-magnetic-electric-seismic inversion of stochastic modeling. It is an uncertainty modeling technology of combining gravity-magnetic-electric-seismic inversion built on the basis of separation of field and source gravity-magnetic difference-value (D-value) trend surface, taking distribution-independent fault system as its unit, depths of seismic and electric interfaces of interests as its corresponding bivariate compound reverse function of gravity-magnetic anomalies and using high order polynomial (high order trigonometric function) approximating to its series distribution. The difference from current dominant inversion techniques is that, first, it does not respectively create gravity-seismic, magnetic-seismic deterministic inversion model from theoretical model, but combines gravity-magnetic-electric-seismic stochastic inversion model from stochastic model; second, after the concept of equivalent geological body being introduced, using feature of independent variable of gravity-magnetic field functions, taking density and susceptibility related to gravity-magnetic function as default parameters of model, the deterministic model is established owing to better solution to the contradiction of difficulty in identifying strata and less test analytical data for density and susceptibility in newly explored area; third, under assumption of independent parent distribution, a real modeling by strata, the problem of difficult plane closure arising in profile modeling is avoided. This technology has richer and more detailed fault and strata information than sparse pattern seismic data in newly explored area, successfully inverses and plots structural map of Indosinian discontinuity in Hefei basin with combined gravity-magnetic-electric-seismic inversion. With development of high precision gravity-magnetic and overall geophysical technology, it is certain for introducing new methods of stochastic modeling and computational intelligence and promoting the development of combined gravity-magnetic-electric-seismic inversion to open a new substantial path.

Study on engineering-hydrogeological problems of hydroelectric project

Engineering-hydrogeological problems arise from the interaction between engineering activities and geological environment, in which rock-soil mass and groundwater are especially important constituents. However, up-to-date research on them is relatively dispersive and simple due to their complexity and lack of comprehensive and systematic study methods. Starting from geological analysis of mechanism to geological model based on geological regularities, the paper predicts the tendency of geological evolvement and puts forward proper measures to solve problems. In this paper, elevated water-sensitive structure in rock-soil mass, which mainly causes engineering hydrogeological problems, and problems in hydropower is discussed based on unique construction in Chinese Western hydropower projects. Engineering hydrogeological problems are reservoir induced earthquakes leakage from reservoir bottom in karst, stability of high slope at reservoir banks, sliding of dam foundation and dam abutment, and confined water at key positions which are introduced and determined by using water-sensitive factors(or structure) according to special hydrogeological conditions.

Quantitative evaluation of lateral sealing of extensional fault by an integral mathematical-geological model

To evaluate the lateral sealing mechanism of extensional fault based on the pressure difference between fault and reservoir, an integral mathematical-geological model of diagenetic time on diagenetic pressure considering the influence of diagenetic time on the diagenetic pressure and diagenetic degree of fault rock has been established to quantitatively calculate the lateral sealing ability of extensional fault. By calculating the time integral of the vertical stress and horizontal in-situ stress on the fault rock and surrounding rock, the burial depth of the surrounding rock with the same clay content and diagenesis degree as the target fault rock was worked out. In combination with the statistical correlation of clay content, burial depth and displacement pressure of rock in the study area, the displacement pressure of target fault rock was calculated quantitatively. The calculated displacement pressure was compared with that of the target reservoir to quantitatively evaluate lateral sealing state and ability of the extensional fault. The method presented in this work was used to evaluate the sealing of F_(1), F_(2) and F_(3) faults in No.1 structure of Nanpu Sag, and the results were compared with those from fault-reservoir displacement pressure differential methods without considering the diagenetic time and simple considering the diagenetic time. It is found that the results calculated by the integral mathematical-geological model are the closest to the actual underground situation, the errors between the hydrocarbon column height predicted by this method and the actual column height were 0–8 m only, proving that this model is more feasible and credible.

Characteristics of dissolved inorganic carbon in produced water from coalbed methane wells and its geological significance

Based on long-term dynamic tracing of dissolved inorganic carbon(DIC)and stable carbon isotope(δ13CDIC)in produced water from 20 coalbed methane(CBM)wells in western Guizhou,the spatial-temporal dynamic variations ofδ13CDIC of the GP well group produced in multi-layer commingled manner were analyzed,and the relationship between the value ofδ13CDIC and CBM productivity was examined.The produced water samples of typical wells in the GP well group were amplified and sequenced using 16S rDNA,and a geological response model ofδ13CDIC in produced water from CBM wells with multi-coal seams was put forward.The research shows that:δ13CDIC in produced water from medium-rank coal seams commonly show positive anomalies,the produced water contains more than 15 species of methanogens,and Methanobacterium is the dominant genus.The dominant methanogens sequence numbers in the produced water are positively correlated withδ13CDIC,and the positive anomaly of v is caused by reduction of methanogens,and especially hydrogenotrophic methanogens.Vertical segmentation of sedimentary facies and lithology in stratum with multi-coal seams will result in permeability and water cut segmentation,which will lead to the segmentation ofδ13CDIC and archaea community in produced water,so in the strata with better permeability and high water cut,theδ13CDIC of the produced water is abnormally enriched,and the dominant archaea is mainly Methanobacterium.In the strata with weak permeability and low water cut,theδ13CDIC of the produced water is small,and the microbial action is weak.The shallow layer close to the coal seam outcrop is likely to be affected by meteoric precipitation,so theδ13CDIC of the produced water is smaller.The geological response model ofδ13CDIC in produced water from multi-coal seams CBM wells in the medium-rank coal reveals the geological mechanism and microbial action mechanism of theδ13CDIC difference in the produced water from the multi-coal seams CBM wells.It also provides effective geochemical evidence for the superimposed fluid system controlled by sedimentary facies,and can also be used for the contribution analysis of the produced gas and water by the multi-layer CBM wells.

Modeling the Sediment Fill of the Upper Troy Pre-Glacial Bedrock Valley, McHenry County, Illinois, USA

The Troy Bedrock Valley (TBV) and its tributary valleys are the principal pre-glacial drainage in southern Wisconsin and northern Illinois, USA. This study focused on the headwaters of a tributary that occurs in McHenry County, IL. Drilling, geophysical surveys, and the analysis of existing geologic and water well data were used to determine the lithologic and geometric characteristics of the sediments that fill the paleovalley. A 3D geologic model of these sediments was then developed in Petrel. More than 65 m of Quaternary sediments filled the paleovalley. The model domain covers approximately 30 km2. The valley drains to the west and meanders, which is distinct from the straight course of the overlying modern Kishwaukee River. The sediments that filled the valley were subdivided into five units. These units include Illinois-age Glasford Formation coarse-grained proglacial outwash and alluvial deposits (GS2, GS1) and fine-grained lacustrine and diamicton deposits (G2 and G1). The Wisconsin-age Henry Formation sand and gravel cap the valley fill, and Cahokia alluvium buries everything.