CN112233235B - Shale gas horizontal well space analysis method, device and storage medium - Google Patents
Shale gas horizontal well space analysis method, device and storage medium Download PDFInfo
- Publication number
- CN112233235B CN112233235B CN202011114911.9A CN202011114911A CN112233235B CN 112233235 B CN112233235 B CN 112233235B CN 202011114911 A CN202011114911 A CN 202011114911A CN 112233235 B CN112233235 B CN 112233235B
- Authority
- CN
- China
- Prior art keywords
- well
- index
- drilled
- drilling
- transformation index
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 230000009466 transformation Effects 0.000 claims abstract description 72
- 238000005553 drilling Methods 0.000 claims abstract description 69
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012732 spatial analysis Methods 0.000 claims abstract description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 26
- 238000013461 design Methods 0.000 claims description 26
- 238000005755 formation reaction Methods 0.000 claims description 26
- 238000013178 mathematical model Methods 0.000 claims description 23
- 230000004048 modification Effects 0.000 claims description 23
- 238000012986 modification Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000006872 improvement Effects 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 37
- 238000010586 diagram Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 208000033830 Hot Flashes Diseases 0.000 description 1
- 206010060800 Hot flush Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application provides a shale gas horizontal well space analysis method, a shale gas horizontal well space analysis device and a shale gas horizontal well space analysis storage medium, and relates to the technical field of geology. The method comprises the following steps: establishing a transformation index three-dimensional space model according to the given area and the depth range; extracting an reconstruction index slice along a planned deployment azimuth of drilling according to the reconstruction index three-dimensional space model; and analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice. In the method, the spatial position of the well drilling is evaluated by establishing a transformation index three-dimensional space model and obtaining a transformation index slice by using the model, so that the accuracy of shale gas horizontal well spatial analysis and lost circulation early warning is improved.
Description
Technical Field
The application relates to the technical field of geology, in particular to a shale gas horizontal well space analysis method, a shale gas horizontal well space analysis device and a storage medium.
Background
Along with the continuous expansion of the social demand for clean energy, the price of natural gas is continuously increased, and the understanding of people on shale gas is rapidly improved. In particular, the horizontal well and the fracturing technology are continuously advanced, and the exploration and development of shale gas by human beings are forming hot flashes.
Hydraulic fracturing is the main form of natural gas exploitation at present, and is to squeeze fracturing fluid with higher viscosity into an oil layer through a shaft by using a ground high-pressure pump. When the rate of injection of the fracturing fluid exceeds the absorption capacity of the reservoir, a high pressure builds up on the reservoir at the bottom of the well and when this pressure exceeds the fracture pressure of the reservoir rock near the bottom of the well, the reservoir will be forced apart and a fracture will develop. At this time, the fracturing fluid is continuously squeezed into the oil layer, and the cracks are continuously expanded into the oil layer.
In the prior art, through some methods, parameters such as seam length, seam width, seam height and the like of single well fracturing transformation can be only simulated, and corresponding methods for space analysis and deployment of shale gas horizontal wells are also lacking.
Disclosure of Invention
In order to solve the problems in the prior art, the application provides a shale gas horizontal well space analysis method, a shale gas horizontal well space analysis device and a storage medium.
In order to achieve the above purpose, the technical scheme adopted in the application is as follows:
the first aspect of the application provides a shale gas horizontal well space analysis method, which comprises the following steps:
establishing a transformation index three-dimensional space model according to the given area and the depth range;
extracting an reconstruction index slice along a planned deployment azimuth of drilling according to the reconstruction index three-dimensional space model;
and analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice.
Optionally, before the reconstruction index three-dimensional space model is built according to the given area and the depth range, the method further comprises:
and establishing an index reconstruction mathematical model by using the spherical wave field according to the related well-drilled space parameters.
Optionally, the analyzing to obtain the lost circulation early warning information according to the relevant spatial parameter of the drilled well and the transformation index slice includes:
calculating and obtaining the transformation index of each point in the stratum by adopting the transformation index mathematical model;
calculating an edge layer arithmetic average slice and a root mean square slice of the transformation index slices according to the transformation index of each point in the stratum;
and analyzing and acquiring well leakage early warning information according to the layer edge arithmetic average slice and the root mean square slice.
Optionally, if there are a plurality of drilled wells that produce a coincidence in the formation, the calculating the modification index for each point in the formation using the modification index mathematical model includes:
and summing the plurality of the drilled modification indexes, and calculating the modification index of each point in the acquired stratum.
Optionally, the method further comprises:
acquiring the design track of the drilling well and the related well-drilled track related to the drilling well;
discretizing the design track of the on-drilling well and the related well-drilled track to obtain the discrete data of the on-drilling well and the discrete data of the well-drilled well;
calculating the shortest distance between the design track of the well and the related well-drilled track according to the well-drilled discrete data and the well-drilled discrete data;
and acquiring well leakage early warning information through the shortest distance analysis.
A second aspect of the present application provides a shale gas horizontal well spatial analysis apparatus, comprising: the device comprises a building unit, an extraction unit and an analysis acquisition unit;
the establishing unit is used for establishing a transformation index three-dimensional space model according to the given area and the depth range;
the extraction unit is used for extracting an index reconstruction slice along the planned deployment azimuth of the well drilling according to the index reconstruction three-dimensional space model;
and the analysis acquisition unit is used for analyzing and acquiring lost circulation early warning information according to the related well-drilled space parameters and the transformation index slice.
Optionally, the establishing unit is configured to establish a reconstruction index mathematical model with a spherical wave field according to the relevant spatial parameters of the drilled well.
Optionally, the analysis obtaining unit is configured to calculate and obtain a transformation index of each point in the stratum by using the transformation index mathematical model;
calculating an edge layer arithmetic average slice and a root mean square slice of the transformation index slices according to the transformation index of each point in the stratum;
and analyzing and acquiring well leakage early warning information according to the layer edge arithmetic average slice and the root mean square slice.
Optionally, if there are a plurality of drilled formations that produce a coincidence in the formation, the analysis acquisition unit is configured to sum the plurality of drilled formations and calculate an improvement index for each point in the formation.
Optionally, the apparatus further comprises: an acquisition unit, a discrete unit, and a calculation unit;
the acquisition unit is used for acquiring the design track of the drilling well and the related drilling well track related to the drilling well;
the discrete unit is used for discretizing the design track of the drilling well and the related track which is drilled well to obtain the discrete data of the drilling well and the discrete data of the drilled well;
the calculating unit is used for calculating the shortest distance between the design track of the well drilling and the related well drilling track according to the well drilling discrete data and the well drilling discrete data;
and the analysis acquisition unit is used for acquiring the lost circulation early warning information through the shortest distance analysis.
A third aspect of the present application provides a shale gas horizontal well spatial analysis apparatus, comprising: a processor, a storage medium, and a bus, the storage medium storing machine-readable instructions executable by the processor, the processor and the storage medium communicating over the bus when the apparatus is running, the processor executing the machine-readable instructions to perform the steps of the method of the first aspect described above.
A fourth aspect of the present application provides a storage medium having stored thereon a computer program which, when executed by a processor, performs a method as provided in the first aspect.
In the shale gas horizontal well space analysis method, the shale gas horizontal well space analysis device and the storage medium, an index three-dimensional space model is established according to a given area and a depth range; extracting an reconstruction index slice along a planned deployment azimuth of drilling according to the reconstruction index three-dimensional space model; and analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice. By establishing a transformation index three-dimensional space model and utilizing the model to obtain transformation index slices, the deployment of the drilling well is evaluated, and the accuracy of shale gas horizontal well space analysis and lost circulation early warning is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered limiting the scope, and that other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic flow chart of a shale gas horizontal well space analysis method according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a simulation of a reconstruction index mathematical model provided in an embodiment of the present application;
FIG. 3 is a schematic flow chart of a method for spatial analysis of shale gas horizontal wells according to another embodiment of the present application;
FIG. 4 is a schematic flow chart of a method for spatial analysis of shale gas horizontal wells according to another embodiment of the present application;
fig. 5 is a schematic structural diagram of a shale gas horizontal well space analysis device according to an embodiment of the present disclosure;
fig. 6 is a schematic structural diagram of a shale gas horizontal well space analysis device according to another embodiment of the present disclosure;
fig. 7 is a schematic structural diagram of a shale gas horizontal well space analysis device according to an embodiment of the present application.
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments herein without making any inventive effort, are intended to be within the scope of the present application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
Furthermore, the terms "first," "second," and the like, if any, are used merely for distinguishing between descriptions and not for indicating or implying a relative importance.
It should be noted that, without conflict, features in embodiments of the present application may be combined with each other.
Aiming at shale gas reservoirs developed by adopting a hydraulic fracturing technology, the life cycle of a shale gas horizontal well is generally short, and due to the fact that a well pattern is too dense, firstly, when a new well is designed, the condition of formation transformation caused by well drilling fracturing needs to be fully considered, reasonable deployment of inter-well relations is realized, effective distribution of the well pattern is guaranteed, and the gas reservoir potential is fully exerted; secondly, the new well drilling and the existing well drilling have complex spatial relationship, mutual interference is prevented in the well drilling process, and early warning is timely provided; thirdly, due to the influence of the artificial joint formed by the natural joint and the fracturing, well leakage analysis and early warning are required to be carried out in the drilling process.
Based on the above, the embodiment of the application provides a shale gas horizontal well space analysis method, which is used for solving the technical problems. Fig. 1 is a schematic flow chart of a shale gas horizontal well space analysis method provided in an embodiment of the present application, as shown in fig. 1, the method includes:
s101, establishing a transformation index three-dimensional space model according to a given area and a depth range.
It should be noted that, in the embodiment of the present application, before the complete transformation index three-dimensional space model is established, a three-dimensional gridded data volume needs to be established in a given area, specifically, a given work area plane and a certain depth range in the plane. The three-dimensional meshed data volume is an initialization model of the transformation index three-dimensional space model, and can be subjected to superposition analysis with the three-dimensional seismic model and an attribute volume of the three-dimensional seismic model, and the plane grid parameters of the three-dimensional meshed data volume are generally required to be the same as those of the three-dimensional seismic model, namely the azimuth, the line interval and the channel interval of the three-dimensional meshed data volume are the same as those of the three-dimensional seismic model.
Illustratively, within a certain planar extent and depth, a (I, J, K) three-dimensional meshed data volume may be formed, typically in units of 1 m. Wherein I, J may represent horizontal position information of a certain position in the three-dimensional meshed data volume, and K may represent depth position information of a certain position in the three-dimensional meshed data volume.
In addition, in the three-dimensional meshed data body, the mesh division is not limited to "1m", but may be 2m, 3m and above, and specifically, the division determination may be performed with a given range interval and the accuracy required by the data, which is not limited in the embodiment of the present application.
After the shale gas horizontal well is fractured, importing parameters such as seam length, seam width, seam height and the like of the shale gas horizontal well after fracturing into a three-dimensional meshed data body so as to update the initial three-dimensional meshed data body and form a final transformation index three-dimensional space model. In the embodiment of the application, the parameters of the introduced shale gas horizontal well such as the seam length, the seam width, the seam height and the like are a certain plane range and all the relevant parameters of the well drilling in the depth.
The fracturing of the horizontal well can adopt the modes of horizontal well flow-limiting fracturing, double-clamping staged fracturing, hydraulic fracturing and the like, and is not limited by the embodiment of the application.
S102, extracting a transformation index slice along a planned deployment azimuth of drilling according to the transformation index three-dimensional space model.
After the reconstruction index three-dimensional space model is built, in order to obtain the formation reconstruction condition of the well drilling on the incompletely developed layer system in the well drilling deployment position, the reconstruction index slice can be extracted from the reconstruction index three-dimensional space model along the planned deployment position of the well drilling.
It should be noted that in the embodiment of the present application, the modification index slice may be a three-dimensional spatial spread of a certain tangential plane in the stratum. By analyzing the three-dimensional spatial spread of the tangential plane, the spatial arrangement of the drilled well in the intended deployment orientation of the drilled well and the modification of the formation around the drilled well by the drilled well can be obtained.
S103, analyzing and obtaining lost circulation early warning information according to the relevant space parameters of the drilled well and the reconstruction index slice.
In embodiments of the present application, it is desirable to find relevant wells in connection with drilling from thousands of horizontal wells drilled prior to analysis of spatial information in the well. The determination of the relevant well may be based on information about the relevant parameters that have been drilled, as well as on the extent of the drilling operation.
Specifically, the construction range of the well may be determined first, the determination of the construction range may be determined according to the conditions of the drilling equipment and the number of constructors, and the like, and the well may be determined to be the construction range within 200m of the well, and when the well head, the target a and the target B3 points of the well fall within the construction range of the well, the well is determined to be the related well related to the well in the well.
In the embodiment of the present application, the target points a and B are predetermined target points of the well track calculated according to certain data, and the final point of the well track. Mainly by means of an orienting instrument and an orienting person. The target A is a preset target point of the well track, and the target b is the final point of the well track.
In this embodiment, after determining that the well is drilled in association with the well being drilled, the lost circulation warning information in the well may be obtained by analysis according to the spatial parameters associated with the well being drilled and the transformation index slice of the planned deployment azimuth in the well.
It should be noted that, in the embodiment of the present application, the spatial parameter related to the drilled hole may be the trajectory parameter such as the length, width, height of the drilled hole.
In the shale gas horizontal well space analysis method provided by the embodiment of the application, an index three-dimensional space model is established according to a given area and a depth range; extracting an reconstruction index slice along a planned deployment azimuth of drilling according to the reconstruction index three-dimensional space model; and analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice. By establishing a transformation index three-dimensional space model and utilizing the model to obtain transformation index slices, the spatial position of the well drilling is evaluated, and the accuracy of shale gas horizontal well spatial analysis and lost circulation early warning is improved.
Optionally, before the reconstruction index three-dimensional space model is built according to the given area and the depth range, the method further comprises: and establishing an index reconstruction mathematical model by using the spherical wave field according to the related well-drilled space parameters.
In the embodiment of the application, engineering and geological parameters such as total liquid amount, total sand amount and displacement of hydraulic fracturing after drilling, formation fracture pressure and the like are applied, simulation of an artificial transformation seam can be realized, and parameters such as seam length, seam width, seam height and the like of the transformation seam along a fracturing section are calculated. Traditionally, these parameters are considered to be only a series of "points" along the horizontal segment, in fact, from the mathematical point of view, the parameters of the seam length, the seam width, the seam height and the like reflect a spherical wave field, and from the energy point of view, the hydraulic fracturing is modified with the distance from the center point to be far and near, and has a damping trend, so that a simulated seam mathematical model can be established by using the spherical wave field with the parameters of the seam length, the seam width, the seam height, the distance from the center point and the like, so as to reflect the modification condition of the single-segment fracturing to the stratum.
Fig. 2 is a schematic diagram of a simulation of a reconstruction index mathematical model provided in an embodiment of the present application. Fig. 3 is a flow chart of a shale gas horizontal well space analysis method according to another embodiment of the present application, as shown in fig. 3, the analyzing to obtain lost circulation early warning information according to the relevant spatial parameters of the drilled well and the transformation index slice includes:
s201, calculating and obtaining the transformation index of each point in the stratum by adopting a transformation index mathematical model.
In embodiments of the present application, a modification index mathematical model may be utilized to calculate a modification index for each point in the formation. As shown in FIG. 2, the half-seam height of the simulated seam in the transformation index mathematical model is set to be H 2 Half seam length L 2 Half seam width W 2 Wherein the coordinates of any point relative to the perforation center point o are (l, w, h), and the transformation index d is:
wherein l<＝L 2 ，w<＝W 2 ，h<＝H 2 。
S202, calculating an edge layer arithmetic average slice and a root mean square slice of the transformation index slice according to the transformation index of each point in the stratum.
In the embodiment of the application, the arithmetic average slice and the root mean square slice of the upper edge layer of the reconstruction index slices extracted from the planned deployment position of the drilling well are determined through a calculation method of the reconstruction index of each point in the stratum.
It should be noted that, in this embodiment, the modification index slice may also be set to a certain thickness according to the design requirement. And calculating an edge layer arithmetic average slice and a root mean square slice of the reconstruction index slice within a certain thickness of the reconstruction index slice.
And (3) performing arithmetic average calculation on the transformation index within a certain thickness of the transformation index slice by using the layer-following arithmetic average slice of the transformation index slice. The layer edge root mean square section of the transformation index section is obtained by performing root mean square calculation on the transformation index within a certain thickness of the transformation index section.
S203, analyzing and acquiring well leakage early warning information according to the edge layer arithmetic average slice and the root mean square slice.
And analyzing and acquiring well leakage early warning information of the drill to be deployed through the acquired layer-following arithmetic average slice and root mean square slice.
Specifically, in the embodiment of the application, when the values of the average slice and the root mean square slice of the along-layer arithmetic are larger, the formation reconstruction at the position is shown to be larger, and the probability of lost circulation of the horizontal well deployed at the position is larger, so that the position where the formation reconstruction is larger should be avoided as much as possible when the horizontal well is deployed in drilling.
Therefore, the deployment scheme in the well can be adjusted according to the numerical values of the edge layer arithmetic average slice and the root mean square slice, and after the deployment scheme is determined, the deployment scheme in the well needs to be adjusted in real time according to the deployment track in the well and other related parameters in the well drilling process because errors may exist in data acquisition.
Optionally, if there are a plurality of drilled wells that produce a coincidence in the formation, the calculating the modification index for each point in the formation using the modification index mathematical model includes:
and summing the plurality of the drilled modification indexes, and calculating the modification index of each point in the acquired stratum.
In this embodiment, when a plurality of manual slits affect the same grid point (i, j, k), a plurality of modification indexes of the grid point may be summed.
Illustratively, when there are N workers stitched at the same grid point (i, j, k) with modification indexes d1, d2 … dN, respectively, the modification index d (i, j, k) of the grid point (i, j, k) is:
d(i,j,k)＝d1+d2+…dN
it can be appreciated that in this embodiment, by summing the multiple transformation indexes, the formation transformation condition of each grid point can be accurately obtained, so that the scientificity of horizontal well deployment and the accuracy of horizontal well spatial analysis are improved.
Fig. 4 is a schematic flow chart of a shale gas horizontal well space analysis method according to another embodiment of the present application, as shown in fig. 4, the method further includes:
s301, acquiring a design track of drilling, and a related well-drilled track related to the drilling.
In the design of a well, analysis of the design trajectory of the well and the trajectory of the associated well is required to avoid the close spacing between the well and the well being drilled and to prevent the intersection of the well being drilled and the well being drilled.
In this embodiment, in order to analyze the distribution of the relevant well trajectory around the well, it is necessary to acquire the relevant well trajectory data and the design trajectory data at the well.
It should be noted that, the relevant well-drilled track data can be obtained after the fracturing of the horizontal well is completed, and the well-drilled design track data can be obtained through simulation according to the drawing of a designer and relevant software.
In addition, the determination of the drilled well related to the drilling is the same as the determination of the above embodiment, and this embodiment will not be repeated.
S302, discretizing the design track of the well drilling and the related track of the well drilling to obtain discrete data of the well drilling and discrete data of the well drilling.
In the present embodiment, in order to improve the efficiency of the horizontal well spatial analysis, the trajectory discretization may be performed on the trajectory that has been drilled at certain intervals, and the interval may be set to 1km, for example. The interval is determined by the interference prevention accuracy, which is not limited in the embodiment of the present application.
S303, calculating the shortest distance between the design track of the well and the related well-drilled track according to the well-drilled discrete data and the well-drilled discrete data.
The shortest distance between the design trajectory of the well and the associated drilled trajectory is calculated using a traversal algorithm from the drilled discrete data and the drilled discrete data.
For example, there may be a plurality of distance data between the trajectory of the well and any associated well trajectory, and the shortest distance between any associated well trajectory and the well may be used as the basis for the final data analysis.
S304, acquiring well leakage early warning information through shortest distance analysis.
In the embodiment of the application, when the shortest distance is smaller, the risk of well track crossing or lost circulation in the drilling is larger, because further adjustment of the design scheme of the drilling is needed, so that the shortest distance between the drilling track and the drilled track is kept within a reasonable safety range.
It should be noted that the above-mentioned safety range may be determined according to the situation of formation reconstruction around the well and the design scheme of the designer, which is not limited in this embodiment.
The embodiment of the application provides a shale gas horizontal well space analysis device which is used for executing the shale gas horizontal well space analysis method. Fig. 5 is a schematic structural diagram of a shale gas horizontal well space analysis device according to an embodiment of the present application, as shown in fig. 5, the device includes: a building unit 401, an extracting unit 402, and an analysis acquiring unit 403;
the establishing unit 401 is configured to establish a transformation index three-dimensional space model according to the given area and the depth range;
the extracting unit 402 is configured to extract, according to the transformation index three-dimensional space model, a transformation index slice along a planned deployment azimuth of the well;
the analysis and acquisition unit 403 is configured to analyze and acquire lost circulation early warning information according to the relevant spatial parameters of the drilled well and the transformation index slice.
Optionally, the establishing unit 401 is configured to establish a reconstruction index mathematical model with a spherical wave field according to the relevant spatial parameters that have been drilled.
Optionally, the analysis obtaining unit 403 is configured to calculate and obtain a modification index of each point in the formation by using the modification index mathematical model; calculating an edge layer arithmetic average slice and a root mean square slice of the transformation index slices according to the transformation index of each point in the stratum; and analyzing and acquiring well leakage early warning information according to the layer edge arithmetic average slice and the root mean square slice.
Alternatively, if there are multiple drilled formations that create a coincidence in the formation, the analysis acquisition unit 403 is configured to sum the multiple drilled formation indices and calculate an obtained formation index for each point in the formation.
Fig. 6 is a schematic structural diagram of a shale gas horizontal well space analysis device according to another embodiment of the present application, as shown in fig. 6, the device further includes: an acquisition unit 404, a discrete unit 405, and a calculation unit 406;
the acquiring unit 404 is configured to acquire a design trajectory of the in-drilling well and the related well-drilled trajectory related to the in-drilling well;
the discrete unit 405 is configured to perform discretization processing on the design trajectory of the in-drilling well and the related trajectory of the drilled well, so as to obtain the in-drilling discrete data and the drilled well discrete data;
the calculating unit 406 is configured to calculate, according to the on-drilling discrete data and the drilled discrete data, a shortest distance between the on-drilling design trajectory and the relevant drilled trajectory;
the analysis obtaining unit 403 is configured to obtain lost circulation early warning information through the shortest distance analysis.
Fig. 7 is a schematic structural diagram of a shale gas horizontal well space analysis device according to an embodiment of the present application, including: processor 610, storage medium 620, and bus 630, storage medium 620 storing machine-readable instructions executable by processor 610, processor 610 executing machine-readable instructions to perform steps of the above-described method embodiments when the electronic device is operating, processor 610 communicating with storage medium 620 over bus 630. The specific implementation manner and the technical effect are similar, and are not repeated here.
The present embodiments provide a storage medium having a computer program stored thereon, which when executed by a processor performs the above method.
In the several embodiments provided by the present invention, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in hardware plus software functional units.
The integrated units implemented in the form of software functional units described above may be stored in a computer readable storage medium. The software functional unit is stored in a storage medium, and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) or a processor (english: processor) to perform some of the steps of the methods according to the embodiments of the invention. And the aforementioned storage medium includes: u disk, mobile hard disk, read-Only Memory (ROM), random access Memory (Random Access Memory, RAM), magnetic disk or optical disk, etc.
The foregoing description is only of the preferred embodiments of the present application and is not intended to limit the same, but rather, various modifications and variations may be made by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims (10)
1. The shale gas horizontal well space analysis method is characterized by comprising the following steps of: establishing a transformation index three-dimensional space model according to the given area and the depth range;
extracting an reconstruction index slice along a planned deployment azimuth of drilling according to the reconstruction index three-dimensional space model;
analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice;
calculating the transformation index of each point in the stratum according to the transformation index mathematical model; let the half seam height of simulation seam in the transformation index mathematical model be H2, half seam length be L2, half seam width be W2, the coordinate of arbitrary point relative to perforation central point o be (L, W, H), then its transformation index d is:
wherein l<＝L 2 ，w<＝W 2 ，h<＝H 2 。
2. The method of claim 1, wherein before establishing the transformation index three-dimensional space model based on the given region and the depth range, further comprises: and establishing an index reconstruction mathematical model by using the spherical wave field according to the related well-drilled space parameters.
3. The method of claim 2, wherein analyzing the acquired lost circulation warning information based on the associated drilled spatial parameters and the retrofit index slices comprises: calculating and obtaining the transformation index of each point in the stratum by adopting the transformation index mathematical model;
calculating an edge layer arithmetic average slice and a root mean square slice of the transformation index slices according to the transformation index of each point in the stratum;
and analyzing and acquiring well leakage early warning information according to the layer edge arithmetic average slice and the root mean square slice.
4. The method of claim 3, wherein if there are a plurality of drilled formations that produce a coincidence in the formation, said calculating an improvement index for each point in the formation using the improvement index mathematical model comprises: and summing the plurality of the drilled modification indexes, and calculating the modification index of each point in the acquired stratum.
5. The method according to claim 1, wherein the method further comprises: acquiring the design track of the drilling well and the related well-drilled track related to the drilling well;
discretizing the design track of the on-drilling well and the related well-drilled track to obtain the discrete data of the on-drilling well and the discrete data of the well-drilled well;
calculating the shortest distance between the design track of the well and the related well-drilled track according to the well-drilled discrete data and the well-drilled discrete data;
and acquiring well leakage early warning information through the shortest distance analysis.
6. A shale gas horizontal well spatial analysis device, comprising: the device comprises a building unit, an extraction unit and an analysis acquisition unit;
the establishing unit is used for establishing a transformation index three-dimensional space model according to the given area and the depth range;
the extraction unit is used for extracting an index reconstruction slice along the planned deployment azimuth of the well drilling according to the index reconstruction three-dimensional space model;
the analysis acquisition unit is used for analyzing and acquiring lost circulation early warning information according to the related space parameters of the drilled well and the transformation index slice;
calculating the transformation index of each point in the stratum according to the transformation index mathematical model; let the half seam height of simulation seam in the transformation index mathematical model be H2, half seam length be L2, half seam width be W2, the coordinate of arbitrary point relative to perforation central point o be (L, W, H), then its transformation index d is:
wherein l<＝L 2 ，w<＝W 2 ，h<＝H 2 。
7. The apparatus according to claim 6, wherein the establishing unit is configured to establish the reconstruction index mathematical model as a spherical wave field based on the associated spatial parameters of the drilled well.
8. The apparatus of claim 6, wherein the apparatus further comprises: an acquisition unit, a discrete unit, and a calculation unit;
the acquisition unit is used for acquiring the design track of the drilling well and the related drilling well track related to the drilling well;
the discrete unit is used for discretizing the design track of the drilling well and the related track which is drilled well to obtain the discrete data of the drilling well and the discrete data of the drilled well;
the calculating unit is used for calculating the shortest distance between the design track of the well drilling and the related well drilling track according to the well drilling discrete data and the well drilling discrete data;
and the analysis acquisition unit is used for acquiring the lost circulation early warning information through the shortest distance analysis.
9. A shale gas horizontal well spatial analysis device, comprising: a processor, a storage medium, and a bus, the storage medium storing machine-readable instructions executable by the processor, the processor and the storage medium in communication over the bus when the apparatus is running, the processor executing the machine-readable instructions to perform the steps of the method of any one of claims 1-5.
10. A computer-readable storage medium, characterized in that the storage medium has stored thereon a computer program which, when executed by a processor, performs the steps of the method according to any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011114911.9A CN112233235B (en) | 2020-10-16 | 2020-10-16 | Shale gas horizontal well space analysis method, device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011114911.9A CN112233235B (en) | 2020-10-16 | 2020-10-16 | Shale gas horizontal well space analysis method, device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112233235A CN112233235A (en) | 2021-01-15 |
| CN112233235B true CN112233235B (en) | 2023-04-28 |
Family
ID=74118892
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011114911.9A Active CN112233235B (en) | 2020-10-16 | 2020-10-16 | Shale gas horizontal well space analysis method, device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112233235B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103699751A (en) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | Sand body reservoir architecture modeling method and system based on space vectors |
| CN106569272A (en) * | 2016-11-11 | 2017-04-19 | 东北石油大学 | Earthquake attribute fusion method based on data property space ascending dimension |
| CN106869790A (en) * | 2017-02-24 | 2017-06-20 | 中石化重庆涪陵页岩气勘探开发有限公司 | A kind of quick fine geology guidance method of shale gas horizontal well |
| CN106988739A (en) * | 2017-05-19 | 2017-07-28 | 中国石油集团川庆钻探工程有限公司 | Shale reservoir fracturing fracture identification and interpretation evaluation method |
-
2020
- 2020-10-16 CN CN202011114911.9A patent/CN112233235B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103699751A (en) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | Sand body reservoir architecture modeling method and system based on space vectors |
| CN106569272A (en) * | 2016-11-11 | 2017-04-19 | 东北石油大学 | Earthquake attribute fusion method based on data property space ascending dimension |
| CN106869790A (en) * | 2017-02-24 | 2017-06-20 | 中石化重庆涪陵页岩气勘探开发有限公司 | A kind of quick fine geology guidance method of shale gas horizontal well |
| CN106988739A (en) * | 2017-05-19 | 2017-07-28 | 中国石油集团川庆钻探工程有限公司 | Shale reservoir fracturing fracture identification and interpretation evaluation method |
Non-Patent Citations (5)
| Title |
|---|
| Slice-level diffusion encoding for motion and distortion correction;Jana Hutter等;《Medical Image Analysis》;20181231;第1-37页 * |
| 四川盆地涪陵焦石坝地区五峰―龙马溪组低序级断层识别技术及应用效果;武加鹤等;《石油实验地质》;20180128(第01期);第55-61页 * |
| 基于压裂监测的致密储层甜点识别;李新发等;《断块油气田》;20200123(第05期);第63-67页 * |
| 昭通页岩气示范区复杂地质条件下的地质导向技术;罗鑫等;《钻采工艺》;20180525(第03期);第6-7、40-43页 * |
| 涪陵页岩气田井地联合微地震监测气藏实例及认识;刘尧文等;《天然气工业》;20161231(第10期);第62-68页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112233235A (en) | 2021-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ren et al. | Modeling and simulation of complex fracture network propagation with SRV fracturing in unconventional shale reservoirs | |
| US9836561B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
| US9284829B2 (en) | Simulating subterranean fracture propagation | |
| US8437962B2 (en) | Generating probabilistic information on subterranean fractures | |
| CN113821953B (en) | Temporary plugging repeated fracturing design method based on dynamic ground stress and residual oil characteristics | |
| CA2945467C (en) | Enhancing reservoir characterization using real-time srv and fracture evolution parameters | |
| EP2904530B1 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| CN107220493A (en) | Shale gas horizontal well network fracture modeling method based on micro-seismic event | |
| CN104992468A (en) | Three-dimensional geological modeling method for fracture-cave carbonate reservoir | |
| EP2831804B1 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
| Yu et al. | Embedded discrete fracture model assisted study of gas transport mechanisms and drainage area for fractured shale gas reservoirs | |
| Rogers et al. | Characterising the in situ fragmentation of a fractured rock mass using a discrete fracture network approach | |
| CN112233235B (en) | Shale gas horizontal well space analysis method, device and storage medium | |
| CN114429054A (en) | Horizontal well staged fracturing simulation method and device, electronic equipment and medium | |
| RU2441145C1 (en) | Method of developing oil deposit with several oil accumulations located one above another | |
| Zeng et al. | Pore pressure disturbance induced by multistage hydraulic fracturing in shale gas: modelling and field application | |
| US12056422B1 (en) | System for designing perforating gun capable of reducing energy consumption | |
| Sepehrnoori et al. | EDFM for field-scale reservoir simulation with complex corner-point grids | |
| CN117910190A (en) | Well position assessment method, device, storage medium and equipment for broken solution oil and gas reservoir | |
| Zeng et al. | Research Article Pore Pressure Disturbance Induced by Multistage Hydraulic Fracturing in Shale Gas: Modelling and Field Application | |
| Wardhana et al. | Optimizing Field Development Strategy of Carbonate Gas Reservoir; An Integrated Subsurface Study of Alur Siwah Field | |
| AU2012396845B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| CN118153354A (en) | Method for optimizing interval between closely-cut volume fracturing clusters of land shale horizontal well | |
| CN118446127A (en) | Unconventional oil gas hydraulic fracturing and drainage integrated simulation method and device and electronic equipment | |
| Firsin et al. | Selection of the Optimal Technology of Re-Fracturing on Horizontal Wells with Multi-Stage Fracturing at the Field M |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20220601 Address after: 100020 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Applicant after: CHINA PETROLEUM & CHEMICAL Corp. Applicant after: SINOPEC CHONGQING FULING SHALE GAS EXPLORATION AND DEVELOPMENT Co.,Ltd. Address before: 408000 No. 6 Hefeng Avenue, Fuling District, Chongqing Applicant before: SINOPEC CHONGQING FULING SHALE GAS EXPLORATION AND DEVELOPMENT Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |