CN114417254A - Method for calculating shale content based on logging data of drilling database - Google Patents
Method for calculating shale content based on logging data of drilling database Download PDFInfo
- Publication number
- CN114417254A CN114417254A CN202111566240.4A CN202111566240A CN114417254A CN 114417254 A CN114417254 A CN 114417254A CN 202111566240 A CN202111566240 A CN 202111566240A CN 114417254 A CN114417254 A CN 114417254A
- Authority
- CN
- China
- Prior art keywords
- shale content
- calculating
- logging data
- database
- target horizon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005553 drilling Methods 0.000 title claims abstract description 32
- 238000011160 research Methods 0.000 claims abstract description 8
- 229910052770 Uranium Inorganic materials 0.000 abstract description 17
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 abstract description 17
- 230000033558 biomineral tissue development Effects 0.000 abstract description 8
- 238000011156 evaluation Methods 0.000 abstract description 4
- 239000010410 layer Substances 0.000 description 12
- 239000004576 sand Substances 0.000 description 6
- 230000002349 favourable effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Agronomy & Crop Science (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Remote Sensing (AREA)
- Evolutionary Biology (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the technical field of uranium mine geological research and uranium resource region evaluation, and particularly relates to a method for calculating shale content based on logging data of a drilling database, which comprises the following steps: step S1: counting the exploration drilling data; step S2: counting the logging data of the target horizon; step S3: calculating the mud content of the target horizon; step S4: and generating a shale content contour map of the target layer. When the shale content in the stratum needs to be counted at any time, the method can quickly, efficiently and accurately count the shale content in the target layer, further identify the control factors of uranium mineralization and effectively overcome the problems of timeliness, low efficiency and low accuracy in the traditional core counting method.
Description
Technical Field
The invention belongs to the technical field of uranium mine geological research and uranium resource region evaluation, and particularly relates to a method for calculating shale content based on logging data of a drilling database.
Background
The formation of sandstone-type uranium mineralization mainly depends on sand bodies with good permeability, the top and the bottom of the sand bodies are provided with mudstones with the function of separation, and the mud content or the sand content can be used as main signs for judging whether the thickness of the sand body of a target layer is developed. The argillaceous content is used as one of key indexes for evaluating the influence of sandstone-type uranium mineralization, and key evaluation parameters can be provided for sandstone-type uranium ore prediction. In previous work, the shale content in the stratum is usually counted by using a core logging method, and the method has certain timeliness, low efficiency and low accuracy. Timeliness, namely, the mudstone content and the sandstone content are recorded and calculated after the drilling engineering is finished, and if the mudstone content is not counted in time, the core cannot be counted after being destroyed; the efficiency is low, namely the method adopts a manual counting mode, so the working efficiency is low; the accuracy is low, i.e. the method needs to ensure that the core sampling rate is high enough, and if the core sampling rate is low, the accuracy of shale content statistics is reduced by deducing the lost lithology.
Disclosure of Invention
The invention aims to provide a method for calculating the shale content based on logging data of a drilling database, which can quickly, efficiently and accurately complete the statistics of the shale content in a target layer when the shale content in a stratum needs to be counted at any time, further identify the control factors of uranium mineralization and effectively overcome the problems of timeliness, low efficiency and low accuracy in the traditional core counting method.
The technical scheme for realizing the purpose of the invention is as follows:
a method of calculating shale content based on borehole database log data, the method comprising the steps of:
step S1: counting the exploration drilling data;
step S2: counting the logging data of the target horizon;
step S3: calculating the mud content of the target horizon;
step S4: and generating a shale content contour map of the target layer.
The step S1 specifically includes: and (4) constructing a drilling database of the research area, and counting coordinates, encountered stratums and comprehensive logging data of different drill holes in the drilling database.
The study area bore database of step S1 includes: a drilling basic information table, a comprehensive record table and a logging record table.
The step S2 specifically includes: and selecting a target horizon, recording the coordinates of the drill holes for all the drill holes of the target horizon revealed by the drilling engineering, and respectively obtaining the drill hole numbers, the drill hole coordinates, the natural gamma values and the like of different drill holes on the target horizon.
The step S3 specifically includes: and summarizing a regression equation of the natural gamma value and the shale content in the comprehensive logging data according to the statistical rule, calculating the shale content of the single drill hole one by one according to the regression equation by using the natural gamma value obtained in the step S2, and calculating the shale content of all drill holes in the target horizon according to all drill holes of the target horizon counted in the drill hole database until the calculation of the shale content of all drill holes in the target horizon is completed.
The regression equation of the natural gamma value and the shale content in the comprehensive logging data in the step S3 is as follows:
y＝a*exp(b*x)，
wherein, y: calculating the mud content result;
x: a natural gamma count API;
a, b: through the correlation coefficient obtained by statistics, the correlation coefficient between the two regions is different.
The step S4 specifically includes: and generating a shale content contour map of the target horizon by adopting a space interpolation method according to the drill hole coordinates of the target horizon revealed by the drilling engineering and the shale content calculated in the step S3.
The invention has the beneficial technical effects that:
1. the invention provides a method for calculating the shale content based on logging data of a drilling database, which is characterized in that a favorable ore-forming structure in an ore-finding target layer is identified through the logging data, a lithologic space combination of mud-sand-mud which is favorable for sandstone-type uranium ore-forming is mainly identified, the thickness of the general upper and lower layers of mudstone is 1-2 m, the lithologic combination relation of the thickness of the middle sand body which is 20-30 m is most favorable for interlayer oxidized zone-type uranium ore-forming, the upper and lower mudstones are barrier layers of confined water, and the middle sand body is a migration space of interlayer water and a favorable space of sandstone-type uranium ore-forming.
2. The method can quickly and effectively acquire the shale content of the target layer, automatically generate the contour map of the shale content of the target layer, effectively identify the control factors of uranium mineralization, has the characteristics of simplicity, quickness and the like, is suitable for screening of sandstone-type uranium mineralization environment, research of uranium mineralization rules and the like, and further lays a foundation for sandstone-type uranium ore prediction and evaluation work.
Drawings
FIG. 1 is a flow chart of a method for calculating shale content based on borehole database logging data according to the present invention;
FIG. 2 shows a first prosodic layer (J) of a first lithologic segment of a Xishan kiln group in the decared beach area of Tuoha basin2X1-1) And (5) a argillaceous content equivalence map.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples.
As shown in fig. 1, the method for calculating the shale content based on the logging data of the borehole database provided by the invention specifically comprises the following steps:
step S1: statistical exploration of borehole data
Constructing a research area drilling database, wherein the database at least comprises the following 3 tables: a drilling basic information table, a comprehensive record table and a logging record table. The drilling basic information table is used for recording the geographical position of drilling and mainly comprises drilling number, drilling coordinates, responsible units, construction dates, operators and other general information; the comprehensive record table is used for recording record data of the drill core, and mainly comprises detailed record information of a drill hole number, a stratum, lithology, granularity and the like; the logging record table is used for recording data acquired by radioactive logging methods, electrical logging methods and other logging methods, and mainly comprises a drilling hole number, a logging point number, an earth surface depth, a natural gamma value (API), resistivity and the like. In the borehole database, the coordinates of different boreholes, the encountered formations, the synthetic log data (natural gamma values (API)), etc. may be counted.
Step S2: statistical target horizon logging data
According to the geological prospecting research result of the research area, selecting a corresponding prospecting target horizon, counting all drilled holes of the target horizon exposed by the drilling engineering in a drilled hole database, and further acquiring data such as corresponding drilled hole numbers, drilled hole coordinates, natural gamma values (API) and the like.
Step S3: calculating the mudness content of the target horizon
Because the natural gamma value (API) and the shale content in the comprehensive logging data have a certain corresponding relation, a regression equation is firstly summarized according to a statistical rule, and the regression equation of the natural gamma value (API) and the shale content in the comprehensive logging data is as follows:
y＝a*exp(b*x)，
wherein, y: calculating the mud content result;
x: a natural gamma count API;
a, b: through the correlation coefficient obtained by statistics, the correlation coefficient between the two regions is different.
And calculating the mud content of each drill hole one by one according to the regression equation by using the natural gamma value (API) obtained in the step S2, and calculating the mud content of all drill holes in the target horizon according to all drill holes of the target horizon counted in the drill hole database until the calculation of the mud content of all drill holes in the target horizon is completed.
Step S4: generating a shale content contour map of a target layer
And (4) according to the drill hole coordinates of the target position revealed by the drilling engineering and the shale content calculated in the step S3, generating a contour map of the shale content of the target position by using a space interpolation method in software by adopting GIS (geographic information System) mapping software.
Taking the 9-15 lines of the decared beach area of the spit-Ha basin as an example, the method comprises the following steps:
step S1: and (4) counting 9-15 line exploration drilling data in the decared beach area of the Tuoha basin.
Step S2: and counting the logging data of the target horizon.
Selecting a first lithologic section first rhythm layer J of the west mountain kiln group in the region2X1-1And (3) as a target horizon development application, counting all drill holes of the stratum revealed by the drilling engineering, recording the coordinates of the drill holes, and respectively obtaining the natural gamma values (API) of different drill holes in the target horizon.
Step S3: computing target horizon J2X1-1The argillaceous content of (a).
Calculating target horizon J using regression equations2X1-1The mud content of (A):
y＝4.66*exp(0.0127*x)
wherein y is the calculated mud content result;
x is natural gamma counting API;
by statistics, the correlation coefficient values of the decared beach area of the Tuhaan basin are respectively 4.66 and 0.0127.
Step S4: and (3) completing all drilling statistics, and generating a shale content contour map of the target layer by adopting ArcGIS software, wherein the contour map reflects that uranium mineralization is mainly positioned in a section with the shale content of less than 15% among 9-15 lines of the decared beach area as shown in figure 2.
The present invention has been described in detail with reference to the drawings and examples, but the present invention is not limited to the examples, and various changes can be made within the knowledge of those skilled in the art without departing from the spirit of the present invention. The prior art can be adopted in the content which is not described in detail in the invention.
Claims (7)
1. A method of calculating shale content based on borehole database log data, the method comprising the steps of:
step S1: counting the exploration drilling data;
step S2: counting the logging data of the target horizon;
step S3: calculating the mud content of the target horizon;
step S4: and generating a shale content contour map of the target layer.
2. The method for calculating the shale content based on the logging data of the borehole database according to claim 1, wherein the step S1 is specifically: and (4) constructing a drilling database of the research area, and counting coordinates, encountered stratums and comprehensive logging data of different drill holes in the drilling database.
3. The method of claim 2, wherein the step of studying the zone borehole database of step S1 comprises: a drilling basic information table, a comprehensive record table and a logging record table.
4. The method for calculating the shale content based on the logging data of the borehole database according to claim 3, wherein the step S2 is specifically: and selecting a target horizon, recording the coordinates of the drill holes for all the drill holes of the target horizon revealed by the drilling engineering, and respectively obtaining the drill hole numbers, the drill hole coordinates, the natural gamma values and the like of different drill holes on the target horizon.
5. The method for calculating the shale content based on the logging data of the borehole database according to claim 4, wherein the step S3 is specifically: and summarizing a regression equation of the natural gamma value and the shale content in the comprehensive logging data according to the statistical rule, calculating the shale content of the single drill hole one by one according to the regression equation by using the natural gamma value obtained in the step S2, and calculating the shale content of all drill holes in the target horizon according to all drill holes of the target horizon counted in the drill hole database until the calculation of the shale content of all drill holes in the target horizon is completed.
6. The method for calculating the shale content based on the logging data of the borehole database as claimed in claim 5, wherein the regression equation of the natural gamma value and the shale content in the synthetic logging data in the step S3 is:
y＝a*exp(b*x)，
wherein, y: calculating the mud content result;
x: a natural gamma count API;
a, b: through the correlation coefficient obtained by statistics, the correlation coefficient between the two regions is different.
7. The method for calculating the shale content based on the logging data of the borehole database according to claim 6, wherein the step S4 is specifically: and generating a shale content contour map of the target horizon by adopting a space interpolation method according to the drill hole coordinates of the target horizon revealed by the drilling engineering and the shale content calculated in the step S3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111566240.4A CN114417254A (en) | 2021-12-20 | 2021-12-20 | Method for calculating shale content based on logging data of drilling database |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111566240.4A CN114417254A (en) | 2021-12-20 | 2021-12-20 | Method for calculating shale content based on logging data of drilling database |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114417254A true CN114417254A (en) | 2022-04-29 |
Family
ID=81268191
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111566240.4A Pending CN114417254A (en) | 2021-12-20 | 2021-12-20 | Method for calculating shale content based on logging data of drilling database |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114417254A (en) |
-
2021
- 2021-12-20 CN CN202111566240.4A patent/CN114417254A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105223612B (en) | A kind of coal mine flood prediction and evaluation method based on earthquake information | |
| CN111090709A (en) | Big data geological analysis method for sandstone-type uranium ore mineralization prediction | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| Harris | The role of geology in reservoir simulation studies | |
| CN116299672B (en) | Fracture-cavity reservoir geomechanical heterogeneity-anisotropy modeling method | |
| Zhang et al. | Application of integrated geology and geomechanics to stimulation optimization workflow to maximize well potential in a tight oil reservoir, Ordos Basin, northern central China | |
| Satur et al. | Characteristics of a wing-like sandstone intrusion, Volund Field | |
| WO2024064095A1 (en) | Automated cell-to-cell calibration of subsidence information map in forward geological models | |
| Dai et al. | Characteristics and formation mechanism of the fractures in Archaean buried hill: A case study in the BZ19‐6 Block, Bohai Bay Basin, China | |
| CN114417254A (en) | Method for calculating shale content based on logging data of drilling database | |
| Olson et al. | Reservoir characterization of the giant Hugoton gas field, Kansas | |
| CN113514879B (en) | Ear layer identification method | |
| CN107664776A (en) | Carbonate rock weathering crust karst reservoir thickness earthquake prediction method and device | |
| Peijs-van Hilten et al. | Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada | |
| CN112394416A (en) | Heterogeneous fracture control reservoir prediction method and device | |
| Li et al. | Multiscale Modeling of Meandering Fluvial Reservoir Architecture Based on Multiple‐Point Geostatistics: A Case Study of the Minghuazhen Formation, Yangerzhuang Oilfield, Bohai Bay Basin, China | |
| Grazulis | Analysis of stress and geomechanical properties in the Niobrara Formation of Wattenberg Field, Colorado, USA | |
| Machado et al. | Characterization of Arbuckle-basement wastewater disposal system, Payne County, Northern Oklahoma | |
| Chapin et al. | Integrated static and dynamic modeling of the Pinedale tight gas field, Wyoming | |
| Wach et al. | Well placement, cost reduction, and increased production using reservoir models based on outcrop, core, well-log, seismic data, and modern analogs: Onshore and offshore Western Trinidad | |
| Budding et al. | Probabilistic modelling of discontinuous reservoirs | |
| Ameri et al. | RESERVOIR CHARACTERIZATION OF UPPER DEVONIAN GORDON SANDSTONE, JACKSONBURG STRINGTOWN OIL FIELD, NORTHWESTERN WEST VIRGINIA | |
| Kulikowski et al. | Exploring for unswept gas: Integrating reservoir models, production data and fault seal analysis in the Northern Carnarvon Basin, Australia | |
| Zeynalov et al. | Petro physical Modeling of the Tertiary Reservoirs in the Small Giant Bai Hassan Oil Field, Northern Iraq | |
| Macaluso et al. | Perched Water-Identification and Production Behavior In A Real Case |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |