Interwell profile interwell data interpolation and interwell profile map generation method and device

Abstract

The application discloses a method and device for well profile interwell data interpolation and well profile drawing. The method comprises the following steps: assigning a probability density function to attribute data to be interpolated in each well section of a well profile, and determining a probability distribution function corresponding to the probability density function of each well section; connecting points with equal probability distribution function values on adjacent well sections to obtain data pair connecting lines of the adjacent well sections; and determining attribute values to be interpolated of points to be interpolated on the data pair connecting lines according to attribute values to be interpolated on the adjacent two well sections corresponding to the data pair connecting lines. The probability distribution function is used to realize interwell data interpolation, and the method does not depend on the spatial distribution of known data, thereby avoiding the adverse effects of uneven and unfixed known data on the interpolation results and realizing accurate control of the interpolation results.

Description

Technical Field

[0001] This invention relates to the field of geological engineering technology, and in particular to a method and apparatus for interpolating well profile data and generating well profile maps. Background Technology

[0002] In geological research, it is common practice to interpolate unknown geological attribute data between wells using known geological attribute data from the well site on a planar surface, thereby obtaining the planar distribution trend of relevant geological attribute data in the well area. Commonly used planar interpolation methods include the minimum curvature algorithm, the Kriging algorithm, the moving average algorithm, and the inverse distance multiplication method. These interpolation methods generally use known data distributed along the well to interpolate unknown data distributed along the strata. Summary of the Invention

[0003] The inventors of this application have discovered that, for inter-well data interpolation in well-connected profiles, due to the stratified characteristics of formation sedimentation, the distribution of various properties (porosity, permeability, etc.) along the layers between wells typically exhibits a much greater lateral continuity than vertical continuity. However, the spatial location of known data shows the opposite trend: the density of vertically distributed known data along the well is much greater than the density of laterally distributed known data along the layers. Furthermore, when the well-connected profile is flattened or its longitudinal and transverse scales are changed, the spatial location of the known data will change drastically. Therefore, existing planar interpolation methods, whose calculation results heavily depend on the spatial distribution of known data, are not suitable for interpolating inter-well data in well-connected profiles.

[0004] In view of the above problems, the present invention is proposed to provide a method and apparatus for well-to-well profile data interpolation and well-to-well profile generation that overcomes or at least partially solves the above problems.

[0005] This invention provides a method for inter-well data interpolation in well profiles, comprising:

[0006] Assign a probability density function to the attribute data to be interpolated for each well segment in the well profile, and determine the probability distribution function corresponding to the probability density function of each well segment.

[0007] Connect the points with equal probability distribution function values ​​on adjacent well sections to obtain the data pair connection lines for adjacent well sections;

[0008] Based on the interpolation attribute values ​​corresponding to the data pairs on two adjacent well sections, determine the interpolation attribute values ​​of the points to be interpolated on the data pairs.

[0009] In some optional embodiments, assigning a probability density function to the interpolated attribute data of each well segment in the well profile, and determining the probability distribution function corresponding to the probability density function of each well segment, includes:

[0010] For the to-be-interpolated attribute of the well section k in the well tie profile, a probability density function f k (x k ) is given, where the value range of x k is the length of the well section;

[0011] The probability density function f k (x k ) of the well section k is integrated to obtain the probability distribution function F k (x k ) of the well section k.

[0012] In some optional embodiments, the connecting of the points with equal probability distribution function values on adjacent well sections to obtain the data pair connecting lines of the adjacent well sections comprises:

[0013] For the well section k and the well section s in the well tie profile, the points satisfying F k (x k ) = F s (x s ) are connected to obtain a plurality of data pair connecting lines between the well section k and the well section s.

[0014] In some optional embodiments, according to the to-be-interpolated attribute values of the data pair connecting lines on the corresponding to-be-interpolated attribute values of the adjacent two well sections, the to-be-interpolated attribute value of the to-be-interpolated point on the data pair connecting line is determined, comprising:

[0015] According to the distance a of the to-be-interpolated point to the well section k, the distance b of the to-be-interpolated point to the well section s, and the probability distribution value m of the to-be-interpolated point on the well section k and the probability distribution value n of the to-be-interpolated point on the well section s, the probability distribution value of the end point on the data pair connecting line is determined.

[0016] According to the probability distribution value of the end point, the position of the end point is determined, and according to the position of the end point, the to-be-interpolated attribute value, and the position of the to-be-interpolated point, the to-be-interpolated attribute value of the to-be-interpolated point is determined.

[0017] In some optional embodiments, the above method further comprises:

[0018] The well tie profile to be interpolated, the horizon, and the to-be-interpolated attribute are determined, the to-be-interpolated attribute data of the well section included in the well tie profile to be interpolated in the horizon range is obtained, and the to-be-interpolated attribute data of each well section in the well tie profile is obtained; the to-be-interpolated attribute data comprises the position coordinates of each position on the well section and the corresponding to-be-interpolated attribute value.

[0019] In some optional embodiments, after obtaining the to-be-interpolated attribute data of each well section in the well tie profile, the method further comprises:

[0020] The pre-processing of the to-be-interpolated attribute data of each well section in the well-to-well profile is performed: for unknown positions on the well section that need to be supplemented with to-be-interpolated attribute values, the to-be-interpolated attribute values of the unknown positions are determined and supplemented based on the to-be-interpolated attribute values of the adjacent positions on both sides.

[0021] The embodiment of the present application also provides a well-to-well profile generation method, which comprises: determining the to-be-interpolated attribute values of each to-be-interpolated point between adjacent well sections in the well-to-well profile by using the above method, and generating a well-to-well profile based on the to-be-interpolated attribute values of each to-be-interpolated point.

[0022] The embodiment of the present application also provides a well-to-well profile interwell data interpolation device, which comprises:

[0023] A pre-processing module is configured to assign a probability density function to the to-be-interpolated attribute data of each well section in the well-to-well profile, and determine a probability distribution function corresponding to the probability density function of each well section, respectively.

[0024] A connecting line module is configured to connect the points with equal probability distribution function values on adjacent well sections to obtain a data pair connecting line of the adjacent well sections.

[0025] An interpolation module is configured to determine the to-be-interpolated attribute values of the to-be-interpolated points on the data pair connecting line according to the to-be-interpolated attribute values of the data pair connecting line on adjacent two well sections.

[0026] The embodiment of the present application also provides a well-to-well profile generation device, which comprises:

[0027] The above well-to-well profile interwell data interpolation device is configured to determine the to-be-interpolated attribute values of each to-be-interpolated point between adjacent well sections in the well-to-well profile.

[0028] A generation module is configured to generate a well-to-well profile based on the to-be-interpolated attribute values of each to-be-interpolated point.

[0029] The embodiment of the present application also provides a computer storage medium, which stores computer executable instructions, and the computer executable instructions are executed by a processor to implement the above well-to-well profile interwell data interpolation method and / or the above well-to-well profile generation method.

[0030] The embodiment of the present application also provides a device, which comprises a memory, a processor and a computer program stored on the memory and executable on the processor, and the processor implements the above well-to-well profile interwell data interpolation method and / or the above well-to-well profile generation method when executing the program.

[0031] The above technical solution provided by the embodiment of the present application has at least the following beneficial effects:

[0032] The well-to-well profile interwell data interpolation method provided by the embodiment of the present application is based on the probability density function and the corresponding probability distribution function of the attribute data to be interpolated of each well section in the well-to-well profile to perform interpolation calculation, and the probability distribution function is used to realize one-to-one correspondence of the interwell data, which essentially overcomes the adverse effects of the uneven distribution of known data and the non-fixed relative position of the well-to-well profile on the interpolation results, and the interpolation results of the well-to-well profile will no longer change due to the flattening of each layer profile and the adjustment of horizontal and vertical scale, so that the interpolation algorithm has the characteristics of stability; and the adjustment of the probability density function greatly realizes the accurate control of the interpolation results and the convenient display of geological understanding, so that the interpolation algorithm has the characteristics of flexibility.

[0033] Other features and advantages of the present application will be set forth in the following description, and in part will become apparent to those skilled in the art from the description, or can be learned by practice of the present application. The objects and other advantages of the present application can be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.

[0034] The technical solutions of the present application will be further described in detail below with the help of the accompanying drawings and embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0035] The accompanying drawings are included to provide a further understanding of the present application, and constitute a part of the specification, illustrate the present application, and are used to explain the present application together with the embodiments of the present application, and do not constitute a limitation on the present application. In the drawings:

[0036] Figure 1 The flow chart of the well-to-well profile interwell data interpolation method in the embodiment one of the present application is shown in the figure;

[0037] Figure 2 The flow chart of the well-to-well profile interwell data interpolation method in the embodiment two of the present application is shown in the figure;

[0038] Figure 3 The principle diagram of the well-to-well profile interwell data interpolation method in the embodiment two of the present application is shown in the figure;

[0039] Figure 4 The flow chart of the well-to-well profile generation method in the embodiment three of the present application is shown in the figure;

[0040] Figure 5 The structure diagram of the well-to-well profile interwell data interpolation device in the embodiment of the present application is shown in the figure;

[0041] Figure 6 The structure diagram of the well-to-well profile generation device in the embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0042] Exemplary embodiments of the present disclosure will be described in greater detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it is understood that the present disclosure can be embodied in various forms and should not be limited by 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 present disclosure to those skilled in the art.

[0043] In order to solve the problem that the interpolation algorithm in the prior art is heavily dependent on the spatial distribution of known data and is not applicable to well-to-well data interpolation of a well profile, an embodiment of the present disclosure provides a well-to-well data interpolation method for a well profile, which realizes one-to-one correspondence of well-to-well data by using a probability distribution function, and essentially solves the problem that the unevenness of known data distribution and the non-fixity of relative positions in the well profile adversely affect the interpolation result.

[0044] Embodiment one

[0045] Embodiment one of the present disclosure provides a well-to-well data interpolation method for a well profile, the flowchart of which is shown in Figure 1

[0046] Step S101: A probability density function is assigned to the attribute data to be interpolated of each well section in the well profile, and a probability distribution function corresponding to the probability density function of each well section is determined.

[0047] The attribute data to be interpolated of the well profile can be various characteristic attributes of each well section in the well profile, such as shale content, saturation, porosity, formation pressure, permeability, etc. The interpolation calculation can be performed for each attribute to be interpolated.

[0048] A probability density function greater than zero is assigned to the data corresponding to all well sections on the well profile to be interpolated, and a formula can be directly given according to needs, or an arbitrary continuous line segment can be hand-drawn as the probability density function, and the probability distribution function is calculated respectively.

[0049] The simplest density function is a constant, as long as the integral of the density function on the corresponding line segment is 1. For example, if the length of a well section is 2 meters, the density function can be set as a constant f(x) = 0.5, that is, 1 divided by the length 2.

[0050] The probability distribution function is the integral of the probability density function, that is, the area between the probability density function and the 0 axis. Using the above example, if the probability density function is f(x) = 0.5, the probability distribution function of the probability distribution function on the well trajectory corresponding line segment should be F(x) = 0.5x, and the range of x is 0 to 2, that is, to the end of the well section included in the layer to be interpolated, the probability distribution function should be 1.

[0051] ​Step S102: connecting the points with equal probability distribution function values on the adjacent well sections to obtain the data pair connection lines of the adjacent well sections.

[0052] Connecting the points with equal probability distribution function values on the adjacent well sections as the corresponding data pairs of the adjacent well sections, so that there is only one data pair connection line passing through any position on the to-be-interpolated continuous well profile.

[0053] Step S103: determining the to-be-interpolated attribute value of the to-be-interpolated point on the data pair connection line according to the to-be-interpolated attribute values of the corresponding points on the two adjacent well sections.

[0054] Based on the data pair connection line, the to-be-interpolated position point on the to-be-interpolated continuous well profile is subjected to an interpolation operation to obtain the to-be-interpolated attribute value of the to-be-interpolated position point, and then the interpolation data of any position on the entire profile is obtained. There are many interpolation operations, and any interpolation method can be selected here, such as a distance weighting method or any other method.

[0055] In the implementation of the above method, the key is how to find the equal probability distribution function connection line passing through the to-be-interpolated position between wells, and then determine the attribute value of the to-be-interpolated position between wells through the attribute values of the points on the two well sections at the ends of the equal probability distribution function connection line.

[0056] The above method of the embodiment realizes the interpolation of the attribute feature data between wells of the continuous well profile based on the probability density function and the probability distribution function, and introduces the probability density function and the probability distribution function to assist the profile interpolation based on the monotone non-decreasing feature of the probability distribution function in probability theory. The specific method is to set the probability density function of the known data (usually the logging curve) on the well in the same formation participating in the interpolation, realize the monotone increasing feature of the corresponding probability distribution function from zero to one, and thus ensure that any position between different well sections in the formation realizes a one-to-one corresponding relationship. Then, the values at the corresponding positions are used to perform linear or other way interpolation operation on any unknown position between wells, so as to realize the interpolation calculation and accurate adjustment of the related data of the entire continuous well profile.

[0057] Embodiment Two

