Three-dimensional stratum anti-drilling attribute volume modeling method and device, medium and equipment

Abstract

The invention relates to a three-dimensional stratum anti-drilling attribute volume modeling method and device, a medium and equipment. The three-dimensional stratum anti-drilling attribute volume modeling method comprises the following steps that different types of data are obtained; sampling different types of data into the stratigraphic grid model; according to data in the stratum grid model, a plane constraint model of the stratum anti-drilling attribute is constructed; according to the data in the stratum grid model, constructing a longitudinal constraint model of the stratum anti-drilling attribute; combining the plane constraint model of the formation anti-drilling attribute with the longitudinal constraint model of the formation anti-drilling attribute to obtain a three-dimensional constraint model of the formation anti-drilling attribute; and constructing based on the three-dimensional constraint model of the formation drilling resistance attribute and different types of data to obtain a formation three-dimensional drilling resistance attribute volume model. According to the method, effective fusion of data of different dimensions is realized, professional technical barriers among different types of data are solved, and convenient and rapid visual support is provided for deep understanding of drilling resistance characteristics of regional strata, implementation of field drilling operation and the like.

Description

Technical Field

[0001] This invention relates to the field of oil and gas field exploration technology, and more specifically, to a method, apparatus, medium, and equipment for modeling three-dimensional formation drill resistance properties. Background Technology

[0002] Currently, in the process of oil and gas field exploration and development, data from different dimensions (such as curve data, seismic data, phase diagram data, etc.) are independent of each other and cannot be effectively integrated. There are professional technical barriers between different types of data, and they cannot be used together, which makes it impossible to effectively carry out work such as gaining in-depth understanding of the drilling resistance characteristics of regional strata and implementing on-site drilling operations. Summary of the Invention

[0003] The technical problem to be solved by the present invention is to provide a method, apparatus, medium and equipment for modeling three-dimensional formation drill resistance properties, addressing the problems existing in the prior art.

[0004] The technical solution adopted by this invention to solve its technical problem is: a three-dimensional formation drill resistance property modeling method, comprising the following steps:

[0005] Acquire different types of data; these different types of data include: curve data, seismic data, and phase diagram data;

[0006] The different types of data are sampled into the formation grid model;

[0007] A planar constraint model of formation drill resistance properties is constructed based on the data in the formation grid model.

[0008] A vertical constraint model of formation drill resistance properties is constructed based on the data in the formation grid model.

[0009] By combining the planar constraint model of the formation's drill resistance property and the vertical constraint model of the formation's drill resistance property, a three-dimensional constraint model of the formation's drill resistance property is obtained.

[0010] A three-dimensional constraint model of the formation's drill resistance properties and the different types of data are used to construct a three-dimensional drill resistance property volume model of the formation.

[0011] In the three-dimensional formation drill resistance property modeling method of the present invention, the different types of data include: well logging datasets, seismic datasets, and phase diagram datasets;

[0012] The well logging dataset includes: natural gamma, sonic transit time, density, resistivity, clay content, and lithology;

[0013] The earthquake dataset includes: two-dimensional and three-dimensional segy seismic data, inverted Poisson's ratio, density, and clay content;

[0014] The phase map dataset includes: layered depositional phase plane maps.

[0015] In the three-dimensional formation drill resistance property modeling method of the present invention, the step of sampling the different types of data into the formation mesh model includes:

[0016] Based on the characteristics of the curve-type data, a corresponding data sampling method is selected for sampling, and the sampled curve-type data is sampled into the stratigraphic grid model;

[0017] The earthquake-type data is directly sampled into the stratigraphic grid model;

[0018] The phase diagram data is sampled into the formation grid model by using a phase boundary regional assignment phase code method.

[0019] In the three-dimensional formation drill resistance property modeling method of the present invention, the data sampling method includes: the maximum value method, the arithmetic mean method, and the harmonic mean method.

[0020] In the three-dimensional formation drill resistance property volume modeling method of the present invention, the step of constructing a planar constraint model of formation drill resistance properties based on the data in the formation mesh model includes:

[0021] Based on the data in the stratigraphic grid model, determine the correlation between the seismic interpretation data and the phase diagram at the corresponding locations of the curve-type data;

[0022] Based on the curve data, the seismic interpretation data, and the correlation of the phase diagram, a planar constraint model for the formation's drill resistance properties is established.

[0023] In the three-dimensional formation drill resistance property modeling method of the present invention, the step of constructing a vertical constraint model of formation drill resistance properties based on the data in the formation mesh model includes:

[0024] The vertical trend of formation drill resistance properties is determined based on the curve-type data in the formation grid model.

[0025] Based on the vertical trend and related curve data, a vertical constraint model for the formation's drill resistance properties is established.

[0026] In the three-dimensional formation drill resistance property volume modeling method described in this invention, the construction of both the planar constraint model and the vertical constraint model of the formation drill resistance property requires the following operations:

[0027] Calculate the range of the variogram;

[0028] The model trend is determined based on the range of the variogram.

[0029] Construct the planar constraint model or the longitudinal constraint model based on the model trend.

[0030] The present invention also provides a three-dimensional formation drill resistance property modeling device, comprising:

[0031] The data acquisition unit is used to acquire different types of data, including: curve data, seismic data, and phase diagram data.

[0032] A data sampling unit is used to sample the different types of data into the formation grid model;

[0033] A planar model construction unit is used to construct a planar constraint model of the formation's drill resistance properties based on the data in the formation grid model.

[0034] The vertical model building unit is used to build a vertical constraint model of the formation's drill resistance properties based on the data in the formation grid model.

[0035] A three-dimensional model construction unit is used to combine the planar constraint model of the formation's drill resistance properties and the vertical constraint model of the formation's drill resistance properties to obtain a three-dimensional constraint model of the formation's drill resistance properties.

[0036] A three-dimensional formation volume construction unit is used to construct a three-dimensional constrained model based on the formation's drill resistance properties and the different types of data to obtain a three-dimensional drill resistance property volume model of the formation.

[0037] The present invention also provides a storage medium storing a computer program adapted for loading by a processor to perform the steps of the three-dimensional formation drill resistance property volume modeling method as described above.

[0038] The present invention also provides an electronic device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the steps of the three-dimensional formation drill resistance property modeling method as described above by calling the computer program stored in the memory.

[0039] The method, apparatus, medium, and equipment for modeling three-dimensional formation drill resistance properties according to the present invention have the following beneficial effects: The method includes the following steps: acquiring different types of data; sampling different types of data into a formation mesh model; constructing a planar constraint model of formation drill resistance properties based on the data in the formation mesh model; constructing a vertical constraint model of formation drill resistance properties based on the data in the formation mesh model; combining the planar constraint model and the vertical constraint model of formation drill resistance properties to obtain a three-dimensional constraint model of formation drill resistance properties; and constructing a three-dimensional formation drill resistance property volume model based on the three-dimensional constraint model of formation drill resistance properties and different types of data. This invention achieves effective fusion of data from different dimensions, overcomes the technical barriers between different types of data, and provides convenient and rapid visualization support for in-depth understanding of regional formation drill resistance characteristics and on-site drilling operations. Attached Figure Description

[0040] The present invention will be further described below with reference to the accompanying drawings and embodiments. In the accompanying drawings:

[0041] Figure 1 This is a flowchart illustrating the three-dimensional formation drill resistance property modeling method provided in this embodiment of the invention;

[0042] Figure 2 This is a schematic diagram of the logic block of the three-dimensional formation anti-drilling attribute body modeling device provided in the embodiment of the present invention. Detailed Implementation

[0043] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0044] This invention provides a three-dimensional formation drill resistance attribute volume modeling method, which is an integrated formation drill resistance attribute volume modeling method based on the joint construction of deterministic modeling technology and stochastic modeling technology. After standardizing and managing various types of formation drill resistance attribute data, it samples them into the formation grid model in different ways. Different modeling techniques are selected according to the data characteristics and attribute modeling control process to achieve effective integration of data from different dimensions. This solves the professional technical barriers between different types of data and provides convenient and fast visualization support for in-depth understanding of regional formation drill resistance characteristics and implementation of field drilling operations.

[0045] refer to Figure 1 In a preferred embodiment, the three-dimensional formation drill resistance property modeling method includes the following steps:

[0046] Step S101: Obtain data of different types.

[0047] In this embodiment of the invention, different types of data include: curve data, seismic data, and phase diagram data.

[0048] In some embodiments, different types of data include: well logging datasets, seismic datasets, and facies map datasets; well logging datasets include: natural gamma, sonic transit time, density, resistivity, shale content, and lithology; seismic datasets include: two-dimensional and three-dimensional Segy seismic data, inverted Poisson's ratio, density, and shale content; facies map datasets include: layered sedimentary facies plane maps.

[0049] Specifically, well logging datasets are data based on field geological well logging data, including but not limited to: natural gamma ray, sonic transit time, density, resistivity, clay content, lithology, etc. Each well constitutes a dataset.

[0050] The seismic dataset mainly consists of two-dimensional and three-dimensional seismic data measured during the exploration and development phases, as well as formation drill resistance properties such as Poisson's ratio, density, and clay content obtained through inversion calculations.

[0051] The phase diagram dataset mainly consists of layered sedimentary phase plane diagrams.

[0052] Step S102: Sample different types of data into the formation grid model.

[0053] In some embodiments, sampling different types of data into the stratigraphic grid model includes: selecting the corresponding data sampling method based on the data characteristics of curve-type data, and sampling the sampled curve-type data into the stratigraphic grid model; directly assigning seismic data to the stratigraphic grid model; and using the phase boundary regional assignment phase code method to sample phase diagram data into the stratigraphic grid model.

[0054] Specifically, for curve-type data, the optimal data sampling method can be selected based on its characteristics. These methods include, but are not limited to, the maximum value method, the arithmetic mean method, and the harmonic mean method. For example, the maximum value method is generally chosen for averaging discrete data types such as sedimentary bodies or lithofacies; the arithmetic mean method is generally chosen for porosity and water saturation; and the harmonic mean method is generally chosen for permeability.

[0055] Step S103: Construct a planar constraint model of formation drill resistance properties based on the data in the formation grid model.

[0056] In some embodiments, constructing a planar constraint model of formation drill resistance properties based on data in the formation grid model includes: determining the correlation between seismic interpretation data and phase diagrams at locations corresponding to curve-type data based on data in the formation grid model; and establishing a planar constraint model of formation drill resistance properties based on the correlation between curve-type data, seismic interpretation data, and phase diagrams.

[0057] Specifically, based on the relationships between curve-type data in the stratigraphic grid model and the corresponding seismic interpretation data and phase diagrams, a plane-constrained model of the formation's drill resistance properties is established. The seismic interpretation data mainly includes post-stack amplitude volumes and inverted stratigraphic property volumes (such as density volumes, wave impedance volumes, or some inverted drill resistance property volumes). After determining the seismic interpretation data at the corresponding locations using the curve-type data, a planar trend surface is formed using the correlation between the seismic interpretation data and the phase diagram. Based on this planar trend surface, a plane interpolation algorithm is used to constrain the interpolation trend, ultimately generating a plane-constrained model of the formation's drill resistance properties.

[0058] The construction of the plane-constrained model for formation drill resistance properties first requires calculating the range of the variogram (in the plane-constrained model construction, this range is the plane range). Then, the model trend is determined based on the range of the variogram, and finally, the plane-constrained model for formation drill resistance properties is constructed based on this trend. Specifically, the plane-constrained model for formation drill resistance properties mainly determines the trend model by analyzing the range of the plane variogram of the input data (curve data, seismic data, phase diagrams, and other sampled data). The core calculation formula for the variogram is:

[0059]

[0060] (1) In the formula, r * (h) represents the planar range when Z(x) and Z(x+h) are two points on the same plane; Z(x) and Z(x+h) are the variances between the two points; N(h) is the total number of points; and h is the increment. The variance between the two points is obtained by subtracting the square of the average value of the regional sampled data from the variances of the curve data, seismic data, phase diagrams, etc., from the variances of the sampled data in adjacent grids.

[0061] Step S104: Construct a vertical constraint model of formation drill resistance properties based on the data in the formation grid model.

[0062] In some embodiments, constructing a vertical constraint model of formation drill resistance properties based on data in the formation grid model includes: determining the vertical trend of formation drill resistance properties based on curve-like data in the formation grid model; and establishing a vertical constraint model of formation drill resistance properties based on the vertical trend and related curve-like data.

[0063] The construction of the vertical constraint model for formation drill resistance properties is as follows: The range of the variogram is calculated (in the construction of the vertical constraint model, this range is the vertical range); the model trend is determined based on the range of the variogram; and a planar constraint model or a vertical constraint model is constructed based on the model trend. Specifically, the vertical trend of formation drill resistance properties is determined based on the sampled relevant curve data (such as density, porosity, permeability, water saturation, sonic transit time, drill resistance, etc.), and then the vertical constraint model of formation drill resistance properties is constructed based on the relationship between the acquired relevant curve data and the vertical trend. For example, the trend of drill resistance property curves already coarsened into the formation grid model along the vertical depth can be statistically analyzed (for example, the variation law of drill resistance curve values ​​at different depths in different wells can be statistically analyzed to guide the vertical interpolation trend of the interpolation algorithm), the vertical attribute variation trend curves at different locations are regressed, and then multiple vertical trend curves are bound to the locations to form a vertical constraint model of formation drill resistance properties.

[0064] The plane-constrained model of formation drill resistance properties mainly determines the trend model by analyzing the range of the plane variation function of the input data (sampling data such as curve data, seismic data, and phase diagrams). The core calculation formula for the variation function is:

[0065]

[0066] (2) In the formula, r * (h) represents the longitudinal range when Z(x) and Z(x+h) are two points in the longitudinal direction; Z(x) and Z(x+h) are the variances between the two points; N(h) is the total number of points; and h is the increment.

[0067] Step S105: Combine the planar constraint model of formation drill resistance properties and the vertical constraint model of formation drill resistance properties to obtain a three-dimensional constraint model of formation drill resistance properties.

[0068] In some embodiments, the combination of the planar constraint model and the vertical constraint model of formation drill resistance attributes involves simultaneously inputting planar and vertical ranges during formation drill resistance attribute modeling to obtain a three-dimensional constraint model of formation drill resistance attributes that shares common characteristics with curve data, formation, and phase diagram. The combination process primarily involves using a spatial interpolation algorithm to interpolate the planar constraint model and the vertical constraint model of formation drill resistance attributes to form a three-dimensional constraint model (a three-dimensional trend volume) of formation drill resistance attributes.

[0069] Step S106: Construct a three-dimensional constraint model based on the formation's anti-drilling properties and different types of data to obtain a three-dimensional anti-drilling property volume model of the formation.

[0070] Specifically, based on the acquired different types of data (curve data, seismic data, phase diagram data) (source data), and under the constraint of the three-dimensional constraint model of formation drill resistance properties, an interpolation algorithm is used to generate the corresponding three-dimensional formation drill resistance property volume model.

[0071] In this embodiment of the invention, when performing attribute modeling, the conditional joint probability model of N random variables in the sequential Gaussian simulation of the stochastic modeling method used is as follows:

[0072] F N [z1,z2,...,z N |(n)]=Prob{Z i ≤z i ,i=1,2,...,N|(n)} (3);

[0073] (3) In the formula, F N [z1,z2,....,z N [|(n)] represents the joint conditional distribution of N points; Z1, Z2, ... Z N , is a random variable; N is the total number; z N |(n) represents the vector of the nth known sample value.

[0074] The two-dimensional indicator random variable in the sequential indicator simulation is:

[0075]

[0076] (4) In the formula, I(x; z0) represents a binary random variable that can only take the value 0 or 1; {Z(x), x∈A} is the random function corresponding to the regional variable, where A represents a certain region and a threshold value of z0 is given for the regional variable.

[0077] The Kriging equations in the sequential instruction are:

[0078]

[0079] (5) In the formula, C I (xx a (x; z0) is the covariance function; (x; z0) is the Lagrange multiplier; β is the βth known sample point; λ β The Kriging weights are assigned to the indicator value of the β-th known sample point; μ(x, z0) is the indicator variation function model; and n is the number of known sample points used as conditional data.

[0080] refer to Figure 2 The present invention also provides a three-dimensional formation anti-drilling property modeling device.

[0081] like Figure 2As shown, in a preferred embodiment, the three-dimensional formation drill resistance property modeling device includes:

[0082] The data acquisition unit 201 is used to acquire different types of data, including curve data, seismic data, and phase diagram data.

[0083] Data sampling unit 202 is used to sample different types of data into the formation grid model;

[0084] Planar model building unit 203 is used to build a planar constraint model of formation drill resistance properties based on data in the formation grid model.

[0085] The vertical model building unit 204 is used to build a vertical constraint model of formation drill resistance properties based on the data in the formation grid model.

[0086] The three-dimensional model building unit 205 is used to combine the planar constraint model of the formation's drill resistance properties and the vertical constraint model of the formation's drill resistance properties to obtain a three-dimensional constraint model of the formation's drill resistance properties.

[0087] The three-dimensional formation volume construction unit 206 is used to construct a three-dimensional constrained model based on the formation's anti-drilling properties and different types of data to obtain a three-dimensional formation anti-drilling property volume model.

[0088] The above principles are illustrated with examples.

[0089] Example 1. The target area is block C35.

[0090] Formation drillability property modeling was conducted in block C35. A three-dimensional formation drillability property model was created by comprehensively applying seismic, geological, and well logging data. This model included structural, facies, rock hardness, compressive strength, drillability, porosity, permeability, and water saturation models. Verification showed that the structural model matched the seismic interpretation horizons and geological strata, and the distribution characteristics of the property models were consistent with sedimentary facies and well logging curves. Details are shown in Table 1.

[0091] Table 1

[0092]

[0093] Example 2. The target area is block J50.

[0094] Formation drillability property modeling was conducted in block J50. A three-dimensional formation drillability property model was created in block J50 by comprehensively applying seismic, geological, and well logging data. This model includes structural models, lithofacies models, Poisson's ratio models, rock hardness models, compressive strength models, drillability models, porosity models, permeability models, and sonic transit time models. Verification showed that the structural model matched the seismic interpretation horizons and geological strata, and the property model distribution at well points was consistent with the well logging data. The model was constrained by the seismic data volume between wells and showed a trend consistent with the seismic inversion data volume. Details are shown in Table 2.

[0095] Table 2

[0096]

[0097]

[0098] Example 3. The target area is block T4.

[0099] Formation drillability property modeling was conducted in Block T4. A three-dimensional formation drillability property model was created in Block T4 by comprehensively applying seismic, geological, and well logging data. This model included structural models, Poisson's ratio models, drillability models, porosity models, and permeability models. Verification showed that the structural model matched the seismic interpretation horizons and geological strata, and the property model distribution was consistent with the well point and well logging data. Details are shown in Table 3.

[0100] Table 3

[0101]

[0102] As can be seen from the examples of the three different blocks mentioned above, this invention achieves effective integration of data from different dimensions, and the established three-dimensional formation anti-drilling attribute model has high reliability. It provides convenient and quick visualization support for in-depth understanding of regional formation anti-drilling characteristics and implementation of on-site drilling operations, improves the efficiency of on-site drilling operation decision-making, significantly reduces costs, and enhances economic benefits.

[0103] Specifically, the specific operational process of the various units in the three-dimensional formation drill resistance property modeling device can be referred to the three-dimensional formation drill resistance property modeling method mentioned above, and will not be repeated here.

[0104] Furthermore, an electronic device of the present invention includes a memory and a processor; the memory is used to store a computer program; the processor is used to execute the computer program to implement the three-dimensional formation drill resistance property volume modeling method as described above. Specifically, according to embodiments of the present invention, the process described above with reference to the flowchart can be implemented as a computer software program. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a computer-readable medium, the computer program containing program code for performing the method shown in the flowchart. In such embodiments, when the computer program is downloaded, installed, and executed by an electronic device, it performs the functions defined above in the method of the embodiments of the present invention. The electronic device in the present invention can be a terminal such as a laptop, desktop computer, tablet computer, or smartphone, or it can be a server.

[0105] Furthermore, one type of storage medium of the present invention stores a computer program thereon, which, when executed by a processor, implements the three-dimensional formation drill resistance property modeling method described above. Specifically, it should be noted that the storage medium described above in the present invention can be a computer-readable signal medium or a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, device, or apparatus. In the present invention, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. The transmitted data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. The computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wires, optical fibers, RF (radio frequency), etc., or any suitable combination thereof.

[0106] The aforementioned computer-readable medium may be included in the aforementioned electronic device; or it may exist independently and not assembled into the electronic device.

[0107] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to the method section.

[0108] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0109] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0110] The above embodiments are only for illustrating the technical concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement it accordingly. They do not limit the scope of protection of the present invention. All equivalent changes and modifications made within the scope of the claims of the present invention should fall within the scope of the claims of the present invention.

