CN116109691A - Reserve calculating method and device based on three-dimensional model and interpolation - Google Patents
Reserve calculating method and device based on three-dimensional model and interpolation Download PDFInfo
- Publication number
- CN116109691A CN116109691A CN202310118882.0A CN202310118882A CN116109691A CN 116109691 A CN116109691 A CN 116109691A CN 202310118882 A CN202310118882 A CN 202310118882A CN 116109691 A CN116109691 A CN 116109691A
- Authority
- CN
- China
- Prior art keywords
- target
- dimensional model
- dimensional
- model
- coordinates
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005553 drilling Methods 0.000 claims abstract description 63
- 238000004364 calculation method Methods 0.000 claims abstract description 49
- 239000011435 rock Substances 0.000 claims abstract description 27
- 230000015572 biosynthetic process Effects 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000000465 moulding Methods 0.000 claims description 3
- 239000003245 coal Substances 0.000 description 45
- 238000005755 formation reaction Methods 0.000 description 36
- 230000008569 process Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 239000002585 base Substances 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000012458 free base Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The disclosure provides a reserve calculation method and device based on a three-dimensional model and interpolation, comprising the following steps: the method comprises the steps of obtaining drilling data of a plurality of drilling holes in a target area, wherein the drilling data comprise drilling start coordinates, ending coordinates and rock stratum characteristic data, constructing a target three-dimensional model of the target area according to the starting coordinates and the ending coordinates, carrying out interpolation operation on the target three-dimensional model based on the target characteristic data in the rock stratum characteristic data to calculate the average quantity of target characteristics, and calculating reserves of the target characteristics in the target area according to the average quantity and the volumes of the target three-dimensional model, so that the accuracy of reserve calculation can be improved, and the reserves can be calculated in real time according to the drilling hole data along with exploitation, so that the timeliness of calculation is improved.
Description
Technical Field
The disclosure relates to the technical field of geological exploration, in particular to a reserve calculation method and device based on a three-dimensional model and interpolation.
Background
In the current coal mine design, boundaries are defined in a plan view by contour lines, fault boundaries, synclines and anticlines according to exploration data, and reserves of pure coal, gas and the like are calculated through an approximation algorithm. However, when reservoir prediction and statistics are performed on a rock and coal seam area within a given range, the rock area and volume of the area can only be estimated by using contour lines, and the calculation result is not accurate. In addition, the coal mine design is full life, along with the construction of a mining area, the mining of the coal field, the drilling holes in the mining area are changed from less to more, and the collected various rock and coal seam characteristic data are also more and more abundant. The coal and rock stratum model which is the basis of the coal mine design needs to be changed along with the change of data. However, in the actual situation, updating of drawings such as exploration geological maps used in coal mine design often has long time delay, so that the calculation timeliness is affected.
Disclosure of Invention
The disclosure provides a reserve calculation method and device based on a three-dimensional model and interpolation, which aim to solve one of the technical problems in the related art at least to a certain extent.
An embodiment of a first aspect of the present disclosure provides a reserve calculating method based on a three-dimensional model and interpolation, including: acquiring drilling data of a plurality of drilling holes in a target area, wherein the drilling data comprise drilling start coordinates, termination coordinates and formation characteristic data; constructing a target three-dimensional model of a target area according to the initial coordinates and the termination coordinates; performing interpolation operation on the target three-dimensional model based on target feature data in the rock stratum feature data to calculate average quantity of the target features; and calculating the reserve of the target feature in the target area according to the average quantity and the volume of the target three-dimensional model.
Embodiments of a second aspect of the present disclosure provide a reserve calculation apparatus based on a three-dimensional model and interpolation, including: the system comprises an acquisition module, a calculation module and a calculation module, wherein the acquisition module is used for acquiring drilling data of a plurality of drilling holes in a target area, and the drilling data comprise drilling start coordinates, termination coordinates and formation characteristic data; the building module is used for building a target three-dimensional model of the target area according to the initial coordinate and the termination coordinate; the first calculation module is used for carrying out interpolation operation on the target three-dimensional model based on target feature data in the rock stratum feature data so as to calculate the average quantity of the target features; and the second calculation module is used for calculating the reserve of the target feature in the target area according to the average quantity and the volume of the target three-dimensional model.
Embodiments of a third aspect of the present disclosure provide a computer device comprising: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to perform a three-dimensional model and interpolation-based formation reservoir calculation method of an embodiment of the present disclosure.
A fourth aspect embodiment of the present disclosure proposes a non-transitory computer-readable storage medium storing computer instructions for causing the computer to perform the formation reserves calculation method based on the three-dimensional model and interpolation disclosed in the embodiments of the present disclosure.
In this embodiment, by acquiring drilling data of a plurality of holes in a target area, where the drilling data includes a drilling start coordinate, a drilling end coordinate, and formation feature data, constructing a target three-dimensional model of the target area according to the start coordinate and the drilling end coordinate, performing interpolation operation on the target three-dimensional model based on the target feature data in the formation feature data to calculate an average amount of target features, and calculating reserves of the target features in the target area according to the average amount and a volume of the target three-dimensional model, the accuracy of reserve calculation can be improved, and the reserves can be calculated in real time according to the drilling data as exploitation proceeds, thereby improving timeliness of calculation.
Additional aspects and advantages of the disclosure will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the disclosure.
Drawings
The foregoing and/or additional aspects and advantages of the present disclosure will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a flow diagram of a reserve calculation method based on a three-dimensional model and interpolation provided in accordance with an embodiment of the present disclosure;
FIG. 2 is a flow diagram of a reserve calculation method based on a three-dimensional model and interpolation provided in accordance with another embodiment of the present disclosure;
FIG. 3 is a schematic diagram of a reserve volume overall calculation flow based on a three-dimensional model and interpolation provided in accordance with an embodiment of the present disclosure;
FIG. 4 is a schematic diagram of a reserve computing device based on a three-dimensional model and interpolation provided in accordance with another embodiment of the present disclosure;
fig. 5 illustrates a block diagram of an exemplary computer device suitable for use in implementing embodiments of the present disclosure.
Detailed Description
Embodiments of the present disclosure are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary only for explaining the present disclosure and are not to be construed as limiting the present disclosure. On the contrary, the embodiments of the disclosure include all alternatives, modifications, and equivalents as may be included within the spirit and scope of the appended claims.
It should be noted that, the execution body of the formation reserve calculation method based on the three-dimensional model and interpolation in this embodiment may be a formation reserve calculation device based on the three-dimensional model and interpolation, where the device may be implemented by software and/or hardware, and the device may be configured in an electronic device, and the electronic device may include, but is not limited to, a terminal, a server, and the like.
FIG. 1 is a flow chart of a method for calculating formation reserves based on a three-dimensional model and interpolation according to one embodiment of the present disclosure, as shown in FIG. 1, the method includes:
s101: drilling data of a plurality of drilling holes in a target area is acquired.
In practical applications, the geological exploration process usually needs to perform perforation sampling, and analyze the sampled data to obtain the drilling data of each drilling hole and store the drilling data in a database. In some embodiments, the borehole data includes, for example, borehole start coordinates, end coordinates, formation characteristic data, and any other possible data, without limitation.
Wherein data related to the formation within the borehole may be referred to as formation characteristic data. In one embodiment, formation characteristic data from geological exploration of a coal mine includes, for example: the rock stratum type (coal seam, gangue layer, mud sand layer, etc.), rock stratum thickness, coal volume weight, gas amount, water content, float coal bonding index, raw coal dry base ash, float coal dry ash-free base volatile matter, raw coal dry base total sulfur, raw coal dry base low-position heating value, poisson ratio, porosity, tensile strength, etc. are not limited thereto.
In the process of predicting and counting reserves of target features (such as coal reserves, gas volumes, water contents and the like) in a coal mine, the embodiment of the disclosure can firstly obtain drilling data (including drilling start coordinates, end coordinates, formation feature data and the like) of a plurality of drilling holes in a target area from a database, that is, the plurality of drilling holes surround to obtain the target area, wherein the range and the shape of the target area can be flexibly set according to practical application scenes, and the target area in the embodiment of the disclosure can be rectangular, for example, without limitation.
S102: and constructing a target three-dimensional model of the target area according to the starting coordinate and the ending coordinate.
After the drilling data of the multiple drilling holes in the target area are obtained, a three-dimensional model can be further constructed for the target area according to the starting coordinate and the ending coordinate of each drilling hole, and the three-dimensional model can be called as a target three-dimensional model. In some embodiments, the start coordinates and the end coordinates of the multiple boreholes are located on the same horizontal plane, and the target area is rectangular, so that the three-dimensional model of the embodiment is a cuboid structure.
In some embodiments, in the process of constructing the target three-dimensional model of the target area, interpolation operations can be performed on the start coordinates and the end coordinates respectively to determine a plurality of model surfaces, wherein the plurality of model surfaces comprise a model top plate (top surface), a model bottom plate (bottom surface) and model side plates (i.e. four side surfaces); further, a plurality of mold surfaces will be combined to construct the target three-dimensional model.
In some embodiments, in the case that multiple rock layers exist in the target area, in order to accurately calculate the reserves of the target features, the embodiment may first construct a three-dimensional model including multiple rock layers according to the start coordinates and the end coordinates, where the three-dimensional model is called a candidate three-dimensional model; further, the present embodiment processes the model height of the candidate three-dimensional model based on the formation characteristic data (e.g., each formation thickness) to obtain a target three-dimensional model of the target formation, that is, removes other formations in the candidate three-dimensional model to obtain a target three-dimensional model of the formation (e.g., coal seam) related to the target characteristic. In practical application, the bottom surface of the candidate three-dimensional model (or the bottom surface without elevation is re-established) is taken as the bottom surface of the target three-dimensional model, and new interpolation is performed according to the spatial information (x, y) of each drilled hole in the coal seam as the elevation, so as to determine the top plate, thereby obtaining the target three-dimensional model. It will be appreciated that the spatial information (x, y) of the borehole at each layer may be determined based on the start coordinates, end coordinates of the borehole, and thickness information for each formation.
S103: interpolation is performed on the target three-dimensional model based on target feature data in the formation feature data to calculate an average amount of the target feature.
The target features are features of interest to the user in the prediction and statistics process, and the target features are as follows: coal, gas, water, etc., while data related to a target feature may be referred to as target feature data, namely: specific values of the target characteristics, for example, the volume weight of coal, the gas amount, the water content, etc., are not limited thereto. While embodiments of the present disclosure may determine target feature data from the layer feature data. It will be appreciated that the analysis of the partial target feature data in the borehole during the sampling process, in this embodiment, requires determining the average amount of the target region based on the partial target feature data during the prediction and statistics process.
In view of this, embodiments of the present disclosure may interpolate a target three-dimensional model from target feature data to calculate an average amount of target features within the target three-dimensional model. For example, where the target characteristic is coal, the average amount is the average bulk weight of the coal within the target three-dimensional model; for another example, when the target feature is gas, the average amount is the average gas content in the target three-dimensional model, which is not limited.
Some embodiments interpolate the target three-dimensional model, for example using a kriging interpolation algorithm, to calculate an average amount of target features, for example: average volume weight of coal.
S104: and calculating the reserve of the target feature in the target area according to the average volume and the volume of the target three-dimensional model.
Specifically, the present embodiment may predefine a calculation formula, for example, calculate a coal reserve, based on the target characteristics, where the calculation formula is set as follows: coal reserves = volume x average volumetric weight of coal, then this example, after determining the average volumetric weight of coal (average volume) and the volume of the target three-dimensional model (bottom area times side height), can calculate the reserves of coal in the target area from the average volume and volume.
It will be appreciated that different calculation formulas may be set for different target features to calculate the reserves of the target feature in the target area, which is not limited.
In this embodiment, by acquiring drilling data of a plurality of holes in a target area, where the drilling data includes a drilling start coordinate, a drilling end coordinate, and formation feature data, constructing a target three-dimensional model of the target area according to the start coordinate and the drilling end coordinate, performing interpolation operation on the target three-dimensional model based on the target feature data in the formation feature data to calculate an average amount of target features, and calculating reserves of the target features in the target area according to the average amount and a volume of the target three-dimensional model, the accuracy of reserve calculation can be improved, and the reserves can be calculated in real time according to the drilling data as exploitation proceeds, thereby improving timeliness of calculation.
FIG. 2 is a flow chart of a method for formation reserve calculation based on a three-dimensional model and interpolation, as shown in FIG. 2, according to another embodiment of the present disclosure, the method comprising:
s201: drilling data of a plurality of boreholes in a target area is obtained, wherein the drilling data comprises drilling start coordinates, drilling end coordinates and formation characteristic data.
S202: and constructing a target three-dimensional model of the target area according to the starting coordinate and the ending coordinate.
The specific description of S201 to S202 is referred to the above embodiments, and is not repeated here.
S203: and determining a plurality of two-dimensional planes projected in the vertical direction of the three-dimensional model of the target according to the contour principle.
In practical applications, there may be cases where the start coordinates or the end coordinates of each borehole are not located on the same horizontal plane, that is: the different regions of the target three-dimensional model are of different heights, in which case it may be difficult to calculate the volume of the target three-dimensional model.
In view of this, the embodiments of the present disclosure may determine a plurality of two-dimensional planes projected in a vertical direction of the target three-dimensional model according to a contour principle, that is, divide a top surface or a bottom surface of the target three-dimensional model into a plurality of two-dimensional planes, each of which is located at the same contour.
S204: and determining a plurality of three-dimensional sub-models included in the target three-dimensional model according to the plurality of two-dimensional planes.
That is, the target three-dimensional model is divided into a plurality of three-dimensional sub-models having different heights according to a plurality of two-dimensional planes.
S205: and respectively carrying out interpolation operation on each three-dimensional sub-model based on target feature data in the rock stratum feature data so as to calculate the average quantity of the target feature in each three-dimensional sub-model.
That is, interpolation operation is performed on each three-dimensional sub-model based on the target feature data, and the average amount of the target feature in each three-dimensional sub-model is calculated.
In addition, the present embodiment calculates the volume of each three-dimensional sub-model based on the bottom area (the area of the two-dimensional plane) and the height of each three-dimensional sub-model.
S206: and calculating the reserves of the target features in each three-dimensional sub-model according to the average volume and the average volume of each three-dimensional sub-model.
That is, reserves of the target feature (coal reserves) within each three-dimensional sub-model are calculated separately.
S207: and accumulating reserves in the plurality of three-dimensional submodels to obtain reserves of the target features in the target area.
Further, the reserves of coal in the plurality of three-dimensional submodels are accumulated to obtain reserves of target features in the target area, such as the reserves of coal.
In this embodiment, by acquiring drilling data of a plurality of holes in a target area, where the drilling data includes a drilling start coordinate, a drilling end coordinate, and formation feature data, constructing a target three-dimensional model of the target area according to the start coordinate and the drilling end coordinate, performing interpolation operation on the target three-dimensional model based on the target feature data in the formation feature data to calculate an average amount of target features, and calculating reserves of the target features in the target area according to the average amount and a volume of the target three-dimensional model, the accuracy of reserve calculation can be improved, and the reserves can be calculated in real time according to the drilling data as exploitation proceeds, thereby improving timeliness of calculation. In addition, the embodiment can divide the target three-dimensional model into a plurality of sub-models according to the contour lines, so that the reserves of the target features can be accurately calculated under the condition that the heights of the target three-dimensional models are inconsistent.
In a specific example, fig. 3 is a schematic diagram of a three-dimensional model and interpolation-based formation reserves overall calculation flow, according to an embodiment of the present disclosure, and as shown in fig. 3, the three-dimensional model and interpolation-based formation reserves overall calculation flow includes the following steps:
step 1: taking out the rock, coal seam roof and floor coordinates required to be calculated by a user from a drilling database, interpolating, establishing models of the rock, coal seam roof and floor surfaces, and further combining to form a space body of the coal and the rock stratum;
step 2: according to the model in step 1, the user is divided into a plurality of different two-dimensional plane area models (i.e., the two-dimensional planes of the above embodiments);
step 3: formulas and areas of reserves requiring statistics, such as calculating pure coal reserves, the volume of the formation requiring pure coal multiplied by the volume weight, are set.
Step 4: according to the model in the step 1 and the formula in the step 3, a user can establish a new bottom plate without elevation according to the bottom plate of the coal and the rock stratum or according to the original bottom plate, and then perform new interpolation according to the spatial information (x, y) of each drilled hole on the horizon and related characteristic information (if pure coal reserves are calculated, pure coal thickness characteristic information is needed to be used, and in other cases, characteristics such as gas content, water content and the like or a calculation result of a calculation formula can also be used) as the elevation of the point to obtain a new top plate, and then the new bottom plate and the top plate enclose a new model body as a statistical basis, namely a calculation model.
Step 5: calculating the average number (i.e., the average amount) in the specified range according to the partial parameters in the formula to calculate the statistical value (loop body in fig. 3) in the specified range:
(1) And (3) forming a new three-dimensional surface according to the two-dimensional plane area model in the step (2) projected to the coal and rock stratum bottom plate in the step (3), and taking the area as the bottom area of the area.
(2) And (3) forming a new three-dimensional body by projecting the two-dimensional plane area model in the step (2) to the coal and rock stratum body in the step (3), and taking the volume of the new three-dimensional body as the volume of the area.
(3) Averaging the formula part parameters in the region generated in step 2 according to the formula in step 3 (e.g. calculating pure coal reserves, calculating the average in this region and bringing the average into the formula to obtain the result)
(4) The bottom surface area and volume obtained above are taken into the formula set in step 3 as the average value of each parameter.
Fig. 4 is a schematic diagram of a reserve calculation device based on a three-dimensional model and interpolation provided in accordance with another embodiment of the present disclosure. As shown in fig. 4, the reserve calculating device 40 based on the three-dimensional model and interpolation includes:
an acquisition module 401, configured to acquire borehole data of a plurality of boreholes in a target area, where the borehole data includes a borehole start coordinate, a borehole end coordinate, and formation characteristic data;
a construction module 402, configured to construct a target three-dimensional model of the target region according to the start coordinate and the end coordinate;
a first calculation module 403, configured to perform interpolation operation on the target three-dimensional model based on target feature data in the formation feature data to calculate an average amount of the target feature; and
a second calculation module 404, configured to calculate a reserve of the target feature in the target area according to the average volume and the volume of the target three-dimensional model.
Some embodiments, the construction module 402 is specifically configured to: performing interpolation operation on the initial coordinates and the final coordinates to determine a plurality of model surfaces, wherein the model surfaces comprise a model top plate, a model bottom plate and model side plates; and combining the plurality of molding surfaces to construct the target three-dimensional model.
Some embodiments, the construction module 402 is specifically configured to: constructing a candidate three-dimensional model according to the initial coordinate and the termination coordinate; and processing the model height of the candidate three-dimensional model based on the rock stratum characteristic data to obtain a target three-dimensional model of the target rock stratum.
In some embodiments, the first computing module 403 is specifically configured to: determining a plurality of two-dimensional planes projected in the vertical direction of the three-dimensional model of the target according to the contour principle; determining a plurality of three-dimensional sub-models included in the target three-dimensional model according to the plurality of two-dimensional planes; respectively carrying out interpolation operation on each three-dimensional sub-model to calculate the average quantity of the target feature in each three-dimensional sub-model;
and, the second calculation module 404 is specifically configured to:
calculating reserves of the target features in each three-dimensional sub-model according to the average quantity and the volume of each three-dimensional sub-model; and accumulating reserves in the plurality of three-dimensional submodels to obtain reserves of the target features in the target area.
In some embodiments, the first computing module 403 is specifically configured to: and carrying out interpolation operation on the target three-dimensional model by using a Kriging interpolation algorithm to calculate the average quantity of the target characteristics.
In this embodiment, by acquiring drilling data of a plurality of holes in a target area, where the drilling data includes a drilling start coordinate, a drilling end coordinate, and formation feature data, constructing a target three-dimensional model of the target area according to the start coordinate and the drilling end coordinate, performing interpolation operation on the target three-dimensional model based on the target feature data in the formation feature data to calculate an average amount of target features, and calculating reserves of the target features in the target area according to the average amount and a volume of the target three-dimensional model, the accuracy of reserve calculation can be improved, and the reserves can be calculated in real time according to the drilling data as exploitation proceeds, thereby improving timeliness of calculation.
According to embodiments of the present disclosure, the present disclosure also provides a computer device, a readable storage medium and a computer program product.
To achieve the above embodiments, the present disclosure also proposes a computer program product which, when executed by an instruction processor in the computer program product, performs a formation reservoir calculation method based on a three-dimensional model and interpolation as proposed in the previous embodiments of the present disclosure.
Fig. 5 illustrates a block diagram of an exemplary computer device suitable for use in implementing embodiments of the present disclosure. The computer device 12 shown in fig. 5 is merely an example and should not be construed as limiting the functionality and scope of use of the disclosed embodiments.
As shown in FIG. 5, the computer device 12 is in the form of a general purpose computing device. Components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, a bus 18 that connects the various system components, including the system memory 28 and the processing units 16.
Computer device 12 typically includes a variety of computer system readable media. Such media can be any available media that is accessible by computer device 12 and includes both volatile and nonvolatile media, removable and non-removable media.
Although not shown in fig. 5, a magnetic disk drive for reading from and writing to a removable non-volatile magnetic disk (e.g., a "floppy disk"), and an optical disk drive for reading from or writing to a removable non-volatile optical disk (e.g., a compact disk read only memory (Compact Disc Read Only Memory; hereinafter CD-ROM), digital versatile read only optical disk (Digital Video Disc Read Only Memory; hereinafter DVD-ROM), or other optical media) may be provided. In such cases, each drive may be coupled to bus 18 through one or more data medium interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules configured to carry out the functions of the various embodiments of the disclosure.
A program/utility 40 having a set (at least one) of program modules 42 may be stored in, for example, memory 28, such program modules 42 including, but not limited to, an operating system, one or more application programs, other program modules, and program data, each or some combination of which may include an implementation of a network environment. Program modules 42 generally perform the functions and/or methods in the embodiments described in this disclosure.
The computer device 12 may also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), one or more devices that enable a user to interact with the computer device 12, and/or any devices (e.g., network card, modem, etc.) that enable the computer device 12 to communicate with one or more other computing devices. Such communication may occur through an input/output (I/O) interface 22. Moreover, the computer device 12 may also communicate with one or more networks such as a local area network (Local Area Network; hereinafter LAN), a wide area network (Wide Area Network; hereinafter WAN) and/or a public network such as the Internet via the network adapter 20. As shown, network adapter 20 communicates with other modules of computer device 12 via bus 18. It should be appreciated that although not shown, other hardware and/or software modules may be used in connection with computer device 12, including, but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, data backup storage systems, and the like.
The processing unit 16 executes various functional applications, such as implementing the formation reserves calculation method based on the three-dimensional model and interpolation mentioned in the foregoing embodiment, by running a program stored in the system memory 28.
Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This disclosure is intended to cover any adaptations, uses, or adaptations of the disclosure following the general principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.
It is to be understood that the present disclosure is not limited to the precise arrangements and instrumentalities shown in the drawings, and that various modifications and changes may be effected without departing from the scope thereof. The scope of the present disclosure is limited only by the appended claims.
It should be noted that in the description of the present disclosure, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present disclosure, unless otherwise indicated, the meaning of "a plurality" is two or more.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and further implementations are included within the scope of the preferred embodiment of the present disclosure in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present disclosure.
It should be understood that portions of the present disclosure may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.
Furthermore, each functional unit in the embodiments of the present disclosure may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present disclosure. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Although embodiments of the present disclosure have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the present disclosure, and that variations, modifications, alternatives, and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the present disclosure.
Claims (10)
1. A reserve calculation method based on a three-dimensional model and interpolation, comprising:
acquiring drilling data of a plurality of drilling holes in a target area, wherein the drilling data comprise drilling start coordinates, ending coordinates and formation characteristic data;
constructing a target three-dimensional model of the target area according to the initial coordinate and the termination coordinate;
performing interpolation operation on the target three-dimensional model based on target feature data in the rock stratum feature data to calculate average quantity of target features; and
and calculating the reserve of the target feature in the target area according to the average volume and the volume of the target three-dimensional model.
2. The method of claim 1, wherein constructing a target three-dimensional model of the target region from the start coordinates and the end coordinates comprises:
performing interpolation operation on the initial coordinates and the final coordinates to determine a plurality of mold surfaces, wherein the mold surfaces comprise a mold top plate, a mold bottom plate and mold side plates; and
combining the plurality of molding surfaces to construct the target three-dimensional model.
3. The method of claim 1, wherein constructing a target three-dimensional model of the target region from the start coordinates and the end coordinates comprises:
constructing a candidate three-dimensional model according to the initial coordinate and the termination coordinate; and
and processing the model height of the candidate three-dimensional model based on the rock stratum characteristic data to obtain a target three-dimensional model of the target rock stratum.
4. The method of claim 1, wherein interpolating the target three-dimensional model to calculate an average amount of target features comprises:
determining a plurality of two-dimensional planes projected in the vertical direction of the target three-dimensional model according to a contour principle;
determining a plurality of three-dimensional sub-models included in the target three-dimensional model according to the plurality of two-dimensional planes;
respectively carrying out interpolation operation on each three-dimensional sub-model to calculate the average quantity of the target feature in each three-dimensional sub-model;
and calculating a reserve of the target feature in the target region based on the average and the volume of the target three-dimensional model, comprising:
calculating reserves of the target features in each three-dimensional sub-model according to the average volume and the average volume of each three-dimensional sub-model; and
and accumulating reserves in the plurality of three-dimensional submodels to obtain reserves of the target features in the target area.
5. The method of claim 1, wherein interpolating the target three-dimensional model to calculate an average amount of target features comprises:
and carrying out interpolation operation on the target three-dimensional model by using a Kriging interpolation algorithm to calculate the average quantity of the target characteristics.
6. A three-dimensional model and interpolation-based formation reserve calculation apparatus, comprising:
the system comprises an acquisition module, a calculation module and a calculation module, wherein the acquisition module is used for acquiring drilling data of a plurality of drilling holes in a target area, and the drilling data comprise drilling start coordinates, termination coordinates and formation characteristic data;
the construction module is used for constructing a target three-dimensional model of the target area according to the initial coordinate and the termination coordinate;
the first calculation module is used for carrying out interpolation operation on the target three-dimensional model based on target feature data in the rock stratum feature data so as to calculate the average quantity of target features; and
and the second calculation module is used for calculating the reserve of the target feature in the target area according to the average quantity and the volume of the target three-dimensional model.
7. The apparatus of claim 6, wherein the building block is specifically configured to:
performing interpolation operation on the initial coordinates and the final coordinates to determine a plurality of mold surfaces, wherein the mold surfaces comprise a mold top plate, a mold bottom plate and mold side plates; and
combining the plurality of molding surfaces to construct the target three-dimensional model.
8. The apparatus of claim 6, wherein the building block is specifically configured to:
constructing a candidate three-dimensional model according to the initial coordinate and the termination coordinate; and
and processing the model height of the candidate three-dimensional model based on the rock stratum characteristic data to obtain a target three-dimensional model of the target rock stratum.
9. The apparatus of claim 6, wherein the first computing module is specifically configured to:
determining a plurality of two-dimensional planes projected in the vertical direction of the target three-dimensional model according to a contour principle;
determining a plurality of three-dimensional sub-models included in the target three-dimensional model according to the plurality of two-dimensional planes;
respectively carrying out interpolation operation on each three-dimensional sub-model to calculate the average quantity of the target feature in each three-dimensional sub-model;
and, the second calculation module is specifically configured to:
calculating reserves of the target features in each three-dimensional sub-model according to the average volume and the average volume of each three-dimensional sub-model; and
and accumulating reserves in the plurality of three-dimensional submodels to obtain reserves of the target features in the target area.
10. A non-transitory computer readable storage medium storing computer instructions for causing the computer to perform the method of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310118882.0A CN116109691B (en) | 2023-02-03 | 2023-02-03 | Reserve calculating method and device based on three-dimensional model and interpolation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310118882.0A CN116109691B (en) | 2023-02-03 | 2023-02-03 | Reserve calculating method and device based on three-dimensional model and interpolation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116109691A true CN116109691A (en) | 2023-05-12 |
| CN116109691B CN116109691B (en) | 2023-11-07 |
Family
ID=86267074
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310118882.0A Active CN116109691B (en) | 2023-02-03 | 2023-02-03 | Reserve calculating method and device based on three-dimensional model and interpolation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116109691B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2713019A1 (en) * | 2009-09-03 | 2011-03-03 | Schlumberger Canada Limited | Gridless geological modeling |
| KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
| CN106529755A (en) * | 2016-08-25 | 2017-03-22 | 中国黄金集团内蒙古矿业有限公司 | Mine geological resource reserve management method |
| WO2022062519A1 (en) * | 2020-09-22 | 2022-03-31 | 上海钛米机器人股份有限公司 | Ground detection method and apparatus, device, and storage medium |
| CN114429529A (en) * | 2022-04-01 | 2022-05-03 | 国家电投集团科学技术研究院有限公司 | Three-dimensional geological model generation method and device, electronic equipment and storage medium |
| CN115205476A (en) * | 2022-06-27 | 2022-10-18 | 中铁第四勘察设计院集团有限公司 | Three-dimensional geological modeling method, device, electronic equipment and storage medium |
| CN115482269A (en) * | 2022-09-22 | 2022-12-16 | 佳都科技集团股份有限公司 | Method and device for calculating earth volume, terminal equipment and storage medium |
-
2023
- 2023-02-03 CN CN202310118882.0A patent/CN116109691B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2713019A1 (en) * | 2009-09-03 | 2011-03-03 | Schlumberger Canada Limited | Gridless geological modeling |
| KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
| CN106529755A (en) * | 2016-08-25 | 2017-03-22 | 中国黄金集团内蒙古矿业有限公司 | Mine geological resource reserve management method |
| WO2022062519A1 (en) * | 2020-09-22 | 2022-03-31 | 上海钛米机器人股份有限公司 | Ground detection method and apparatus, device, and storage medium |
| CN114429529A (en) * | 2022-04-01 | 2022-05-03 | 国家电投集团科学技术研究院有限公司 | Three-dimensional geological model generation method and device, electronic equipment and storage medium |
| CN115205476A (en) * | 2022-06-27 | 2022-10-18 | 中铁第四勘察设计院集团有限公司 | Three-dimensional geological modeling method, device, electronic equipment and storage medium |
| CN115482269A (en) * | 2022-09-22 | 2022-12-16 | 佳都科技集团股份有限公司 | Method and device for calculating earth volume, terminal equipment and storage medium |
Non-Patent Citations (3)
| Title |
|---|
| 刘晓宁;罗金辉;马晓理;董志博;: "基于3DMine软件的煤炭资源数据建模与分析", 煤田地质与勘探, no. 02 * |
| 李春生;李兰勇;邓连利;: "基于钻孔数据的矿区地质模型构建", 西部探矿工程, no. 07 * |
| 陈东越;陈建平;陈三明;李长平;: "辽东白云金矿地质体三维模型的构建与储量估算", 桂林理工大学学报, no. 01 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116109691B (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2375412C (en) | An improved simulation method and apparatus | |
| US8589135B2 (en) | Systems, computer implemented methods, and computer readable program products to compute approximate well drainage pressure for a reservoir simulator | |
| US9026420B2 (en) | Method and system of updating a geocellular model | |
| CN111079217B (en) | BIM-based geotechnical engineering comprehensive investigation information interpretation method and system | |
| CN104884974A (en) | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks | |
| CN111175818A (en) | CO2Method for judging position of gas drive front edge and method and device for training model of gas drive front edge | |
| CA2691192C (en) | Correction of velocity cubes for seismic depth modeling | |
| CN111414657B (en) | Method for determining coal seam bifurcation position and establishing refined three-dimensional solid model | |
| CN114898057A (en) | Geological profile three-dimensional space information extraction method and device and storage medium | |
| CN115131486A (en) | Engineering exploration data acquisition system and method | |
| CN116109691B (en) | Reserve calculating method and device based on three-dimensional model and interpolation | |
| CN112147684A (en) | Method and device for characterizing activity intensity of co-deposition fault | |
| US7788069B2 (en) | System and method for generating grading transitions | |
| CN111413741A (en) | Sandstone-type uranium ore resource amount calculation method and device | |
| CN115880449A (en) | Method and system for establishing three-dimensional geological model based on tunnel advanced drilling data | |
| CN113742962B (en) | Shale reservoir attribute three-dimensional modeling method based on horizontal well | |
| CN112347707B (en) | Method and device for characterizing flow heterogeneity of oil reservoir | |
| CN103632396A (en) | Method for modeling complex geology strata sequence | |
| CN111765870A (en) | Earth volume calculation method based on oblique photography technology and regional accumulation | |
| CN116305451B (en) | Continuous-discontinuous geological model building method and device | |
| EP3918381B1 (en) | Hydrocarbon flow simulation | |
| US11808148B2 (en) | Systems and methods for back-allocation of oil produced by waterflooding | |
| CN114114390B (en) | Oil and gas exploration method, device, equipment and storage medium | |
| CN108537470B (en) | Fitting method of ionic rare earth ore body | |
| US20230043782A1 (en) | Layering For Geomodeling |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |