High-precision three-dimensional geological modeling method for horizontal well area

Abstract

The invention provides a horizontal well area high-precision three-dimensional geological modeling method which comprises the steps that on the basis of manual interpretation of geological stratification, a horizontal section to be subjected to geological stratification treatment of a horizontal well in a target area is selected, and stratification treatment is conducted on the horizontal section of the horizontal well; based on the vertical well, the horizontal well and the seismic data of the target area, establishing a three-dimensional grid structure model of the target area; and establishing a three-dimensional geological attribute model of the target area. According to the method, the three-dimensional geological modeling effect can be improved, and the modeling precision is improved.

Description

Technical Field

[0001] This application relates to the field of shale oil and gas reservoirs, and in particular to a high-precision three-dimensional geological modeling method for horizontal well areas. Background Technology

[0002] In shale oil and gas reservoir development, horizontal wells are commonly used. Compared to conventional vertical wells, horizontal wells provide less information, such as formation thickness. Therefore, 3D geological modeling of shale reservoirs differs from conventional oil reservoir modeling. Currently, geological stratification interpretation of horizontal wells is typically done manually, comparing horizontal well characteristics with those of vertical wells and using the vertical well as a standard for stratification. 3D structural modeling employs a virtual well method, assuming stable shale reservoir depositional thickness with uniform lateral variation in 3D space. At each horizontal well stratification point, a virtual vertical well is assumed. Based on the thickness information of the smaller layers on the vertical well, other vertical geological stratification points are calculated at vertical depth. This method converts the depth information of the horizontal well into several vertical wells, making it applicable to conventional structural modeling methods. Then, in structural modeling, a method of hard-constraining seismic horizons at well stratification points is used to establish 3D structural horizons, ensuring consistency between the modeled horizons and the well stratification points while preserving the trend of the seismic horizons.

[0003] However, current methods for interpreting geological stratification in horizontal wells suffer from insufficient accuracy. This is primarily because as the thickness of the strata decreases, the characteristic curves of adjacent strata become more similar. In long horizontal sections, it becomes increasingly difficult to manually distinguish curve characteristics. Furthermore, the horizontal well trajectory more frequently cuts through various sub-strata, making it challenging to manually determine the relationship between strata encountered by the drill bit. This lack of accuracy results in larger intervals between adjacent geological stratification points in horizontal well sections, which is insufficient to constrain seismic horizons, leading to incorrect strata crossings in 3D modeling. Simultaneously, insufficient accuracy in interpreting stratum thickness leads to insufficient vertical accuracy in 3D geological attribute modeling. Summary of the Invention

[0004] This application provides a high-precision three-dimensional geological modeling method for horizontal well areas, aiming to solve the problem of insufficient vertical accuracy in existing three-dimensional geological attribute modeling.

[0005] The technical solution of this application is:

[0006] A high-precision three-dimensional geological modeling method for horizontal well areas includes the following steps:

[0007] S1, Based on the artificial interpretation of geological stratification, select the horizontal segment of the horizontal well in the target area to be geologically stratified, and perform stratification processing on the horizontal segment of the horizontal well;

[0008] S2, based on vertical wells, horizontal wells and seismic data of the target area, establish a three-dimensional mesh construction model of the target area;

[0009] S3, Establish a three-dimensional geological attribute model of the target area.

[0010] As a technical solution of this application, the automatic layering process for the horizontal segment in step S1 includes the following steps:

[0011] S11, in the conventional well logging curve sequence, select or fit a layered curve to characterize different formation characteristics;

[0012] S12, Automatic picking of layered curve feature points: Select a vertical well section on the vertical well that corresponds to the horizontal section. Identify multiple first inflection point peaks of the layered curve of the vertical well section by the magnitude of the deviation of the layered curve of the vertical well section from the mean line, and record the multiple first inflection point peaks as the first feature points of the layered curve of the vertical well section. Identify multiple second inflection point peaks of the layered curve of the horizontal section on the horizontal well by the magnitude of the deviation of the layered curve of the horizontal section from the mean line, and record the multiple second inflection point peaks as the second feature points of the layered curve of the horizontal section.

[0013] S13, perform layered calibration on multiple second feature points picked up on the horizontal segment using full path search and minimum error method;

[0014] S14, set a longitudinal accuracy value; when the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is greater than the set longitudinal accuracy value, repeat the operations in steps S12 to S13 on the horizontal section until the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is less than the set longitudinal accuracy value; map the two adjacent automatically stratified second feature points in step S13 to points A and B on the corresponding formations on the vertical well section respectively; if the longitudinal length of the vertical well section with points A and B as the top and bottom is greater than the set longitudinal accuracy value, continue to repeat the operations in steps S12 to S13 on the horizontal section determined with the two adjacent automatically stratified second feature points as the top and bottom until the longitudinal length of the vertical well section with points A and B as the top and bottom is less than the set longitudinal accuracy value.

[0015] As a technical solution of this application, step S13, performing layered calibration on the multiple second feature points picked on the horizontal segment includes:

[0016] S131, preset a percentage error σ;

[0017] S132, the second feature point in the horizontal section is denoted as T(T1, ..., Tm), and the first feature point in the vertical well section is denoted as S(S1, ..., Sn);

[0018] S133, the simulated drill bit starts from point T1 at the starting point and drills forward sequentially in the horizontal section. The feature value of the second feature point is compared with the feature value of a corresponding first feature point on the vertical well: if the feature value of the second feature point is the same as the feature value of a first feature point on the vertical well or the difference percentage is less than σ, then the second feature point is marked as a possible point; if the feature value of the second feature point, after being multiplied by the coefficients 1+σ and 1-σ respectively, is within the feature value range of two adjacent first feature points on the vertical well, then the second feature point is marked as a possible point between formation sections.

[0019] S134, based on the marking status of each of the second feature points, the simulated drill bit searches forward sequentially; when the simulated drill bit can reach Tm and Tm can be marked as Sn, it indicates that the path searched by the simulated drill bit is feasible, and the percentage error σ1 between the second feature point of the horizontal segment and the corresponding layer curve value at the first feature point in the vertical well segment is calculated; after searching all feasible paths, the path with the smallest percentage error σ1 is taken as the optimal path, and the automatic layer calibration result on the horizontal segment is obtained.

[0020] As a technical solution of this application, in step S2, establishing the three-dimensional mesh construction model of the target area includes:

[0021] S21, Establish a virtual vertical well: Create a virtual vertical well at the planar position of the horizontal well stratification point, and the trajectory of the virtual vertical well passes through the target shale reservoir;

[0022] S22, Calculate the thickness of the sub-layers on each of the virtual vertical wells;

[0023] S23, Correction of seismic interpretation depth domain horizon;

[0024] S24, Calculate the depth level of each sub-layer of each of the virtual vertical wells;

[0025] S25 uses vertical proportional subdivision and orthogonal regular meshing methods to establish a three-dimensional mesh construction model of the target area.

[0026] As a technical solution of this application, in step S22, the dip angle of the layer is calculated based on the seismic interpretation depth domain layer plane, and the dip angle of the formation of the vertical well is extracted from it, and the apparent thickness of the small layer at the vertical well is converted into the true thickness of the small layer;

[0027] If the number of vertical wells is equal to 1, then the true thickness of each sub-layer on the virtual vertical well is equal to the true thickness of the sub-layer of the vertical well;

[0028] If the number of vertical wells in the target area is greater than 1, the true thickness spatial variation surface of each sub-layer is calculated by using the Kriging interpolation method on the plane, and the true thickness of each sub-layer on the virtual vertical well is extracted from it.

[0029] The formation dip angle of the virtual vertical well is extracted by the dip angle surface of the formation, and the true thickness of each sub-layer of the virtual vertical well is converted into apparent thickness and marked as D.

[0030] As a technical solution of this application, in step S23, the layer points on the vertical well are sequentially marked as B1, ..., Bk, ..., Bp; based on the apparent thickness D value of each sublayer of each virtual vertical well, the layer point Bk corresponding to the virtual vertical well is moved to a set unified sublayer interface; the first depth average value of the corresponding layer points on each virtual vertical well after the movement is obtained, and the second depth average value of the corresponding layer points on each horizontal well is obtained on the seismic depth domain layer plane; the difference between the first depth average value and the second depth average value is corrected to the seismic interpretation depth domain layer plane, thus completing the correction of the horizontal well depth and seismic depth information; using the cokriging method, the seismic layer plane is constrained by the point depth value of the moved virtual vertical well, and a well-seismic constrained layer plane is established and marked as surface Y.

[0031] As a technical solution of this application, in step S24, the depth and stratum plane of each layer of the virtual vertical well are calculated based on the apparent thickness D of each layer of the virtual vertical well and the combined surface Y.

[0032] As a technical solution of this application, in step S3, based on the three-dimensional mesh construction model, well logging information and seismic data are combined, and the cokriging interpolation method is used to establish a three-dimensional geological attribute model of the target area.

[0033] The beneficial effects of this application are:

[0034] This application provides a high-precision three-dimensional geological modeling method for horizontal well areas. It densifies the number of well-logging geological stratification points in space by automatically and quantitatively identifying well-logging geological stratification points, and further constrains the three-dimensional structural modeling using stratification points, thereby improving the spatial accuracy of the structural model. As a result, when carrying out geological attribute modeling, it can more precisely fill the geological attribute model with external input attribute information such as well logging and seismic data, thereby improving the three-dimensional geological modeling effect and increasing the accuracy of three-dimensional geological attribute modeling. Attached Figure Description

[0035] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained from these drawings without creative effort.

[0036] Figure 1 A flowchart illustrating the high-precision three-dimensional geological modeling method for horizontal well areas provided in this application embodiment;

[0037] Figure 2 This is a schematic diagram of automatic picking of the first feature point in a vertical well layered structure, provided in an embodiment of this application.

[0038] Figure 3 This is a schematic diagram of automatic picking of second feature points for horizontal well stratification provided in an embodiment of this application;

[0039] Figure 4 A cross-sectional view of a well structure model established using automatic stratification data of a horizontal well, provided in an embodiment of this application;

[0040] Figure 5 A three-dimensional geological model cross-section provided for an embodiment of this application. Detailed Implementation

[0041] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can typically be arranged and designed in various different configurations.

[0042] Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0043] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0044] In the description of this application, it should be noted that the terms "upper" and "lower" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings, or the orientation or positional relationship in which the product of the invention is usually placed when in use. They are only used to facilitate the description of this application and to simplify the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this application.

[0045] Furthermore, in this application, unless otherwise expressly specified and limited, "above or below" the first feature may include direct contact between the first and second features, or contact between the first and second features through another feature between them. Moreover, "above," "over," and "on" the first feature includes the first feature directly above or diagonally above the second feature, or simply indicates that the first feature is at a higher horizontal level than the second feature. "Below," "below," and "under" the first feature includes the first feature directly below or diagonally below the second feature, or simply indicates that the first feature is at a lower horizontal level than the second feature.

[0046] Furthermore, terms such as "horizontal" and "vertical" do not imply that components must be absolutely horizontal or suspended, but rather that they can be slightly tilted. For example, "horizontal" simply means that its direction is more horizontal than "vertical," and does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0047] In the description of this application, it should also be noted that, unless otherwise expressly specified and limited, the terms "set up," "connected," and "linked" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal connection of two components. Those skilled in the art can understand the specific meaning of the above terms in this application based on the specific circumstances.

[0048] Example:

[0049] Please refer to Figure 1 (Refer to) Figures 2 to 5 This application provides a high-precision three-dimensional geological modeling method for horizontal well areas, including the following steps:

[0050] S1. Based on the artificial interpretation of geological stratification, select the horizontal section of the horizontal well in the target area that needs to be geologically stratified, and perform stratification processing on the horizontal section of the horizontal well.

[0051] S2, based on vertical wells, horizontal wells and seismic data of the target area, establish a three-dimensional mesh construction model of the target area;

[0052] S3, Establish a three-dimensional geological attribute model of the target area.

[0053] In step S1, the automatic layering of horizontal segments includes the following steps:

[0054] S11, in the conventional well logging curve sequence, select or fit a layered curve to characterize different formation characteristics;

[0055] S12, Automatic Picking of Layered Curve Feature Points: In a vertical well, select a vertical well segment corresponding to the horizontal segment. Identify multiple first inflection point peaks of the layered curve of the vertical well segment by the magnitude of the deviation of the layered curve from the mean line, and record these multiple first inflection point peaks as the first feature points of the layered curve of the vertical well segment. Identify multiple second inflection point peaks of the layered curve of the horizontal segment by the magnitude of the deviation of the layered curve from the mean line, and record these multiple second inflection point peaks as the second feature points of the layered curve of the horizontal segment.

[0056] S13, the multiple second feature points picked on the horizontal segment are hierarchically labeled by full path search and minimum error method;

[0057] S14, Set a longitudinal accuracy value; when the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is greater than the set longitudinal accuracy value, repeat the operations in steps S12 to S13 on the horizontal section until the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is less than the set longitudinal accuracy value; map the two adjacent automatically stratified second feature points in step S13 to points A and B on the corresponding formations on the vertical well section respectively; if the longitudinal length of the vertical well section with points A and B as the top and bottom is greater than the set longitudinal accuracy value, continue to repeat the operations in steps S12 to S13 on the horizontal section determined with the two adjacent automatically stratified second feature points as the top and bottom until the longitudinal length of the vertical well section with points A and B as the top and bottom is less than the set longitudinal accuracy value.

[0058] It should be noted that, in step S13, the layered labeling of multiple second feature points picked on the horizontal segment includes:

[0059] S131, preset a percentage error σ;

[0060] S132, the second feature point in the horizontal section is denoted as T(T1, ..., Tm), and the first feature point in the vertical section is denoted as S(S1, ..., Sn);

[0061] S133, the simulated drill bit starts from point T1 at the starting point in the horizontal section and drills forward sequentially. The characteristic value of the second characteristic point is compared with the characteristic value of a corresponding first characteristic point on the vertical well: if the characteristic value of the second characteristic point is the same as the characteristic value of a first characteristic point on the vertical well, or the percentage difference is less than σ, then the second characteristic point is marked as a possible point; if the characteristic value of the second characteristic point, after being multiplied by the coefficients 1+σ and 1-σ respectively, is within the characteristic value range of two adjacent first characteristic points on the vertical well, then the second characteristic point is marked as a possible point between formation sections.

[0062] S134, based on the marking status of each second feature point, the simulated drill bit searches forward sequentially; when the simulated drill bit can reach Tm and Tm can be marked as Sn, it indicates that the path searched by the simulated drill bit is feasible, and the percentage error σ1 between the second feature point of the horizontal section and the corresponding layer curve value at the first feature point in the vertical section is calculated; after searching all feasible paths, the path with the smallest percentage error σ1 is taken as the optimal path, and the automatic layer calibration result on the horizontal section is obtained.

[0063] It should be noted that in step S2, establishing the three-dimensional mesh construction model of the target area includes:

[0064] S21, Establish a virtual vertical well: Create a virtual vertical well at the planar location of the horizontal well stratification point, and the trajectory of the virtual vertical well passes through the target shale reservoir;

[0065] S22, calculate the thickness of the sub-layers on each virtual vertical well;

[0066] S23, Correction of seismic interpretation depth domain horizon;

[0067] S24, calculate the depth and stratigraphic plane of each sub-layer in each virtual vertical well;

[0068] S25 uses vertical proportional subdivision and orthogonal regular meshing methods to establish a three-dimensional mesh construction model of the target area.

[0069] It should be noted that in step S22, the dip angle of the layer is calculated based on the seismic interpretation depth domain layer plane, and the dip angle of the formation of the vertical well is extracted from it. The apparent thickness of the small layer at the vertical well is converted into the true thickness of the small layer.

[0070] If the number of vertical wells is equal to 1, then the true thickness of each sub-layer on the virtual vertical well is equal to the true thickness of the sub-layer of the vertical well.

[0071] If the number of vertical wells in the target area is greater than 1, the true thickness spatial variation surface of each sub-layer is calculated by using the Kriging interpolation method on the plane, and the true thickness of each sub-layer on the virtual vertical well is extracted from it.

[0072] The formation dip angle of the virtual vertical well is extracted by the dip angle surface of the formation, and the true thickness of each sub-layer of the virtual vertical well is converted into the apparent thickness and marked as D.

[0073] It should be noted that in step S23, the layer points on the vertical well are sequentially marked as B1, ..., Bk, ..., Bp; based on the apparent thickness D value of each sublayer of each virtual vertical well, the layer point Bk corresponding to the virtual vertical well is moved to a set unified sublayer interface; the first depth average value of the corresponding layer points on each virtual vertical well after the movement is obtained, and the second depth average value of the corresponding layer points on each horizontal well is obtained on the seismic depth domain layer plane; the difference between the first depth average value and the second depth average value is corrected to the seismic interpretation depth domain layer plane, completing the correction of the horizontal well depth and seismic depth information; using the cokriging method, the seismic layer plane is constrained by the point depth value of the moved virtual vertical well, establishing a well-seismic constrained layer plane, and marked as surface Y.

[0074] It should be noted that in step S24, the depth and stratum planes of each sub-layer of the virtual vertical well are calculated based on the apparent thickness D of each sub-layer of the virtual vertical well and the combined surface Y.

[0075] It should be noted that in step S3, based on the three-dimensional grid construction model, well logging information and seismic data are combined, and the cokriging interpolation method is used to establish a three-dimensional geological attribute model of the target area.

[0076] It should be noted that, Figure 2 It is a vertical well used to depict the stratification. Three geological stratification points, 1, 2, and 3, are manually marked on the map. Two stratification feature points, 1.1 and 1.2, are also identified by the high-precision three-dimensional geological modeling method for horizontal well areas provided in this embodiment. Figure 3 It is a horizontal well awaiting automatic stratification. Figure 3 Two geological stratification markers, 1 and 2, were manually marked in the horizontal well. Using the high-precision 3D geological modeling method for the horizontal well area provided in this embodiment, nine feature points were identified. A total of 144 feasible paths were obtained using the minimum error method and full path search. Among them, four stratification points, 1.1, 1.1, 1.1, and 1.2, were identified as the optimal paths. The average error between the optimal path stratification points and the vertical well stratification points was 5.6927%. The vertical accuracy was set to 10 meters, meaning only one loop was run. Therefore, this method can effectively improve the efficiency and accuracy of 3D geological modeling.

[0077] In summary, this application provides a high-precision three-dimensional geological modeling method for horizontal well areas. It densifies the number of well-logging geological stratification points in space by automatically and quantitatively identifying well-logging geological stratification points, and further constrains the three-dimensional structural modeling using these stratification points, thereby improving the spatial accuracy of the structural model. As a result, when carrying out geological attribute modeling, it can more precisely fill the geological attribute model with external input attribute information such as well logging and seismic data, improving the three-dimensional geological modeling effect and increasing the accuracy of three-dimensional geological attribute modeling.

[0078] The above description is merely a preferred embodiment of this application and is not intended to limit this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A high-precision three-dimensional geological modeling method for horizontal well areas, characterized in that, Includes the following steps: S1, Based on the artificial interpretation of geological stratification, select the horizontal segment of the horizontal well in the target area to be geologically stratified, and perform stratification processing on the horizontal segment of the horizontal well; S2, based on vertical wells, horizontal wells and seismic data of the target area, establish a three-dimensional mesh construction model of the target area; S3, Establish a three-dimensional geological attribute model of the target area.

2. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 1, characterized in that, In step S1, the automatic stratification of the horizontal segment includes the following steps: S11, in the conventional well logging curve sequence, select or fit a layered curve to characterize different formation characteristics; S12, Automatic picking of layered curve feature points: Select a vertical well section on the vertical well that corresponds to the horizontal section. Identify multiple first inflection point peaks of the layered curve of the vertical well section by the magnitude of the deviation of the layered curve of the vertical well section from the mean line, and record the multiple first inflection point peaks as the first feature points of the layered curve of the vertical well section. Identify multiple second inflection point peaks of the layered curve of the horizontal section on the horizontal well by the magnitude of the deviation of the layered curve of the horizontal section from the mean line, and record the multiple second inflection point peaks as the second feature points of the layered curve of the horizontal section. S13, perform layered calibration on multiple second feature points picked up on the horizontal segment using full path search and minimum error method; S14, set a longitudinal accuracy value; when the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is greater than the set longitudinal accuracy value, repeat the operations in steps S12 to S13 on the horizontal section until the longitudinal length between the top and bottom of the horizontal section and the top and bottom of the vertical well section is less than the set longitudinal accuracy value; map the two adjacent automatically stratified second feature points in step S13 to points A and B on the corresponding formations on the vertical well section respectively; if the longitudinal length of the vertical well section with points A and B as the top and bottom is greater than the set longitudinal accuracy value, continue to repeat the operations in steps S12 to S13 on the horizontal section determined with the two adjacent automatically stratified second feature points as the top and bottom until the longitudinal length of the vertical well section with points A and B as the top and bottom is less than the set longitudinal accuracy value.

3. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 2, characterized in that, In step S13, the layered calibration of the multiple second feature points picked on the horizontal segment includes: S131, preset a percentage error σ; S132, the second feature point in the horizontal section is denoted as T(T1, ..., Tm), and the first feature point in the vertical well section is denoted as S(S1, ..., Sn); S133, the simulated drill bit starts from point T1 at the starting point and drills forward sequentially in the horizontal section. The feature value of the second feature point is compared with the feature value of a corresponding first feature point on the vertical well: if the feature value of the second feature point is the same as the feature value of a first feature point on the vertical well or the difference percentage is less than σ, then the second feature point is marked as a possible point; if the feature value of the second feature point, after being multiplied by the coefficients 1+σ and 1-σ respectively, is within the feature value range of two adjacent first feature points on the vertical well, then the second feature point is marked as a possible point between formation sections. S134, based on the marking status of each of the second feature points, the simulated drill bit searches forward sequentially; when the simulated drill bit can reach Tm and Tm can be marked as Sn, it indicates that the path searched by the simulated drill bit is feasible, and the percentage error σ1 between the second feature point of the horizontal segment and the corresponding layer curve value at the first feature point in the vertical well segment is calculated; after searching all feasible paths, the path with the smallest percentage error σ1 is taken as the optimal path, and the automatic layer calibration result on the horizontal segment is obtained.

4. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 1, characterized in that, In step S2, establishing the three-dimensional mesh construction model of the target area includes: S21, Establish a virtual vertical well: Create a virtual vertical well at the planar position of the horizontal well stratification point, and the trajectory of the virtual vertical well passes through the target shale reservoir; S22, Calculate the thickness of the sub-layers on each of the virtual vertical wells; S23, Correction of seismic interpretation depth domain horizon; S24, Calculate the depth level of each sub-layer of each of the virtual vertical wells; S25 uses vertical proportional subdivision and orthogonal regular meshing methods to establish a three-dimensional mesh construction model of the target area.

5. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 4, characterized in that, In step S22, based on the seismic interpretation depth domain horizon, the horizon dip surface is calculated, and the formation dip angle of the vertical well is extracted from it. The apparent thickness of the sub-layer at the vertical well is converted into the true thickness of the sub-layer. If the number of vertical wells is equal to 1, then the true thickness of each sub-layer on the virtual vertical well is equal to the true thickness of the sub-layer of the vertical well; If the number of vertical wells in the target area is greater than 1, the true thickness spatial variation surface of each sub-layer is calculated by using the Kriging interpolation method on the plane, and the true thickness of each sub-layer on the virtual vertical well is extracted from it. The formation dip angle of the virtual vertical well is extracted by the dip angle surface of the formation, and the true thickness of each sub-layer of the virtual vertical well is converted into apparent thickness and marked as D.

6. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 5, characterized in that, In step S23, the layered points on the vertical well are sequentially marked as B1, ..., Bk, ..., Bp; based on the apparent thickness D value of each sublayer of each virtual vertical well, the layered point Bk corresponding to the virtual vertical well is moved to a set unified sublayer interface; the first average depth of the corresponding layered points on each virtual vertical well after the movement is calculated, and the second average depth of the corresponding layered points on each horizontal well is calculated on the seismic depth domain horizon; the difference between the first average depth and the second average depth is corrected to the seismic interpretation depth domain horizon, completing the correction of the horizontal well depth and seismic depth information; using the cokriging method, the seismic horizon is constrained by the point depth values ​​of the moved virtual vertical well, establishing a well-seismic constrained horizon, and marked as surface Y.

7. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 6, characterized in that, In step S24, based on the apparent thickness D of each sub-layer of the virtual vertical well and the combined surface Y, the depth horizon of each sub-layer of the virtual vertical well is calculated.

8. The high-precision three-dimensional geological modeling method for horizontal well areas according to claim 1, characterized in that, In step S3, based on the three-dimensional grid construction model, well logging information and seismic data are combined, and the cokriging interpolation method is used to establish a three-dimensional geological attribute model of the target area.

Images