Claims

1. A method for modeling three-dimensional formation drill resistance properties, characterized in that, Includes the following steps: Acquire different types of data; these different types of data include: curve data, seismic data, and phase diagram data; The different types of data are sampled into the formation grid model; A planar constraint model of formation drill resistance properties is constructed based on the data in the formation grid model. A vertical constraint model of formation drill resistance properties is constructed based on the data in the formation grid model. By combining the planar constraint model of the formation's drill resistance property and the vertical constraint model of the formation's drill resistance property, a three-dimensional constraint model of the formation's drill resistance property is obtained. A three-dimensional constraint model of the formation's drill resistance properties and the different types of data are used to construct a three-dimensional drill resistance property volume model of the formation.

2. The three-dimensional formation drill resistance property modeling method according to claim 1, characterized in that, The different types of data include: well logging datasets, seismic datasets, and phase diagram datasets; The well logging dataset includes: natural gamma, sonic transit time, density, resistivity, clay content, and lithology; The earthquake dataset includes: two-dimensional and three-dimensional segy seismic data, inverted Poisson's ratio, density, and clay content; The phase map dataset includes: layered depositional phase plane maps.

3. The three-dimensional formation drill resistance property modeling method according to claim 1, characterized in that, The step of sampling the different types of data into the formation grid model includes: Based on the characteristics of the curve-type data, a corresponding data sampling method is selected for sampling, and the sampled curve-type data is sampled into the stratigraphic grid model; The earthquake-type data is directly sampled into the stratigraphic grid model; The phase diagram data is sampled into the formation grid model by using a phase boundary regional assignment phase code method.

4. The three-dimensional formation drill resistance property modeling method according to claim 3, characterized in that, The data sampling methods include: the maximum value method, the arithmetic mean method, and the harmonic mean method.

5. The three-dimensional formation drill resistance property modeling method according to claim 1, characterized in that, The step of constructing a planar constraint model of formation drill resistance properties based on data from the formation grid model includes: Based on the data in the stratigraphic grid model, determine the correlation between the seismic interpretation data and the phase diagram at the corresponding locations of the curve-type data; Based on the curve data, the seismic interpretation data, and the correlation of the phase diagram, a planar constraint model for the formation's drill resistance properties is established.

6. The three-dimensional formation drill resistance property modeling method according to claim 1, characterized in that, The process of constructing a vertical constraint model of formation drill resistance properties based on data from the formation grid model includes: The vertical trend of formation drill resistance properties is determined based on the curve-type data in the formation grid model. Based on the vertical trend and related curve data, a vertical constraint model for the formation's drill resistance properties is established.

7. The three-dimensional formation drill resistance property modeling method according to claim 5 or 6, characterized in that, The construction of both the planar constraint model and the vertical constraint model of the formation's drill resistance properties requires the following operations: Calculate the range of the variogram; The model trend is determined based on the range of the variogram. Construct the planar constraint model or the longitudinal constraint model based on the model trend.

8. A three-dimensional formation drill resistance property modeling device, characterized in that, include: The data acquisition unit is used to acquire different types of data; The different types of data include: curve data, seismic data, and phase diagram data; A data sampling unit is used to sample the different types of data into the formation grid model; A planar model construction unit is used to construct a planar constraint model of the formation's drill resistance properties based on the data in the formation grid model. The vertical model building unit is used to build a vertical constraint model of the formation's drill resistance properties based on the data in the formation grid model. A three-dimensional model construction unit is used to combine the planar constraint model of the formation's drill resistance properties and the vertical constraint model of the formation's drill resistance properties to obtain a three-dimensional constraint model of the formation's drill resistance properties. A three-dimensional formation volume construction unit is used to construct a three-dimensional constrained model based on the formation's drill resistance properties and the different types of data to obtain a three-dimensional drill resistance property volume model of the formation.

9. A storage medium, characterized in that, The storage medium stores a computer program adapted for loading by a processor to perform the steps of the three-dimensional formation anti-drilling property volume modeling method as described in any one of claims 1 to 7.

10. An electronic device, characterized in that, It includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the steps of the three-dimensional formation drill resistance property volume modeling method as described in any one of claims 1 to 7 by calling the computer program stored in the memory.

Images