CN118036232A - Method and device for performing variogram analysis based on horizontal well data - Google Patents
Method and device for performing variogram analysis based on horizontal well data Download PDFInfo
- Publication number
- CN118036232A CN118036232A CN202211408698.1A CN202211408698A CN118036232A CN 118036232 A CN118036232 A CN 118036232A CN 202211408698 A CN202211408698 A CN 202211408698A CN 118036232 A CN118036232 A CN 118036232A
- Authority
- CN
- China
- Prior art keywords
- data
- horizontal well
- variogram
- thinning
- reservoir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004458 analytical method Methods 0.000 title claims abstract description 28
- 230000008859 change Effects 0.000 claims abstract description 49
- 230000000704 physical effect Effects 0.000 claims abstract description 47
- 230000006870 function Effects 0.000 claims description 42
- 238000004590 computer program Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 6
- 239000000243 solution Substances 0.000 description 5
- 230000002776 aggregation Effects 0.000 description 4
- 238000004220 aggregation Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 235000021184 main course Nutrition 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- -1 carbonate hydrocarbon Chemical class 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Landscapes
- Complex Calculations (AREA)
Abstract
The invention discloses a method and a device for performing variogram analysis based on horizontal well data. The method comprises the following steps: subdividing the reservoir along the direction from the top interface to the bottom interface based on the longitudinal change of the reservoir physical properties, and dividing the reservoir into a plurality of small layers so that the reservoir physical property change among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers; carrying out standardized processing on the horizontal well data of the line segments corresponding to each small layer; performing repeated thinning operation on the horizontal well data according to the direction perpendicular to the line segment; and carrying out variogram solving according to the data picked up by the thinning operation, and calculating a variogram corresponding to the small layer according to the result of the multiple variogram solving. The method can remove the influence of the physical property change of the longitudinal multi-period reservoir on the solution of the integral variation function, effectively remove the influence of the horizontal well data on the analysis variation function, and accurately calculate the variation function representing the development rule of the reservoir.
Description
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to a method and a device for performing variation function analysis based on horizontal well data.
Background
The horizontal well is an advanced technology process applied to oil extraction, and is a well with a certain length drilled along the horizontal direction, and the well bore of the horizontal well passes through an oil layer, so that the contact surface area of the well and the oil layer is increased, effective linear displacement is formed near the well bore of the horizontal well, the displacement effect is improved, and the horizontal well has the characteristics of high yield, low oil extraction cost and the like, and therefore, the international large carbonate hydrocarbon reservoir mainly selects a horizontal well water injection development mode at present. Unlike a vertical well, the horizontal well data has the characteristic of high aggregation, can effectively improve the inter-well prediction precision, and effectively represents the regularity of the geologic body. Because of the specificity of horizontal well data, sampling points are often concentrated in a dominant reservoir, and the situation of a real reservoir cannot be reflected if the sampling points participate in modeling, and the application of the horizontal well data in oil reservoir description and high-precision geological modeling is always a leading-edge subject in the world. The difficulty is that in the process of carrying out high-precision geologic modeling based on horizontal well data, the asymmetry of the medium around a well shaft caused by the horizontal well data and the characteristic that a large amount of data are gathered in a specific direction lead to that geologic statistics data often do not accord with geologic knowledge, and the error of reservoir variation function analysis is caused.
Disclosure of Invention
In view of the foregoing, the present invention has been made to provide a method and apparatus for performing a variogram analysis based on horizontal well data that overcomes or at least partially solves the foregoing problems.
In a first aspect, an embodiment of the present invention provides a method for performing a variogram analysis based on horizontal well data, including:
Further subdividing the reservoir along the direction from the top interface to the bottom interface, and dividing the reservoir into a plurality of small layers so that the physical property change of the reservoir among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
for each small layer, carrying out standardization processing on the horizontal well data of the line segment corresponding to the small layer;
performing multiple thinning operations on the horizontal well data on the plane of each small layer according to the direction perpendicular to the line segment;
And respectively carrying out variogram solving according to the data picked up by the multiple thinning operation, and calculating to obtain a variogram corresponding to the small layer according to the result of the multiple variogram solving.
In one embodiment, according to each small layer obtained by division, dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layer comprises:
and inserting each small layer trend surface according to the distribution trend of the layering depth of each small layer along the top interface to the bottom interface of the reservoir, and dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layer by utilizing each small layer trend surface.
In one embodiment, the normalizing the horizontal well data of the horizontal well original physical property curve segment corresponding to the small layer includes:
Establishing a reservoir physical field by utilizing interpolation of a vertical well reservoir physical curve, and generating a virtual physical curve along a horizontal well track grid;
fitting the virtual physical property curve with the original physical property curve line segment;
And adjusting the original physical property curve line segments by using the fitted parameters to obtain the standardized horizontal well data.
In one embodiment, performing a plurality of thinning operations on the horizontal well data according to a direction perpendicular to the line segment includes:
By arranging a plurality of parallel equidistant data thinning lines, the parallel equidistant data thinning lines are positioned on the plane of the small layer and are perpendicular to the line segments;
taking the width of the grid unit in the small layer as a moving step length, and gradually moving the data thinning lines; acquiring intersection points of the data thinning lines and the horizontal well data as data picked up by each thinning operation; until all grids in the small layer have been subjected to a thinning operation.
In one embodiment, the number of the data thinning lines is determined by the pitch of the data thinning lines and the width of the grid of the small layer, so as to satisfy the range that the plurality of parallel equidistant data thinning lines cover the grid of the small layer.
In one embodiment, the number of the thinning operations is obtained by the following formula:
A＝H/h；
a is the number of the thinning operation; h is the interval of the thinning lines; h is the mesh width.
In one embodiment, the method includes respectively performing solution of a variogram according to data picked up by multiple thinning operations, and calculating to obtain a variogram corresponding to the small layer according to a result of the solution of the variogram, where the method includes:
According to the data picked up each time, calculating a corresponding main change path and a corresponding secondary change path, determining a main change path direction, and determining a direction perpendicular to the main change path direction as a secondary change path direction;
Respectively solving a corresponding main variable path arithmetic average value and a corresponding secondary variable path arithmetic average value of the main variable path and the secondary variable path obtained by multiple times of calculation;
And characterizing a variation function corresponding to the small layer by using the main variation arithmetic average value, the secondary variation arithmetic average value, the main variation direction and the secondary variation direction.
In one embodiment, calculating the corresponding main course, sub-course, and main course direction, determining the main course direction, and determining a direction perpendicular to the main course direction as the sub-course direction includes:
determining the direction of the horizontal well as a main change path direction, and taking the direction perpendicular to the main change path direction as a secondary change path direction;
Counting the picked data according to the main change path direction and the secondary change path direction respectively, generating a point pair number-distance variation function histogram corresponding to the main change path and a point pair number-distance variation function histogram corresponding to the secondary change path, and drawing a half variance curve representing the logarithm of each point respectively;
Taking distance data corresponding to a point with half variance of 1 in a point-distance variation function histogram corresponding to a main variation as the main variation;
And taking distance data corresponding to a point with half variance of 1 in the point-to-number-distance variation function histogram corresponding to the secondary variation as the secondary variation.
In a second aspect, an embodiment of the present invention provides an apparatus for performing a variogram analysis based on horizontal well data, including:
the division module is used for further subdividing the reservoir along the direction from the top interface to the bottom interface and dividing the reservoir into a plurality of small layers so that the reservoir physical property change among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
the standardized processing module is used for carrying out standardized processing on the horizontal well data of the line segment corresponding to each small layer;
the thinning module is used for carrying out repeated thinning operation on the horizontal well data on the plane of each small layer according to the direction perpendicular to the line segment;
And the variogram analysis module is used for respectively carrying out variogram solving according to the data picked up by the multiple thinning operation and calculating to obtain the variogram corresponding to the small layer according to the result of the multiple variogram solving.
In a third aspect, an embodiment of the present invention provides a server, including: the system comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the method for carrying out the analysis of the variogram based on the horizontal well data when executing the program.
In a fourth aspect, an embodiment of the present invention provides a computer readable storage medium, where the computer readable storage medium stores a computer program, where the computer program when executed by a processor implements a method for performing a variogram analysis based on horizontal well data as described above.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
According to the method and the device for analyzing the variation function based on the horizontal well data, variation of physical properties of the reservoir in the longitudinal direction is reduced through the reservoir subdivision layers, the horizontal well data are subjected to standardized processing, the horizontal segment data of each small layer are regularly thinned and picked up, and the variation function of the physical properties of the reservoir is calculated by solving the variation function solving result after the data are picked up for each time; the method removes the influence of the physical property change of the longitudinal multi-period reservoir on the integral variation function, avoids the clustering effect caused by the high aggregation of the horizontal segment data, maximally uses the horizontal segment data, effectively characterizes the law of the geologic body, and reduces the uncertainty for the interwell simulation of the geologic model.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention. In the drawings:
FIG. 1 is a flow chart of a method for performing a variogram analysis based on horizontal well data in an embodiment of the present invention;
FIG. 2 is a schematic view of various small-scale trend surfaces for constructing a reservoir according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an acoustic wave field built using vertical well data interpolation in an embodiment of the present invention;
FIG. 4 is a schematic diagram of a fitting relationship between a horizontal well original property curve and field generated virtual property curve data in an embodiment of the present invention;
FIGS. 5-7 are schematic diagrams illustrating a data thinning operation using a data thinning line according to an embodiment of the present invention;
FIG. 8 is a graph illustrating a variation function histogram and a half variance curve in a thinning operation according to an embodiment of the present invention;
Fig. 9 is a block diagram of an apparatus for performing a variogram analysis based on horizontal well data according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
The inventors of the present invention have found that the variational function is an important concept characterizing the spatial correlation and the spatial variability of the variables. When the spatial distribution characteristics of the logging data are analyzed, parameters such as the variation range, the direction and the like of the variation function can provide key knowledge on aspects such as the direction of an object source, the sand spreading characteristics and the like. Because of the characteristics of the existing horizontal well data, when the analysis of the variation function is carried out by utilizing the horizontal well data, the clustering effect possibly caused by the high aggregation of the horizontal segment data can not objectively reflect the situation of a real reservoir.
In order to solve the problems in the prior art, an embodiment of the present invention provides a method for performing a variogram analysis based on horizontal well data, where an implementation process of the method may be shown in fig. 1, and the method includes the following steps:
S11, further subdividing the reservoir along the direction from the top interface to the bottom interface, and dividing the reservoir into a plurality of small layers so that the reservoir physical property change among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
s12, carrying out standardization processing on the horizontal well data of the line segments corresponding to each small layer;
s13, performing multiple thinning operations on the horizontal well data on the plane of each small layer according to the direction perpendicular to the line segment;
S14, respectively carrying out variational function solving according to the data picked up by the multiple thinning operation, and calculating to obtain a variational function corresponding to the small layer according to the result of the multiple variational function solving.
The steps S12-S14 are analysis processes of the variation function in each small layer, and the processing process is similar for each small layer, and the corresponding variation function needs to be obtained for each small layer.
Further, in the step S11, the reservoir is subdivided into small layers according to the characteristics of the deposition loops and the like by using the vertical well information, so as to obtain relatively homogeneous small layers, that is, ensure that the physical property change of the reservoir in each small layer is smaller as much as possible, so as to meet the preset conditions, for example, the parameter value of the physical property change of the reservoir in the small layer is in a preset smaller range, etc., which condition is not limited in the embodiment of the present invention; and then a small-layer trend surface is constructed according to the layering depth of the small layer and the spreading trend along the top-bottom interface of the reservoir, and the original physical property curve of the horizontal well is divided into at least one line segment corresponding to the small layer by utilizing the small-layer trend surface.
The original physical property curve of the horizontal well can be divided into a plurality of sections by each trend surface, the solution of the variation function is carried out in each section by the same method, and the three-dimensional body is converted into a plane body for processing, so that the longitudinal interference can be effectively removed.
Referring to fig. 2, there are a plurality of lines from the reservoir top interface to the reservoir bottom interface, which represent different small-layer trend surfaces, respectively, and the small-layer trend surfaces define each small layer, so that the original physical property curve of the horizontal well is divided into a plurality of line segments, and is distributed in 1-3 small layers.
Further, because the properties of the stratum on the upper side and the lower side of the horizontal well are different and are influenced by the stratum interface and the rock debris at the bottom of the well during drilling, the data of the original physical property curve are required to be subjected to standardized processing so as to prepare the data for the subsequent step, and the standardized processing ensures the accuracy of distinguishing good and bad reservoirs in the original physical property curve and also ensures that the well logging data of the horizontal well are influenced by the track as little as possible.
In the step S12, the standardization processing is performed on the horizontal well data of the horizontal well original physical property curve segment corresponding to the small layer, which may be specifically implemented by the following steps:
1) Establishing a reservoir physical field by utilizing interpolation of a vertical well reservoir physical curve, and generating a virtual physical curve along a horizontal well track grid;
2) Fitting the virtual physical property curve with the original physical property curve line segment;
3) And adjusting the original physical property curve line segments by using the fitted parameters to obtain the standardized horizontal well data.
Referring to FIG. 3, a schematic diagram of a reservoir property field-acoustic wave field is shown using vertical well reservoir property curve interpolation.
Fig. 4 is a schematic diagram showing a fitting relationship obtained by fitting a horizontal well original physical property curve and a virtual physical property curve generated by an acoustic wave field to each small layer data, wherein the abscissa dt_c in fig. 4 is an original acoustic wave curve value, and dtn_c is an acoustic wave curve value generated by a difference field. The lines in fig. 4 represent the fit obtained for each of the different microlayers.
An example of a fitting relationship (represented by a fitting formula) for one different microlayer is shown in table 1 below:
Table 1 different small-layer fitting formulas
| Small layer | Fitting relation |
| XC-1 | DTN_c＝+0.4382277*DT_c+41.48344 |
| XC-2 | DTN_c＝+0.247536*DT_c+55.98665 |
| XC-3 | DTN_c＝-0.2980585*DT_c+94.90803 |
| XC-4 | DTN_c＝+0.1014429*DT_c+67.94511 |
| XC-5 | DTN_c＝+0.0976647*DT_c+67.36352 |
| XC-6 | DTN_c＝+0.03686608*DT_c+76.00446 |
| XC-7 | DTN_c＝+0.1808889*DT_c+66.74947 |
| XC-8 | DTN_c＝+0.1523001*DT_c+67.94797 |
Parameters of the formulas in Table 1 above, such as +0.4382277 and 41.48344 for XC-1 small layers, were obtained by curve fitting.
Based on a fitting formula of the fitting relation of each small layer, substituting the data in the original physical property curve as the value of DT_c to obtain the value of DTN_c, namely the corresponding standardized data.
Further, in the step S13, a plurality of parallel equidistant data thinning lines are correspondingly arranged on the small layer, and the plurality of parallel equidistant data thinning lines are positioned on the plane of the small layer and are perpendicular to the line segments;
Taking the width of the grid unit in the small layer as a moving step length, and gradually moving the data thinning lines; each time a move is made, the following operations are performed: acquiring intersection points of the data thinning lines and the horizontal well data as data picked up by each thinning operation; until all grids in the small layer have been subjected to a thinning operation.
The number of the data thinning lines is determined by the distance between the data thinning lines and the width of the grids of the small layer so as to meet the requirement that the plurality of parallel equidistant data thinning lines cover the range of the grids of the small layer.
The thinning operation can ensure that all the pickup point data from the first time to the last time are not repeated;
in one embodiment, the number of the thinning operations is obtained by the following formula:
A＝H/h；
a is the number of the thinning operation; h is the interval of the thinning lines; h is the mesh width.
Referring to fig. 5-7, the thinning line and the horizontal well data are perpendicular to each other in the horizontal plane.
The numbers 1,2, 3, 4, 5 and 6 shown in fig. 7 represent the positions where the same thinning line moves 6 times, respectively, and fig. 4 shows 3 thinning lines, and each time the thinning lines move, the thinning lines are kept to move in parallel, the distance between the thinning lines is not changed, and the thinning lines move synchronously by the same step length (i.e. the width of one network element).
Referring to fig. 5, according to the direction of vertical horizontal well data, selecting a horizontal well pitch 300m as a thinning distance, setting parallel equidistant data thinning lines, performing horizontal well data thinning, wherein the intersection point of the data thinning lines and a horizontal well track is a thinned data pick-up point, so that the data of the pick-up points are uniformly distributed;
referring to fig. 6, the width of the grid cells of each small layer is 50m×50m, the thinning line interval is 300m, and the number of times of thinning should be calculated to be 6 times by the foregoing formula.
In one embodiment, in the step S14, the solution of the variogram is performed according to the data picked up by the multiple thinning operations, and the variogram corresponding to the small layer is calculated according to the result of the solution of the variogram, which is specifically implemented by the following manner:
According to the data picked up each time, calculating a corresponding main change path and a corresponding secondary change path, determining a main change path direction, and determining a direction perpendicular to the main change path direction as a secondary change path direction;
Respectively solving a corresponding main variable path arithmetic average value and a corresponding secondary variable path arithmetic average value of the main variable path and the secondary variable path obtained by multiple times of calculation;
and characterizing a variation function corresponding to the small layer by using the main variation arithmetic average value, the secondary variation arithmetic average value, the main variation direction and the secondary variation direction.
Taking the thinning operation of the previous example as 6 times as an example, calculating a main transformation path and a secondary transformation path according to data obtained by each operation, and carrying out arithmetic average on the main transformation paths obtained by the 6 times of calculation to obtain a main transformation path arithmetic average value; and carrying out arithmetic average on the secondary variation obtained by 6 times of calculation to obtain a secondary variation arithmetic average value.
Further, the step of calculating the corresponding main transformer path, secondary transformer path and main transformer path direction, determining the main transformer path direction and determining the direction perpendicular to the main transformer path direction as the secondary transformer path direction is realized by the following steps:
Determining the direction of the horizontal well as a main change path direction, and taking the direction vertical to the main change path direction as a secondary change path direction;
Counting the picked data according to the main change path direction and the secondary change path direction respectively, generating a point pair number-distance variation function histogram corresponding to the main change path and a point pair number-distance variation function histogram corresponding to the secondary change path, and drawing a half variance curve representing the logarithm of each point respectively;
Taking distance data corresponding to a point with half variance of 1 in a point-distance variation function histogram corresponding to a main variation as the main variation;
And taking distance data corresponding to a point with half variance of 1 in the point-to-number-distance variation function histogram corresponding to the secondary variation as the secondary variation.
Referring to fig. 8, the upper graph in fig. 8 is a histogram of the point-to-number-distance variation function corresponding to the main variation, the lower graph is a histogram of the point-to-number-distance variation function corresponding to the secondary variation, and half variance curves representing the pairs of points are respectively drawn, and for the upper graph, the distance corresponding to a point with a half variance of 1 on the half variance curve is the value of the main variation (for example, 2364 m); for the following graph, the distance corresponding to the point with half variance of 1 on the half variance curve is the value of the secondary range (for example 1896 m).
Taking the data shown in table 2 as an example, the data picked up for 6 times are calculated as the main variation and the sub-variation respectively, and the last column in table 2 is the arithmetic average, that is, the main variation calculated for 6 times is calculated and then arithmetically averaged to obtain the main variation of the small layer, and the sub-variation calculated for 6 times is calculated and then arithmetically averaged to obtain the sub-variation of the small layer.
TABLE 2
| Sequence number | Main transformer/m | Minor pass/m | Main transformer range direction/° |
| 1 | 2364 | 1896 | -65 |
| 2 | 2467 | 1946 | -65 |
| 3 | 2533 | 1460 | -65 |
| 4 | 2299 | 1739 | -65 |
| 5 | 2653 | 2018 | -65 |
| 6 | 2619 | 1920 | -65 |
| Average of | 2489 | 1830 | -65 |
The variational function can be characterized by the primary, secondary, primary and secondary directions of each minor layer.
Based on the same inventive concept, the embodiment of the invention also provides a device for performing variogram analysis based on the horizontal well data, and because the principles of the problems solved by the devices are similar to the method for performing variogram analysis based on the horizontal well data, the implementation of the device can be referred to the implementation of the method, and the repetition is omitted.
The embodiment of the invention provides a device for performing variogram analysis based on horizontal well data, which is shown by referring to fig. 9 and comprises:
The dividing module 91 is configured to further divide the reservoir into a plurality of small layers along a direction from the top interface to the bottom interface, so that a physical property change of the reservoir between the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
The normalization processing module 92 is configured to perform normalization processing on the horizontal well data of the line segment corresponding to each small layer;
The thinning module 93 is configured to perform a plurality of thinning operations on the horizontal well data on a plane of each small layer according to a direction perpendicular to the line segment;
The variogram analysis module 94 is configured to perform variogram solving according to the data picked up by the multiple thinning operations, and calculate a variogram corresponding to the small layer according to the result of the multiple variogram solving.
The embodiment of the invention provides a server, which comprises: the system comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize a method for performing variogram analysis based on horizontal well data.
Embodiments of the present invention provide a computer readable storage medium storing a computer program which, when executed by a processor, implements a method such as performing a variogram analysis based on horizontal well data. According to the method and the device for analyzing the variation function based on the horizontal well data, variation of physical properties of the reservoir in the longitudinal direction is reduced through the reservoir subdivision layers, the horizontal well data are subjected to standardized processing, the horizontal segment data of each small layer are regularly thinned and picked up, and the variation function of the physical properties of the reservoir is calculated by solving the variation function solving result after the data are picked up for each time; the method can remove the influence of the physical property change of the longitudinal multi-period reservoir on the whole variation function, avoid the clustering effect caused by the high aggregation of the horizontal segment data, furthest apply the horizontal segment data, effectively characterize the law of the geologic body and reduce the uncertainty for the interwell simulation of the geologic model.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, magnetic disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (11)
1. A method for performing a variogram analysis based on horizontal well data, comprising:
Further subdividing the reservoir along the direction from the top interface to the bottom interface, and dividing the reservoir into a plurality of small layers so that the physical property change of the reservoir among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
for each small layer, carrying out standardization processing on the horizontal well data of the line segment corresponding to the small layer;
performing multiple thinning operations on the horizontal well data on the plane of each small layer according to the direction perpendicular to the line segment;
And respectively carrying out variogram solving according to the data picked up by the multiple thinning operation, and calculating to obtain a variogram corresponding to the small layer according to the result of the multiple variogram solving.
2. The method of claim 1, wherein dividing the horizontal well original physical property curve into at least one line segment corresponding to the small layers according to each small layer obtained by division, comprises:
and inserting each small layer trend surface according to the distribution trend of the layering depth of each small layer along the top interface to the bottom interface of the reservoir, and dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layer by utilizing each small layer trend surface.
3. The method of claim 1, wherein the normalizing the horizontal well data of the horizontal well primitive property curve line segment corresponding to the small layer comprises:
Establishing a reservoir physical field by utilizing interpolation of a vertical well reservoir physical curve, and generating a virtual physical curve along a horizontal well track grid;
fitting the virtual physical property curve with the original physical property curve line segment;
And adjusting the original physical property curve line segments by using the fitted parameters to obtain the standardized horizontal well data.
4. The method of claim 1, wherein performing a plurality of thinning operations on the horizontal well data in a direction perpendicular to the line segment comprises:
By arranging a plurality of parallel equidistant data thinning lines, the parallel equidistant data thinning lines are positioned on the plane of the small layer and are perpendicular to the line segments;
taking the width of the grid unit in the small layer as a moving step length, and gradually moving the data thinning lines; acquiring intersection points of the data thinning lines and the horizontal well data as data picked up by each thinning operation; until all grids in the small layer have been subjected to a thinning operation.
5. The method of claim 4, wherein the number of data thinning-out lines is determined by a pitch of the data thinning-out lines and a width of the grid of the small layer to satisfy a range of the plurality of parallel equidistant data thinning-out lines covering the grid of the small layer.
6. The method of claim 4, wherein the number of the thinning operations is obtained by the following formula:
A＝H/h；
a is the number of the thinning operation; h is the interval of the thinning lines; h is the mesh width.
7. The method of claim 1, wherein the calculating the corresponding variogram of the small layer according to the result of the variogram solution performed on the data picked up by the multiple thinning operations, includes:
According to the data picked up each time, calculating a corresponding main change path and a corresponding secondary change path, determining a main change path direction, and determining a direction perpendicular to the main change path direction as a secondary change path direction;
Respectively solving a corresponding main variable path arithmetic average value and a corresponding secondary variable path arithmetic average value of the main variable path and the secondary variable path obtained by multiple times of calculation;
And characterizing a variation function corresponding to the small layer by using the main variation arithmetic average value, the secondary variation arithmetic average value, the main variation direction and the secondary variation direction.
8. The method of claim 7, wherein calculating the corresponding primary and secondary passes, determining the primary pass direction, and determining the direction perpendicular to the primary pass direction as the secondary pass direction, comprises:
determining the direction of the horizontal well as a main change path direction, and taking the direction perpendicular to the main change path direction as a secondary change path direction;
Counting the picked data according to the main change path direction and the secondary change path direction respectively, generating a point pair number-distance variation function histogram corresponding to the main change path and a point pair number-distance variation function histogram corresponding to the secondary change path, and drawing a half variance curve representing the logarithm of each point respectively;
Taking distance data corresponding to a point with half variance of 1 in a point-distance variation function histogram corresponding to a main variation as the main variation;
And taking distance data corresponding to a point with half variance of 1 in the point-to-number-distance variation function histogram corresponding to the secondary variation as the secondary variation.
9. An apparatus for performing a variogram analysis based on horizontal well data, comprising:
the division module is used for further subdividing the reservoir along the direction from the top interface to the bottom interface and dividing the reservoir into a plurality of small layers so that the reservoir physical property change among the small layers is smaller than a preset condition; dividing the original physical property curve of the horizontal well into at least one line segment corresponding to the small layers according to the divided small layers;
the standardized processing module is used for carrying out standardized processing on the horizontal well data of the line segment corresponding to each small layer;
the thinning module is used for carrying out repeated thinning operation on the horizontal well data on the plane of each small layer according to the direction perpendicular to the line segment;
And the variogram analysis module is used for respectively carrying out variogram solving according to the data picked up by the multiple thinning operation and calculating to obtain the variogram corresponding to the small layer according to the result of the multiple variogram solving.
10. A server, comprising: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method of performing a variogram analysis based on horizontal well data according to any one of claims 1 to 8 when the program is executed.
11. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program which, when executed by a processor, implements the method of performing a variogram analysis based on horizontal well data according to any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211408698.1A CN118036232A (en) | 2022-11-11 | 2022-11-11 | Method and device for performing variogram analysis based on horizontal well data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211408698.1A CN118036232A (en) | 2022-11-11 | 2022-11-11 | Method and device for performing variogram analysis based on horizontal well data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118036232A true CN118036232A (en) | 2024-05-14 |
Family
ID=90982715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211408698.1A Pending CN118036232A (en) | 2022-11-11 | 2022-11-11 | Method and device for performing variogram analysis based on horizontal well data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118036232A (en) |
-
2022
- 2022-11-11 CN CN202211408698.1A patent/CN118036232A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7873476B2 (en) | Well log correlation | |
| DE69725468T2 (en) | EXECUTION OF SCIENTIFIC EVALUATIONS WITH SIMULATED DATA | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| CN105184864A (en) | Site stratum three-dimensional geological structure model generation method for natural foundation replacement quantities calculation | |
| CN108561127A (en) | A kind of Formation pressure prediction method based on stochastic simulation | |
| CN110927793A (en) | Reservoir prediction method and system based on sequential random fuzzy simulation | |
| AU2021270865B2 (en) | Image-comparison based analysis of subsurface representations | |
| CN118036232A (en) | Method and device for performing variogram analysis based on horizontal well data | |
| CN111677486A (en) | Compact oil carbon dioxide huff and puff simulation method and device and storage medium | |
| CN104462792B (en) | Log data lithological stratum numerical value reduction method | |
| CN116187001A (en) | Automatic extraction method for effective leaching range of in-situ leaching uranium mining | |
| CN105989541B (en) | A kind of processing method of the well-log information based on big data | |
| CN112016956B (en) | Ore grade estimation method and device based on BP neural network | |
| CN114236624A (en) | Method and system for estimating fracturing modification space volume based on electromagnetic method | |
| GB2518172A (en) | Improvements in or relating to optimisation techniques | |
| CN105718705A (en) | Oil and gas field dynamic analysis method | |
| CN112859172A (en) | Longitudinal and transverse wave micro-logging data processing method and device | |
| CN112346117B (en) | Reservoir characteristic prediction method and device based on seismic attribute fusion | |
| CN115576011B (en) | Geological model-based phased reservoir prediction method, device and system | |
| CN113376697B (en) | Static correction method and device based on three-dimensional near-surface model | |
| CN108875163A (en) | A kind of method and system for evaluating three-dimensional cracking network connectivty | |
| CN112561111A (en) | Oilfield water flooding wave and coefficient prediction method and device | |
| CN112782757B (en) | Carbonate reservoir medium random modeling method | |
| CN112392469A (en) | Reservoir characteristic classification determination method and device | |
| CN117420599A (en) | River simulation and deep learning method for building river sand stratum grillwork |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination |