Horizontal well geosteering projection profile control method and device and storage medium

Abstract

The invention provides a horizontal well geosteering projection profile control method and device and a storage medium, and belongs to the technical field of horizontal well development. The method comprises the following steps: obtaining inflection point coordinates according to an obtained horizontal well trajectory curve, calculating an azimuth angle of a connecting line of two adjacent inflection points, extracting a plurality of single-direction sections from a geologic body model or an earthquake body according to the azimuth angle and the inflection point coordinates, and sequentially connecting and tiling the plurality of single-direction sections to form a multi-direction section. The method combines the geologic body model and the horizontal well trajectory curve, automatically calculates the inflection point of the fitting trajectory, uses the inflection point as the surface extraction basis, replaces the existing manual surface extraction positioning process, and uses the coordinates of the inflection point to quickly position the section position to complete the section extraction work along the trajectory. The working process is optimized, and the fitting property of the section and the track is improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of horizontal well development, and in particular relates to a method, a device and a storage medium for controlling a projection profile of a horizontal well geosteering. Background Art

[0002] Geosteering technology provides important support for fine trajectory control of horizontal well drilling and is one of the key technologies for the development of unconventional and complex oil and gas reservoirs. The geosteering process is generally divided into three steps, namely pre-drilling geological modeling, projection profile extraction, and trajectory and layer optimization adjustment. Among them, projection profile extraction realizes the transformation from a three-dimensional model to a two-dimensional plane. In the actual operation process, it can greatly improve the visibility of the formation information where the drill bit is located, and provide unified azimuth parameters for the calculation of values ​​such as formation apparent inclination and trajectory slope.

[0003] Currently, the projection profile methods provided by geological guidance software on the market are all single-azimuth profile extraction. This method is not applicable when the horizontal well azimuth turns greatly. Guidance personnel have to change the projection azimuth multiple times according to the actual drilling trajectory, which not only affects the continuity of guidance operations, but also causes display conflicts of interpretation results in different azimuths. Summary of the invention

[0004] The purpose of the present invention is to solve the above-mentioned problems existing in the prior art and to provide a method, device and storage medium for controlling the projection profile of a horizontal well geosteering, so as to realize the automatic generation of multi-directional profiles of the entire trajectory.

[0005] One of the purposes of the present invention is to provide a method for controlling a projection profile of a horizontal well geosteering.

[0006] A second object of the present invention is to provide a projection profile control device for horizontal well geosteering.

[0007] A third object of the present invention is to provide a computer storage medium.

[0008] The present invention is achieved through the following technical solutions:

[0009] The first aspect of the present invention provides a projection profile control method for horizontal well geo-steering, inflection point coordinates are obtained based on the acquired horizontal well trajectory curve, the azimuth of the line connecting two adjacent inflection points is calculated, and multiple single-azimuth profiles are extracted from the geological body model or seismic body based on the azimuth and inflection point coordinates, and the multiple single-azimuth profiles are connected in sequence and tiled to form a multi-azimuth profile.

[0010] A further improvement of the present invention is:

[0011] The method comprises the following steps:

[0012] Obtain horizontal well trajectory curves and geological body models and seismic body data of the operating block;

[0013] Calculate the coordinates of the inflection point of the horizontal well trajectory curve;

[0014] Calculate the azimuth of the line between two adjacent inflection points;

[0015] Extract multiple single-azimuth sections from geological volume models or seismic volume data;

[0016] Multiple single-azimuth sections are connected in sequence according to the horizontal well trajectory curve and tiled to form a multi-azimuth section.

[0017] A further improvement of the present invention is:

[0018] The calculating of the inflection point coordinates of the horizontal well trajectory curve specifically includes:

[0019] Step 21: Determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and set the landing point as the first inflection point;

[0020] Step 22: Smoothing the horizontal well trajectory curve to eliminate small-scale trajectory jitter;

[0021] Step 23: Taking the first inflection point as the starting point, thinning the smoothed horizontal well trajectory curve to obtain several inflection point coordinates.

[0022] A further improvement of the present invention is:

[0023] In step 23, the specific operation of thinning includes:

[0024] 1) Setting thresholds;

[0025] 2) Set the end point of the horizontal well trajectory curve as the new inflection point, and connect all current inflection points to form a thinning line;

[0026] 3) Connect two adjacent inflection points on the horizontal well trajectory curve to form a thinning line, determine the point on the horizontal well trajectory curve between the two adjacent inflection points that is farthest from the thinning line, and if the distance from the point to the thinning line is greater than a threshold, set the point as a new inflection point;

[0027] 4) Repeat 3) until the distance between the point between two adjacent inflection points on the horizontal well trajectory curve and the thinning line is no greater than the threshold, then end.

[0028] A further improvement of the present invention is:

[0029] 3) Connect two adjacent inflection points on the horizontal well trajectory curve to form a thinning line, determine the point on the horizontal well trajectory curve between the two adjacent inflection points that is farthest from the thinning line, and if the distance from the point to the thinning line is greater than a threshold, set the point as a new inflection point;

[0030] Specific operations include:

[0031] Determine the coordinates of two adjacent inflection points of the horizontal well trajectory curve;

[0032] The equation of the thinning line formed by the two inflection points is:

[0033] y=kx+b

[0034] Where x and y are x and y coordinates respectively, and k and b are constants;

[0035] Assume that the coordinates of any point between two adjacent inflection points on the horizontal well trajectory curve are (X t ,Y t ), then the distance from this point to the thinning line is:

[0036]

[0037] The maximum distance value obtained by the above calculation is found, and it is determined whether the maximum distance value is greater than a threshold value. If so, the point corresponding to the maximum distance value is set as a new inflection point.

[0038] The improvement of the present invention is:

[0039] The specific operation of calculating the azimuth of the line connecting two adjacent inflection points includes:

[0040] The azimuth angle α of the line connecting adjacent inflection points is calculated by the following formula:

[0041]

[0042] Among them, x i is the x coordinate of the first inflection point between the two connected inflection points; i+1 is the x coordinate of the latter of the two connected inflection points, y i is the y coordinate of the first inflection point between the two connected inflection points, y i+1 It is the y coordinate of the latter inflection point between two connected inflection points.

[0043] A further improvement of the present invention is:

[0044] The specific operations of extracting multiple single-azimuth profiles from the geological body model or seismic body data include:

[0045] Sampling is performed between two adjacent inflection points in the horizontal well trajectory curve to obtain multiple sampling points, which are marked as 1, 2, 3, ..., N in sequence;

[0046] According to the azimuth of the line between two adjacent inflection points, the coordinates of the sampling points between two adjacent inflection points are calculated in sequence:

[0047]

[0048] Where n is the coordinate point number, 1, 2, 3...N, x new_n ,y new_n are the x and y coordinates of the nth coordinate point between two adjacent inflection points, respectively. i ,y i are the x and y coordinates of the previous inflection point between two adjacent inflection points, relv is the sampling interval, and α is the azimuth;

[0049] According to the coordinates of the calculated sampling points, the attribute data of the corresponding coordinate points are extracted from the geological body model or seismic body data to form a profile.

[0050] A second aspect of the present invention provides a projection profile control device for horizontal well geosteering, comprising:

[0051] An acquisition unit, used to acquire the horizontal well trajectory curve and the geological body model and seismic body data of the operating block;

[0052] An inflection point coordinate determination unit, used to calculate the inflection point coordinates of the horizontal well trajectory curve;

[0053] An azimuth angle calculation unit, used to calculate the azimuth angle of a line between two adjacent inflection points;

[0054] A single-azimuth profile extraction unit is used to extract multiple single-azimuth profiles from a geological body model or seismic body data;

[0055] The multi-azimuth profile generation unit is used to sequentially connect multiple single-azimuth profiles according to the horizontal well trajectory curve and lay them out to form the final multi-azimuth profile.

[0056] A further improvement of the present invention is:

[0057] The inflection point coordinate determining unit comprises:

[0058] The first inflection point setting subunit is used to determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and the landing point is set as the first inflection point;

[0059] The smoothing subunit is used to smooth the horizontal well trajectory curve and eliminate small-scale trajectory jitter;

[0060] The thinning subunit is used to thin the smoothed horizontal well trajectory curve with the first inflection point as the starting point to obtain the coordinates of several inflection points.

[0061] The third aspect of the present invention provides a computer storage medium, which stores at least one computer-executable program, and when the at least one program is executed by the computer, the computer executes the steps in the above-mentioned horizontal well geo-steering projection profile control method.

[0062] Compared with the prior art, the present invention has the following beneficial effects:

[0063] The present invention combines the geological body model and the horizontal well trajectory curve, automatically calculates the inflection point of the fitting trajectory, uses the inflection point as the basis for surface extraction, replaces the existing manual positioning and surface extraction process, and simultaneously uses the inflection point coordinates to quickly locate the section position to complete the section extraction along the trajectory, thereby optimizing the workflow and improving the fit between the section and the trajectory.

[0064] The present invention performs automatic extraction and synthesis of multi-directional projection profiles according to the horizontal well trajectory curve, quickly fits the appropriate projection orientation of each well section by analyzing the trajectory trend and forms a synthetic profile, provides complete operation information for guidance personnel on a two-dimensional plane, and improves the efficiency of geological guidance interpretation and information integrity. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 A flow chart of a projection profile control method for geosteering of a horizontal well is provided for the present invention;

[0066] Figure 2 This is an inflection point effect diagram generated based on the horizontal well trajectory curve;

[0067] Figure 3 It is a multi-directional section horizontal plane tiling effect diagram (horizontal projection diagram);

[0068] Figure 4 This is a multi-directional cross-sectional projection rendering. DETAILED DESCRIPTION

[0069] The present invention is further described in detail below in conjunction with the accompanying drawings:

[0070] The present invention provides a method for controlling a projection profile of a horizontal well geo-steering. The method obtains the coordinates of the trajectory inflection point according to the acquired horizontal well trajectory curve, calculates the azimuth of the line connecting two adjacent inflection points, extracts a profile from a geological body model or a seismic body according to the azimuth angle and the inflection point coordinates, and connects and flattens all the profiles in sequence to form a final profile.

[0071] The embodiments of the method of the present invention are as follows:

[0072] [Example 1]

[0073] like Figure 1 As shown, the method comprises the following steps:

[0074] Step 1, obtaining the horizontal well trajectory curve and the geological body model and seismic body data of the operating block;

[0075] Step 2, calculating the coordinates of the inflection point of the horizontal well trajectory curve;

[0076] Preferably, step 2 comprises the following steps:

[0077] Step 21: Determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and set the landing point as the first inflection point.

[0078] Specifically, the position where the trajectory curve is inserted into the upper interface of the target formation is taken as the landing point (the upper interface of the target formation comes from formation modeling), and the landing point is the point where the trajectory curve enters the target formation.

[0079] Step 22: Smoothing the horizontal well trajectory curve to eliminate small-scale trajectory jitter;

[0080] The smoothing process may adopt any existing smoothing process method or algorithm, such as five-point smoothing, Gaussian smoothing, etc., which will not be described in detail here.

[0081] Step 23: Taking the first inflection point as the starting point, thin out the smoothed horizontal well trajectory curve to obtain several inflection point coordinates. The results are as follows: Figure 2 As shown, Figure 2 The midline is the trajectory, and the points on the line are the inflection points;

[0082] The thinning process is based on the Douglas-Peucker thinning theory, and the specific operations include:

[0083] 1) Set the threshold E. Generally, the threshold E is set within 10;

[0084] 2) Set the end point of the horizontal well trajectory curve as a new inflection point, and connect all current inflection points to form a thinning line (for example, there are only two inflection points, namely the starting point and the end point);

[0085] 3) Connect two adjacent inflection points on the horizontal well trajectory curve to form a thinning line, determine the point on the horizontal well trajectory curve between the two adjacent inflection points that is farthest from the thinning line, and if the distance from the point to the thinning line is greater than a threshold value E, set the point as a new inflection point;

[0086] Specific operations include:

[0087] Determine two adjacent inflection points of the horizontal well trajectory curve;

[0088] The equation of the thinning line formed by the two inflection points is:

[0089] y=kx+b

[0090] Where x and y are x and y coordinates respectively, and k and b are constants;

[0091] Assume that the coordinates of any point between two adjacent inflection points on the horizontal well trajectory curve are (X t ,Y t ), then the distance from this point to the thinning line is:

[0092]

[0093] Among them, k and b are two constants in the rarefaction line equation;

[0094] Find the maximum distance value obtained by the above calculation, and determine whether the maximum distance value is greater than a threshold value. If so, set the point corresponding to the maximum distance value as a new inflection point;

[0095] 4) Repeat 3) until the distance between the point between two adjacent inflection points on the horizontal well trajectory curve and the thinning line is no greater than the threshold E, then end, then all the inflection points obtained in the process are the inflection points obtained in this step.

[0096] Step 3, calculating the azimuth of the line between two adjacent inflection points;

[0097] Specific operations include:

[0098] The azimuth angle α of the line connecting adjacent inflection points is calculated by the following formula:

[0099]

[0100] Among them, x i is the x coordinate of the first inflection point between the two connected inflection points; i+1 is the x coordinate of the latter of the two connected inflection points, y i is the y coordinate of the first inflection point between the two connected inflection points, y i+1 It is the y coordinate of the latter inflection point between two connected inflection points.

[0101] Step 4, extracting multiple single-azimuth sections from the geological body model or seismic body data;

[0102] Among them, whether it is a geological body model or a seismic body, the process of extracting the surface is the same, specifically:

[0103] First, sampling is performed between two adjacent inflection points in the horizontal well trajectory curve to obtain multiple sampling points, which are marked as 1, 2, 3, ..., N in sequence;

[0104] According to the azimuth between two adjacent inflection points, the coordinates of the sampling points between two adjacent inflection points are calculated in sequence:

[0105]

[0106] Where n is the coordinate point number, 1, 2, 3...N, x new_n ,y new_n is the x and y coordinates of the nth coordinate point between two adjacent inflection points, x i ,y i are the x and y coordinates of the previous inflection point between two adjacent inflection points, relv is the sampling interval (usually the default is 20 meters per point), and α is the azimuth;

[0107] According to the coordinates of the calculated sampling points, the attribute data of the corresponding coordinate points are extracted from the geological body model or seismic body data to form a profile.

[0108] The section formed between two adjacent inflection points is a unidirectional section. According to the order of inflection points, the sections between two adjacent inflection points on the horizontal well trajectory curve are extracted in sequence to form multiple unidirectional sections.

[0109] like Figure 3 The dotted line between two adjacent inflection points is the section extracted from the geological body model, and the sections between the inflection points are connected in sequence to form a multi-directional section (horizontal projection diagram).

[0110] Step 5: Connect multiple single-direction sections in sequence according to the trajectory curve and lay them out to form the final section. Figure 4 The final result in the cross section.

[0111] The present invention combines the geological body model and the horizontal well trajectory curve, automatically calculates the inflection point of the fitting trajectory, uses the inflection point as the basis for surface extraction, replaces the existing manual positioning and surface extraction process, and simultaneously uses the inflection point coordinates to quickly locate the section position to complete the section extraction along the trajectory, thereby realizing the automatic generation of multi-directional sections of the entire trajectory, optimizing the workflow, and improving the fit between the section and the trajectory.

[0112] The present invention performs automatic extraction and synthesis of multi-directional projection profiles according to the horizontal well trajectory curve, quickly fits the appropriate projection orientation of each well section by analyzing the trajectory trend and forms a synthetic profile, provides complete operation information for guidance personnel on a two-dimensional plane, and improves the efficiency of geological guidance interpretation and information integrity.

[0113] The present invention also provides a projection profile control device for horizontal well geosteering, and the embodiments of the device are as follows:

[0114] [Example 2]

[0115] The device comprises:

[0116] An acquisition unit, used to acquire the horizontal well trajectory curve and the geological body model and seismic body data of the operating block;

[0117] An inflection point coordinate determination unit, used to calculate the inflection point coordinates of the horizontal well trajectory curve;

[0118] An azimuth angle calculation unit, used to calculate the azimuth angle of a line between two adjacent inflection points;

[0119] A single-azimuth profile extraction unit is used to extract multiple single-azimuth profiles from a geological body model or seismic body data;

[0120] The multi-azimuth profile generation unit is used to sequentially connect multiple single-azimuth profiles according to the horizontal well trajectory curve and lay them out to form the final multi-azimuth profile.

[0121] Specifically, the inflection point coordinate determining unit includes:

[0122] The first inflection point setting subunit is used to determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and the landing point is set as the first inflection point;

[0123] The smoothing subunit is used to smooth the horizontal well trajectory curve and eliminate small-scale trajectory jitter;

[0124] The thinning subunit is used to thin the smoothed horizontal well trajectory curve with the first inflection point as the starting point to obtain the coordinates of several inflection points.

[0125] The present invention also provides a computer storage medium, and the embodiments of the computer storage medium are as follows:

[0126] [Example 3]

[0127] The computer storage medium stores at least one computer-executable program, and when the at least one program is executed by the computer, the computer executes the steps in the above-mentioned horizontal well geosteering projection profile control method.

[0128] The above technical solution is only one implementation mode of the present invention. For those skilled in the art, it is easy to make various types of improvements or modifications based on the principles disclosed in the present invention, and it is not limited to the technical solution described in the above specific embodiments of the present invention. Therefore, the above description is only preferred and does not have a restrictive meaning.

Claims

1. A method for controlling a projection profile of a horizontal well geosteering, characterized in that: The inflection point coordinates are obtained according to the acquired horizontal well trajectory curve, the azimuth of the line connecting two adjacent inflection points is calculated, and multiple single-azimuth sections are extracted from the geological body model or seismic body according to the azimuth and inflection point coordinates. The multiple single-azimuth sections are connected in sequence and tiled to form a multi-azimuth section.

2. The method for controlling the projection profile of horizontal well geosteering according to claim 1, characterized in that: The method comprises the following steps: Obtain horizontal well trajectory curves and geological body models and seismic body data of the operating block; Calculate the coordinates of the inflection point of the horizontal well trajectory curve; Calculate the azimuth of the line between two adjacent inflection points; Extract multiple single-azimuth sections from geological volume models or seismic volume data; Multiple single-azimuth sections are connected in sequence according to the horizontal well trajectory curve and tiled to form a multi-azimuth section.

3. The method for controlling the projection profile of horizontal well geosteering according to claim 2, characterized in that: The calculating of the inflection point coordinates of the horizontal well trajectory curve specifically includes: Step 21: Determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and set the landing point as the first inflection point; Step 22: Smoothing the horizontal well trajectory curve; Step 23: Taking the first inflection point as the starting point, thinning the smoothed horizontal well trajectory curve to obtain several inflection point coordinates.

4. The method for controlling the projection profile of horizontal well geosteering according to claim 3, characterized in that: In step 23, the specific operation of thinning includes: 1) Setting thresholds; 2) Set the end point of the horizontal well trajectory curve as the new inflection point, and connect all current inflection points to form a thinning line; 3) Connect two adjacent inflection points on the horizontal well trajectory curve to form a thinning line, determine the point on the horizontal well trajectory curve between the two adjacent inflection points that is farthest from the thinning line, and if the distance from the point to the thinning line is greater than a threshold, set the point as a new inflection point; 4) Repeat 3) until the distance between the point between two adjacent inflection points on the horizontal well trajectory curve and the thinning line is no greater than the threshold, then end.