[0058] The embodiment two of the present application provides a specific implementation process of the above-mentioned well-to-well profile data interpolation method, and the flowchart is as shown in Figure 2 , and the principle is as shown in Figure 3 , and the method comprises the following steps:

[0059] Step S201: determining the to-be-interpolated continuous well profile, horizon and to-be-interpolated attribute, obtaining the to-be-interpolated attribute data of the well sections included in the to-be-interpolated continuous well profile in the to-be-interpolated horizon range, and obtaining the to-be-interpolated attribute data of each well section in the continuous well profile.

[0060] The to-be-interpolated attribute data includes the position coordinates of each position on the well interval and the corresponding to-be-interpolated attribute value. The to-be-interpolated attribute data on each well interval on the to-be-interpolated horizon on the acquired cross-well profile is taken as initial known data. The to-be-interpolated attribute can be any attribute of the well interval, the cross-well profile is a two-dimensional plane, each position corresponds to a position coordinate (X, Y), the to-be-interpolated attribute is Z, and each position has a to-be-interpolated attribute Z value. However, only the Z value on the well track on the cross-well profile is known, and the Z value of other positions needs to be obtained through interpolation. The present application solves the problem of how to calculate the Z value of other positions. The known to-be-interpolated attribute data on the well interval is measured by logging and other means. Since the corresponding attribute value of each position on the well track can be measured, the position corresponding to the measured attribute value is also determined. Therefore, the position coordinates are also included in the above obtained known data, that is, the known data is attribute data including the position coordinates, that is, (X, Y, Z) are all known data.

[0061] Step S202: Preprocessing the to-be-interpolated attribute on each well interval in the cross-well profile: for the unknown position on the well interval that needs to be supplemented with the to-be-interpolated attribute value, based on the to-be-interpolated attribute values of the two adjacent positions, the to-be-interpolated attribute value of the unknown position is determined and supplemented.

[0062] The to-be-interpolated attribute data of each well interval is known. These data can be measured data, but the measured to-be-interpolated attribute data of each well interval cannot guarantee that it is continuous and uninterrupted. If the initial known data is non-continuous point data, the data is preprocessed by interpolation to make it continuous in the to-be-interpolated well interval. The continuity does not mean that the attribute values of all points on the well interval need to be listed or calculated by interpolation. In actual application, if the two adjacent points of the well interval have attribute values and the distance is not greater than a certain distance, the data can be considered to be continuous. If the distance is greater than a certain distance, the data can be supplemented by interpolation. In addition, a continuous straight line or curve of the to-be-measured attribute value can be fitted according to the measured point data. In addition, during the continuity processing, the attribute value of the position point required for interpolation calculation can be determined.

[0063] The continuity processing of the known data on the well track can ensure that the Z value of any position on the well track curve in the to-be-interpolated calculation area is known. If the Z value of some positions on the well track curve is unknown, the interpolation preprocessing can be performed by difference. The most commonly used method of interpolation preprocessing is to calculate the Z value of the unknown position point on the well track curve according to the Z values of the two nearest known position points above and below the unknown position point, and to weight the distance from the unknown position.

[0064] Step S203: Assign a probability density function to the interpolation attribute data of each well segment in the well profile, and determine the probability distribution function corresponding to the probability density function of each well segment.

[0065] For the interpolation attribute of well section k in the well profile, its probability density function is assigned as f. k (x k ), x k The value of f is within the length of the well section. The probability density function f for well section k... k (x k By performing integration, the probability distribution function F of the i-th well section is obtained. k (x k ).

[0066] like Figure 3 As shown, the two wells in the well-connected profile are Well 1 and Well 2, which together form a simple well-connected profile. For any interpolation point 3 in the well-connected profile, interpolation is performed. For the interpolated attribute, a probability density function f1(x1) is assigned to Well 1, and a probability density function f2(x2) is assigned to Well 2. Integral operations are performed on the probability density functions of Well 1 and Well 2 respectively to obtain the probability distribution function F1(x1) of Well 1 and the probability distribution function F2(x2) of Well 2.

[0067] Step S204: Connect the points with equal probability distribution function values ​​on adjacent well sections to obtain data pair connections between adjacent well sections.

[0068] For well sections k and s in the well profile, F k (x k ) = F s (x s Connect the points to obtain multiple data pairs between well segment k and well segment s.

[0069] Using the example above, based on the probability density function f1(x1) and probability distribution function F1(x1) of well 1, and the probability density function f2(x2) and probability distribution function F2(x2) of well 2, we determine the line 9 of the equal probability distribution function of the interpolation point 3, which is a line connecting a data pair through the interpolation point 3.

[0070] Step S205: Determine the interpolation attribute values ​​of the points to be interpolated on the data pair connection line based on the corresponding interpolation attribute values ​​of the data pair connection line on two adjacent well sections.

[0071] According to the distance a of the interpolation point to the well segment k, the distance b of the interpolation point to the well segment s, and the probability distribution value m corresponding to the interpolation point on the well segment k and the probability distribution value n corresponding to the interpolation point on the well segment s, the probability distribution value of the end point of the data pair connecting line is determined. The position of the end point is determined according to the probability distribution value of the end point, and the interpolation attribute value of the interpolation point is determined according to the position of the end point, the interpolation attribute value of the interpolation point, and the position of the interpolation point.

[0072] In the above example, the interpolation operation on the interpolation point 3 is implemented by using the attribute values corresponding to the position i and the position j at the two ends of the equal probability distribution function connecting line 9 passing through the interpolation point 3. There are many methods for the interpolation operation, for example, the following method can be used:

[0073] For the convenience of understanding, the auxiliary line 8 perpendicular to the well 1 and the well 2 is drawn through the interpolation point 3, and the distance a of the interpolation point 3 to the well 1 and the corresponding probability distribution function value m, the distance b of the interpolation point 3 to the well 2 and the corresponding probability distribution function value n can be conveniently obtained based on the auxiliary line 8, m is the probability distribution function value F1(c) corresponding to the foot position c of the perpendicular line of the interpolation point 3 to the well 1, and n is the probability distribution function value F2(d) corresponding to the foot position d of the perpendicular line of the interpolation point 3 to the well 2. According to the equal ratio of the corresponding sides of similar triangles, the following can be obtained:

[0074]

[0075] Since the probability distribution values corresponding to the two ends of the equal probability distribution function connecting line 9 are equal, i.e. F1(i)=F2(j).

[0076] From the above two formulas, the following can be obtained:

[0077] According to the inverse function F1 -1 (x1) of the probability distribution function F1(x1), the position i can be determined, and the position i is

[0078] According to the inverse function F2 -1 (x2) of the probability distribution function F2(x2), the position j can be determined, and the position j is

[0079] After the positions i and j are determined, the attribute value of the interpolation point 3 can be calculated by using the interpolation algorithm according to the attribute values of the positions i and j. The simplest algorithm is to weight according to the distance of the interpolation point 3 to the positions i and j. More complex algorithms can also be used to generate different interpolation attribute fields.

[0080] In some optional embodiments, the well-to-well data interpolation method provided in Embodiments 1 and 2 above, after obtaining the data pairs connecting adjacent well sections, can determine the rationality of the connection positions of each data pair based on the characteristics of the attributes to be interpolated and related geological concepts, and adjust the corresponding positions of the data pairs by adjusting the probability density function of the attribute data of each well section. Specific judgment methods can be selected as needed; examples are given below:

[0081] Example 1: If the attribute to be interpolated is clay content, the stacking pattern of the sand body is determined according to the sedimentary environment of the layer to be interpolated. If the stacking pattern of the sand body is channel sand, the probability density function is usually a flat-topped and convex-bottom shape. When setting the probability density function, it can be set to a convex-bottom curve. The actual value of the probability density at each position on the well section is 1 minus the probability density function value. The actual value of the probability density can correspond to the shape of the sand body. If the stacking pattern of the sand body is estuary bar sand, the probability density function is usually a flat-bottomed and convex-top shape. When setting the probability density function, it can be set to an convex-top curve.

[0082] Example 2: If the attribute to be interpolated is saturation, then the oil-water interface is usually used as the boundary. Below the oil-water interface is pure water, so the probability density function below the oil-water interface needs to be adjusted to be consistent. Above the oil-water interface, it is related to the sand body or physical property distribution characteristics, so the probability density function can be set and adjusted accordingly with reference to Example 1.

[0083] The principle for adjusting the probability density function is that the resulting attribute interpolation results conform to the geological understanding of geologists, or that the technical personnel believe, based on their experience in the field, that the attribute interpolation results have been adjusted to meet the requirements. Alternatively, an adjustment experience value or some experience data can be set, and the probability density function can be adjusted to make the attribute interpolation results conform to the corresponding experience data.

[0084] In some optional embodiments, the well-to-well data interpolation method for well profiles provided in Embodiments 1 and 2 above can, after determining the interpolation attribute values ​​of the points to be interpolated, judge the rationality of the obtained interpolation attribute values ​​based on experience. If there are obvious unreasonable aspects, the probability density function can be adjusted until the interpolation attribute values ​​of each point to be interpolated obtained after interpolation based on the adjusted probability density function are reasonable.

[0085] The well tie profile attribute interpolation method provided by the first and second embodiments above uses a probability distribution function to achieve one-to-one correspondence of data between wells, and essentially solves the problem that the uneven distribution of known data and the relative position of the known data in the well tie profile adversely affect the interpolation result. The interpolation result of the well tie profile will no longer change due to profile flattening of each layer and adjustment of horizontal and vertical scale, so that the interpolation algorithm has the characteristic of stability. The above method uses adjustment of the probability density function to achieve accurate control of the interpolation result and convenient display of geological understanding to the greatest extent, so that the interpolation algorithm has the characteristic of flexibility. Through the above characteristics, the two-dimensional geological analysis and prediction technology based on profile interpolation can be greatly optimized. Not only can it assist in deepening geological understanding, but also can be used as a two-dimensional model to verify a three-dimensional geological model. Whether as an analysis tool or a display tool, it can greatly improve the analysis depth and research efficiency of geological work.

[0086] Embodiment three

[0087] The embodiment three of the present application provides a well tie profile generation method, the flowchart of which is shown in Figure 4 , and includes the following steps:

[0088] Step S301: using the well tie profile interwell data interpolation method described above to determine the attribute values to be interpolated of each point to be interpolated between adjacent well sections in the well tie profile.

[0089] Step S302: generating a well tie profile based on the attribute values to be interpolated of each point to be interpolated.

[0090] After the attribute values of all points to be interpolated in the well tie profile are determined, the well tie profile can be easily generated according to these attribute values. Different interpolation calculations can be performed for different attributes to be interpolated to obtain well tie profile graphs of different attributes.

[0091] Based on the same inventive concept, the embodiment of the present application also provides a well tie profile interwell data interpolation device, which can be arranged in a device with interpolation calculation capability. The structure of the device is shown in Figure 5 , and includes a preprocessing module 11, a connecting line module 12 and an interpolation module 13.

[0092] The preprocessing module 11 is used to assign a probability density function to the attribute data to be interpolated of each well section in the well tie profile, and determine the probability distribution function corresponding to the probability density function of each well section.

[0093] The connecting line module 12 is used to connect the points with equal probability distribution function values on adjacent well sections to obtain the data pair connecting line of the adjacent well sections.

[0094] An interpolation module 13 is configured to determine the to-be-interpolated attribute value of the to-be-interpolated point on the data pair connection according to the to-be-interpolated attribute values of the data pair connection corresponding to the adjacent two well sections.

[0095] Optionally, the preprocessing module 11 is specifically configured to assign the probability density function f k (x k ) to the to-be-interpolated attribute of the well section k in the well-to-well profile, where x k is in the range of the length of the well section; and perform integral operation on the probability density function f k (x k ) of the well section k to obtain the probability distribution function F k (x k ) of the well section k.

[0096] Optionally, the connection module 12 is specifically configured to connect the points of F k (x k ) = F s (x s ) for the well section k and the well section s in the well-to-well profile to obtain a plurality of data pair connections between the well section k and the well section s.

[0097] Optionally, the interpolation module 13 is specifically configured to determine the probability distribution value of the end point on the data pair connection according to the distance a of the to-be-interpolated point to the well section k, the distance b of the to-be-interpolated point to the well section s, and the probability distribution value m corresponding to the to-be-interpolated point on the well section k and the probability distribution value n corresponding to the to-be-interpolated point on the well section s; determine the position of the end point according to the probability distribution value of the end point; and determine the to-be-interpolated attribute value of the to-be-interpolated point according to the position of the end point, the to-be-interpolated attribute value, and the position of the to-be-interpolated point.

[0098] In some optional embodiments, the preprocessing module 11 is further configured to determine the to-be-interpolated well-to-well profile, horizon, and to-be-interpolated attribute, obtain to-be-interpolated attribute data of the well section in the to-be-interpolated well-to-well profile within the horizon range, and obtain the to-be-interpolated attribute data of each well section in the well-to-well profile; the to-be-interpolated attribute data includes the position coordinates and the corresponding to-be-interpolated attribute value of each position on the well section.

[0099] Optionally, after obtaining the to-be-interpolated attribute data of each well section in the well-to-well profile, the preprocessing module 11 is further configured to preprocess the to-be-interpolated attribute data on each well section in the well-to-well profile: determine and supplement the to-be-interpolated attribute value of an unknown position on the well section that needs to supplement the to-be-interpolated attribute value based on the to-be-interpolated attribute values of the adjacent positions on both sides.

[0100] Based on the same inventive concept, the embodiments of the present application further provide a device for generating a well-to-well profile, which can be arranged in a device having interpolation calculation capability and well-to-well profile generation capability, and the structure of the device is as follows:Figure 6 As shown, comprising:

[0101] The well profile interwell data interpolation device 1 is used for determining the attribute values to be interpolated of each point to be interpolated between adjacent well sections in the well profile.

[0102] The generation module 14 is configured to generate a well profile map based on the attribute values to be interpolated of each point to be interpolated.

[0103] Based on the same inventive concept, the embodiments of the present application further provide a computer storage medium, which stores computer executable instructions, and the computer executable instructions are executed by a processor to implement the well profile interwell data interpolation method and / or the well profile map generation method.

[0104] Based on the same inventive concept, the embodiments of the present application further provide a device, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor executes the program to implement the well profile interwell data interpolation method and / or the well profile map generation method.

[0105] As to the method and device in the above embodiments, the related content has been described in one part and will not be described in detail in another part, and the related content in each part can be referred to each other.

[0106] The above method and device are suitable for interpolation of data distribution characteristics in the well profile, so as to solve the contradiction between the known data distribution density and data continuity, and the interpolation result is not affected by the flattening, changing of the vertical and horizontal scale and other operations of the well profile. The two-dimensional analysis means for the well profile interwell interpolation becomes a powerful tool for assisting geological understanding.

[0107] Unless specifically stated otherwise, terms such as processing, computing, calculating, determining, displaying, and the like, can refer to an action and / or process of one or more processing or computing systems, or similar devices, that manipulate and / or transform data represented as physical (e.g., electronic) quantities within the processing system's registers and / or memories into other data similarly represented as physical quantities within the processing system's memories, registers or other such information storage, transmission or display devices. Information and signals can be represented using any of a variety of different technologies and techniques. For example, data, instructions, commands, information, signals, bits, symbols, and chips that can be referenced throughout the above description can be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, or any combination thereof.

[0108] It should be understood that the particular order in which the steps of processes presented in the disclosure have been presented can be rearranged. Furthermore, various aspects of the disclosure can be used alone or in various combinations. In the appended claims, means-plus-function or step-plus-function clauses can have their application interpreted merely as open-ended limits of a particular combination of steps that are done on the apparatus, as specified in 35 U.S.C. § 112(f). However, the steps of the claims are not to be interpreted as requiring their performance in the order recited.

[0109] In the detailed description above, various features are grouped together in single embodiments for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting a necessity to more features than are expressly identified in each claim. Rather, inventive functionality can be more effectively disclosed by grouping features in a single claim. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate preferred embodiment.

[0110] Those skilled in the art will further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein can be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans can implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.

[0111] The steps of a method or algorithm described in connection with the embodiments disclosed herein can be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium can be integral to the processor. The processor and the storage medium can reside in an ASIC. The ASIC can reside in a user terminal. In the alternative, the processor and the storage medium can reside as discrete components in a user terminal.

[0112] For a software implementation, the techniques described herein can be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes can be stored in memory units and executed by processors. The memory unit can be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.

[0113] The above description includes examples of one or more embodiments. Of course, not all possible combinations of components or methods described above can be claimed as embodiments. One of ordinary skill in the art can recognize that modifications and variations of the described embodiments can be practiced that are still within the scope of the present application. Therefore, the described embodiments are intended to cover all such modifications and variations as are within the scope of the following claims. Further, the description herein is intended to enable practitioners of the art to practice the application as claimed. Also, the use of any term "means" or "machine" or "apparatus" or "device" or "component" or "step" or "element" or "step" or "member" or "block" or "module" or "circuit" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor" or "processor"

Claims

1. A method for inter-well data interpolation in well profiles, characterized in that, include: The well profile, stratigraphic level, and interpolation attribute to be interpolated are determined. The interpolation attribute data of the well segments included in the well profile to be interpolated within the stratigraphic level are obtained, resulting in the interpolation attribute data of each well segment in the well profile. The interpolation attribute data includes the position coordinates of each location on the well segment and the corresponding interpolation attribute value. Assign a probability density function to the attribute data to be interpolated for each well segment in the well profile, and perform integration on the probability density function of each well segment to obtain the probability distribution function corresponding to the probability density function of each well segment. Connecting points with equal probability distribution function values ​​on adjacent well segments yields multiple data pairs between adjacent well segments; the adjacent well segments are well segment k and well segment s. Based on the interpolation attribute values ​​corresponding to the data pair connection on two adjacent well segments, the interpolation attribute values ​​of the points to be interpolated on the data pair connection are determined, including: performing interpolation operations on the interpolation location points on the data pair connection of the well profile to be interpolated to obtain the interpolation attribute values ​​of the locations to be interpolated; determining the probability distribution value of the endpoint on the data pair connection based on the distance a from the interpolation point to well segment k, the distance b from the interpolation point to well segment s, and the probability distribution value m and the probability distribution value n corresponding to the interpolation point in well segment k and well segment s; determining the position of the endpoint based on the probability distribution value of the endpoint; and determining the interpolation attribute value of the interpolation point based on the position of the endpoint, the interpolation attribute value, and the position of the interpolation point.

2. The method as described in claim 1, characterized in that, Assigning a probability density function to the interpolated attribute data of each well segment in the well profile, and determining the probability distribution function corresponding to the probability density function of each well segment, includes: For the interpolation attribute of well section k in the well profile, the probability density function is assigned as follows: , The value range is the length of the well section; The probability density function for well section k By performing integration, the probability distribution function of well section k is obtained. .

3. The method as described in claim 2, characterized in that, The step of connecting points with equal probability distribution function values ​​on adjacent well sections to obtain data pair connections between adjacent well sections includes: For well sections k and s in the well profile, = Connect the points to obtain multiple data pairs between well segment k and well segment s.

4. The method according to any one of claims 1-3, characterized in that, After obtaining the interpolation attribute data for each well segment in the well profile, the following is also included: Preprocessing is performed on the interpolated attribute data of each well section in the well profile: for unknown locations in the well section where interpolated attribute values ​​need to be added, the interpolated attribute values ​​of the unknown locations are determined and added based on the interpolated attribute values ​​of the adjacent locations on both sides.

5. A method for generating a well-connected profile diagram, characterized in that, include: The interpolation attribute values ​​of each interpolation point between adjacent well sections in the well-connected profile are determined using the method described in any one of claims 1-4, and a well-connected profile map is generated based on the interpolation attribute values ​​of each interpolation point.

6. A well-to-well profile data interpolation device, characterized in that, include: The preprocessing module is used to determine the well profile, stratigraphic level, and interpolation attribute to be interpolated; to obtain the interpolation attribute data of the well segments included in the well profile within the stratigraphic level; and to obtain the interpolation attribute data of each well segment in the well profile. The interpolation attribute data includes the position coordinates of each location on the well segment and the corresponding interpolation attribute value. A probability density function is assigned to the interpolation attribute data of each well segment in the well profile, and the probability density function of each well segment is integrated to obtain the probability distribution function corresponding to the probability density function of each well segment. The connection module is used to connect points with equal probability distribution function values ​​on adjacent well segments to obtain multiple data pairs between adjacent well segments; the adjacent well segments are well segment k and well segment s; The interpolation module is used to determine the interpolation attribute value of the point to be interpolated on the data pair connection line based on the corresponding interpolation attribute values ​​on two adjacent well segments. This includes: performing an interpolation operation on the point to be interpolated on the data pair connection line to obtain the interpolation attribute value of the point to be interpolated; determining the probability distribution value of the endpoint on the data pair connection line based on the distance a from the point to be interpolated to well segment k, the distance b from the point to be interpolated to well segment s, and the probability distribution value m and n corresponding to the point to be interpolated on well segment k and well segment s; determining the position of the endpoint based on the probability distribution value of the endpoint; and determining the interpolation attribute value of the point to be interpolated based on the position of the endpoint, the interpolation attribute value, and the position of the point to be interpolated.

7. A device for generating well profile diagrams, characterized in that, include: The well-connected profile data interpolation device as described in claim 6 is used to determine the interpolation attribute values ​​of each interpolation point between adjacent well sections in the well-connected profile. The generation module is used to generate a well-connected profile based on the interpolation attribute values ​​of each interpolation point.

8. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the well-to-well data interpolation method of any one of claims 1-4 and / or the well-to-well profile generation method of claim 5.

9. A computer device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the well-to-well profile data interpolation method of any one of claims 1-4 and / or the well-to-well profile generation method of claim 5.

Images