5. The method for controlling the projection profile of horizontal well geosteering according to claim 4, characterized in that: Connect two adjacent inflection points on the horizontal well trajectory curve to form a thinning line, determine the point on the horizontal well trajectory curve between the two adjacent inflection points that is farthest from the thinning line, and if the distance from the point to the thinning line is greater than a threshold, set the point as a new inflection point. The specific operations include: Determine the coordinates of two adjacent inflection points of the horizontal well trajectory curve; The equation of the thinning line formed by the two inflection points is: y=kx+b Where x and y are x and y coordinates respectively, and k and b are constants; Assume that the coordinates of any point between two adjacent inflection points on the horizontal well trajectory curve are (X t ,Y t ), then the distance from this point to the thinning line is: The calculated maximum distance value is determined, and it is judged whether the maximum distance value is greater than a threshold value. If so, the point corresponding to the maximum distance value is set as a new inflection point.

6. The method for controlling the projection profile of horizontal well geosteering according to claim 2, characterized in that: The specific operation of calculating the azimuth of the line connecting two adjacent inflection points includes: The azimuth angle α of the line connecting adjacent inflection points is calculated by the following formula: Among them, x i is the x coordinate of the first inflection point between the two connected inflection points; i+1 is the x coordinate of the latter of the two connected inflection points, y i is the y coordinate of the first inflection point between the two connected inflection points, y i+1 It is the y coordinate of the latter inflection point between two connected inflection points.

7. The method for controlling the projection profile of horizontal well geosteering according to claim 2, characterized in that: The specific operations of extracting multiple single-azimuth profiles from the geological body model or seismic body data include: Sampling is performed between two adjacent inflection points in the horizontal well trajectory curve to obtain multiple sampling points, which are marked as 1, 2, 3, ..., N in sequence; According to the azimuth of the line between two adjacent inflection points, the coordinates of the sampling points between two adjacent inflection points are calculated in sequence: Where n is the coordinate point number, 1, 2, 3...N, x new_n ,y new_n are the x and y coordinates of the nth coordinate point between two adjacent inflection points, respectively. i ,y i are the x and y coordinates of the previous inflection point between two adjacent inflection points, relv is the sampling interval, and α is the azimuth; According to the coordinates of the calculated sampling points, the attribute data of the corresponding coordinate points are extracted from the geological body model or seismic body data to form a profile.

8. A projection profile control device for horizontal well geosteering, characterized in that: include: An acquisition unit, used to acquire the horizontal well trajectory curve and the geological body model and seismic body data of the operating block; An inflection point coordinate determination unit, used to calculate the inflection point coordinates of the horizontal well trajectory curve; An azimuth angle calculation unit, used to calculate the azimuth angle of a line between two adjacent inflection points; A single-azimuth profile extraction unit is used to extract multiple single-azimuth profiles from a geological body model or seismic body data; The multi-azimuth profile generation unit is used to connect multiple single-azimuth profiles in sequence according to the horizontal well trajectory curve, and lay them out to form the final multi-azimuth profile.

9. The projection profile control device for horizontal well geosteering according to claim 8, characterized in that: The inflection point coordinate determining unit comprises: The first inflection point setting subunit is used to determine the landing point of the horizontal well according to the horizontal well trajectory curve and the formation modeling in the geological body model, and the landing point is set as the first inflection point; The smoothing subunit is used to smooth the horizontal well trajectory curve and eliminate small-scale trajectory jitter; The thinning subunit is used to thin the smoothed horizontal well trajectory curve with the first inflection point as the starting point to obtain the coordinates of several inflection points.

10. A computer storage medium, characterized in that: The computer storage medium stores at least one computer-executable program, and when the at least one program is executed by the computer, the computer executes the steps in the horizontal well geosteering projection profile control method as described in any one of claims 1-7.

Images