A new method for establishing top structure map based on well-seismic combination

Abstract

A new method for establishing top surface structure map based on well-seismic combination belongs to the technical field of oil and gas exploration and development, and comprises the following steps: step 1, understanding the basic geological profile of the research area; step 2, pre-processing the seismic data; step 3, top surface horizon interpretation; step 4, generating horizon interpretation three-dimensional model; step 5, identifying the fault in the above three-dimensional model and inserting the fault rod to form the fault model to generate the three-dimensional fault curved surface model; step 6, based on the three-dimensional fault curved surface model, carrying out fault interpretation; step 7, after completing the horizon and fault interpretation, editing the structure map by using the drawing software. The application describes the representation method of the fault top surface structure map, which is based on the conventional structure interpretation process and the master of the structure characteristics of the research area, and different faults are formed by using the fault model function in the Petrel software. The representation method is simple and clear, and has strong operability and very strong practicability in the structure map compilation.

Description

Technical Field

[0001] This invention belongs to the field of oilfield seismic data interpretation technology, specifically relating to a new method for establishing top surface structural maps based on well-seismic integration. Background Technology

[0002] With the increasing demand for oil and gas energy, the exploration depth of oil and gas strata is gradually increasing, and Cretaceous oil and gas exploration has become a key focus in many regions. Based on current reservoir analysis, which primarily consists of structural oil and gas reservoirs, structural faults are the main controlling factor in hydrocarbon accumulation. Therefore, accurate and precise structural interpretation is the most fundamental and crucial task in oil and gas exploration.

[0003] The strata in the study area have undergone multiple tectonic movements. Faults mainly formed during the subsidence stage of the Songliao Basin tectonic movement. During this period, crustal movement was characterized by relatively rapid and stable subsidence, accompanied by intermittent uplift. The subsidence rate exceeded the deposition rate for a considerable period, resulting in large areas of deep lakes within the basin, forming sedimentary deposits dominated by non-compensating sedimentation, and depositing the main source and reservoir rock series. The depositional rate was relatively high during this period, ranging from 120 to 290 m / Ma for the Quantou Formation, decreasing to 30 to 50 m / Ma from the Qingshankou Formation to the Nenjiang Formation. Around the end of the Nenjiang period, the Sea of ​​Japan began to expand, and the westward compressional force affected the basin, known as the Nenjiang Movement. This created a compressive-torsional stress field, leading to folding and widespread uplift of the basin, more pronounced in the eastern region. Local tectonic and secondary tectonic zones formed, ending this evolutionary history.

[0004] In oil and gas exploration, top surface structural maps are fundamental maps. For faults that have no displacement at the top of the strata but play an important role in hydrocarbon accumulation, there is no clear method for representing them on structural maps. This has brought many inconveniences to oil and gas exploration work. To address this, we have invented a new method for representing top surface structural maps of faults, which solves the above-mentioned technical problems for top surface structural maps of strata formed by multiple tectonic movements.

[0005] Existing technology CN104849762B discloses a technique for calculating the true dip angle of formations in horizontal well development reservoirs, thereby correcting the microstructure of the reservoir top surface. Specifically, it involves utilizing the delayed response characteristics of natural gamma-ray logging curves crossing sandstone-mudstone interfaces to determine the length of the delay interval. Combined with the drilling inclination and azimuth angles interpreted during logging, a relatively accurate method for calculating the true dip angle of formations at the horizontal well entry point is derived and calculated. A formation dip angle plan view is then drawn, and based on this, the microstructure of the horizontal well reservoir top surface is corrected. Contour maps with spacing less than 2 meters are drawn to guide the rational development of the reservoir and the tapping of remaining oil potential. This invention overcomes the limitations of current methods that only utilize well point elevation and seismic data for top surface microstructure correction. It proposes a method for correcting the microstructure of the horizontal well reservoir top surface, fully utilizing the information from the horizontal well entry point to obtain the true formation dip angle, thus enabling more accurate research on the microstructure of the reservoir top surface.

[0006] The traditional methods used in the existing technology CN104849762B (such as using horizontal well entry point information to correct reservoir top surface microstructure) are difficult to accurately depict the geometry of faults and their relationship with surrounding strata in complex geological environments (such as normal fault areas). The lack of sufficient three-dimensional visualization means results in an inaccurate depiction of fault geometry. Summary of the Invention

[0007] To address the aforementioned problems, this invention proposes a novel method for establishing a top surface structural map based on well-seismic analysis, comprising the following steps:

[0008] Step 1: Understand the basic geological conditions of the study area;

[0009] Step 2: Preprocess the seismic data;

[0010] Step 3: Perform top-level stratigraphic interpretation;

[0011] Step 4: Generate a 3D model for layer interpretation;

[0012] Step 5: Identify faults in the above three-dimensional model and insert fault rods to form a fault model, generating a three-dimensional fault surface model;

[0013] Step 6: Based on the three-dimensional fault surface model, perform fault interpretation;

[0014] Step 7: After completing the interpretation of the stratigraphy and faults, edit the structural diagram using drawing software.

[0015] Furthermore, in step 1, the basic geological conditions of the study area are understood, relevant data on tectonic evolution are collected, and the tectonic movements and stress directions experienced by the study area are familiarized.

[0016] Furthermore, in step 2, for pre-stack or post-stack seismic data with poor signal-to-noise ratio, filtering and denoising processing is first performed on them, using structurally guided filtering or dip-guided filtering methods.

[0017] Furthermore, in step 2, the seismic data is subjected to coherent enhancement processing to obtain a coherent enhanced data volume. The coherent volume is calculated using the filtered seismic data. Based on the changes in tectonic amplitude and the degree of fault development in the study area, the calculation parameters are tested to obtain the coherent enhanced data volume.

[0018] Furthermore, in step 3, based on the well logging data from the drilled wells, the top interface is interpreted according to the conventional stratigraphic interpretation process. First, filtered seismic data is applied to find marker layers with strong seismic reflection amplitudes and relatively continuous phase axes near the target stratigraphic segment for structural interpretation, based on the well-seismic calibration results. Point-to-point or dip-based automatic tracking technology is used to interpret the stratigraphic position on the marker layer. For areas with local tracking anomalies, manual interactive modification is performed. SMI and GeoEast interpretation software are used to perform structural interpretation of the top surface of the stratigraphic layer. Specifically, the sonic transit time (AC) and density (DEN) logging curves are imported into the seismic interpretation software—Syntool module—to synthesize seismic records for stratigraphic calibration and stratigraphic interpretation.

[0019] Furthermore, in step 4, the interpreted top surface layer data is imported into Petrel, and the make surface module in Petrel is used to generate a three-dimensional model.

[0020] Furthermore, in step 5, the fault is interpreted to generate a fault model. The generated top-level stratigraphic model is displayed in three-dimensional space, and the area of ​​a certain fault to be interpreted is magnified. Next, the normal section perpendicular to the fault strike is opened to most clearly show the cross-sectional features of the fault. The fault is interpreted at intervals along the fault strike. After the interpretation is completed, linear interpolation technology is used to generate the fault model.

[0021] Further, in step 6, two sets of fault bars, namely the main survey line and the connecting line, are extracted based on the fault model. The control points on the fault bars are fitted to obtain spline fault bars, and spline mesh interpolation is performed on the two sets of spline fault bars to generate a three-dimensional fault surface model. The two sets of fault bars, namely the main survey line and the connecting line, are extracted from the generated fault model, and the control points on the fault bars are fitted using B-spline curves to obtain spline fault bars. The interval of the spline fault bars can be set equally according to the interpretation requirements. Then, spline mesh interpolation is performed on the two sets of spline fault bars to generate a three-dimensional fault surface model.

[0022] The equation of the B-spline curve used is:

[0023]

[0024] Among them, P i (i = 0, 1…n) are the control vertices (coordinates), B i,k (i = 0, 1, ..., n) are k-order normalized B-spline basis functions; the basis function is a sequence U of non-decreasing parameter u: u0 ≤ u1 ≤ ... ≤ u n+k+1 The determined k-th degree piecewise polynomial.

[0025] Furthermore, in step 6, based on the three-dimensional fault surface model, the spline fault bars are adjusted along any direction of the main survey line and the connecting line until the entire fault space is automatically closed for interpretation. On the three-dimensional surface model, spline fault bars in any direction of the main survey line and the connecting line are selected for adjustment. After adjustment, a new spline network surface model is generated based on the spline fault bars in that direction. The entire spline network surface model is updated in real time. This means that fault interpretation is performed by adjusting the three-dimensional surface instead of adjusting the conventional two-dimensional section fault bars, ensuring that the fault in this direction is closed for interpretation while also ensuring that the spline fault bars in the other direction are adjusted in real time. The projection points of the same fault on the same section and the spline fault bars are strictly closed.

[0026] Furthermore, in step 7, the faults are classified according to whether there is a fault displacement at the top interface, and then the faults are graded according to the size of the fault displacement and the fault length, and the faults are represented by different line types and line numbers.

[0027] The beneficial effects of this invention are as follows: The method for representing the top surface structure of a fault in this invention is mainly aimed at the representation method on the top interface structure map of a main normal fault. This representation method is based on understanding the structural characteristics of the study area and, according to the conventional structural interpretation process, utilizes the fault model function in Petrel software to form different faults. Moreover, this representation method is simple, clear, and highly operable, making it very practical in the compilation of structural maps. Attached Figure Description

[0028] Figure 1 This is a flowchart illustrating a specific embodiment of the method for representing a fault top surface structure diagram according to the present invention;

[0029] Figure 2 This is a comparison diagram of seismic profiles before constructing guided filtering in a specific embodiment of the present invention;

[0030] Figure 3 This is a comparison diagram of seismic profiles after constructing guided filtering in a specific embodiment of the present invention;

[0031] Figure 4 This is a schematic diagram illustrating the establishment of a top-level layer model in a specific embodiment of the present invention;

[0032] Figure 5This is a schematic diagram of establishing a three-dimensional fault surface model in a specific embodiment of the present invention;

[0033] Figure 6 This is a top surface structure diagram of a specific embodiment of the present invention. Detailed Implementation

[0034] Example 1

[0035] To make the technical means and objectives of this invention easier to understand, the invention is further described below with reference to specific embodiments. A new method for establishing a top surface structural map based on well-seismic integration includes the following steps:

[0036] Step 1: Understand the basic geological conditions of the study area;

[0037] Step 2: Preprocess the seismic data;

[0038] Step 3: Perform top-level stratigraphic interpretation;

[0039] Step 4: Generate a 3D model for layer interpretation;

[0040] Step 5: Identify faults in the above three-dimensional model and insert fault rods to form a fault model, generating a three-dimensional fault surface model;

[0041] Step 6: Based on the three-dimensional fault surface model, perform fault interpretation;

[0042] Step 7: After completing the interpretation of the stratigraphy and faults, edit the structural diagram using drawing software.

[0043] like Figure 1 As shown, the method for representing a fault top surface structure diagram according to the present invention includes:

[0044] Step 1 involves understanding the basic geological conditions of the study area and collecting relevant data on tectonic evolution, as well as familiarizing oneself with the tectonic movements and stress directions experienced by the area. The main purpose of this step is to gain an understanding of the basic geological framework of the study area, ensuring that the subsequent tectonic interpretation is consistent with the geological conditions of the region and making the tectonic interpretation more reasonable.

[0045] Step 2: Preprocess the seismic data to improve the overall quality of the seismic data.

[0046] Step 3: Interpret the top surface strata of the work area. Based on the logging data from the drilled wells, interpret the top interface according to the conventional stratigraphic interpretation process. This step is the foundation of this invention. It mainly involves using interpretation software such as SMI and GeoEast to perform structural interpretation of the top surface strata according to the conventional stratigraphic interpretation process. Specifically, the steps include importing AC (acoustic transit time), DEN (density), and other logging curves into the Syntool module of the seismic interpretation software to synthesize seismic records, calibrate the stratigraphy, and interpret the strata. During the strata interpretation process, it is necessary to check whether the interpreted strata are closed.

[0047] Step 4: Import the interpreted stratigraphic data into the Petrel software and use the make surface module to build a 3D model from the imported stratigraphic data.

[0048] Step 5: Identify faults in the aforementioned 3D model and insert fault rods to form a fault model. Based on the fault model, extract two sets of fault rods: the main survey line and the connecting line. Fit the control points on the fault rods to obtain spline fault rods, and perform spline mesh interpolation on the two sets of spline fault rods to generate a 3D fault surface model. Extract two sets of fault rods: the main survey line and the connecting line, from the generated fault model, and use B-spline curves to fit the control points on the fault rods to obtain spline fault rods. The interval of these spline fault rods can be set equally according to interpretation requirements. Then, perform spline mesh interpolation on the two sets of spline fault rods to generate a 3D fault surface model.

[0049] Step 6: Based on the three-dimensional fault surface model, select any direction of the main survey line and the connecting line to adjust the spline fault rod until the space of the entire fault is automatically closed and interpreted.

[0050] This step is the key to the invention. The result is shown in the diagram. Figure 5 As shown, fault interpretation is both a key point and a challenge. When interpreting faults, we must first pay attention to the inheritance of faults in the adjacent strata above and below. Secondly, we must refer to the tectonic evolution history of the study area to determine whether the fault is a normal fault or a reverse fault, whether it conforms to the tectonic stress direction of that period, and whether the fault cutting relationship is correct.

[0051] Step 7: After interpreting the stratigraphy and faults, the structural diagram needs to be edited using drawing software. The main task of editing the structural diagram is to edit the faults, which involves classifying and grading them. Faults are classified based on whether they have a displacement at the top interface, and then graded according to the size of the displacement and the fault length. Different line types and line numbers are then used to represent the faults.

[0052] Example 2

[0053] In a specific embodiment of the application of the present invention, such as Figure 1 As shown, Figure 1This is a flowchart of the accurate top surface mapping method of the present invention. The accurate top surface mapping method includes:

[0054] Step 101: Understand the basic geological conditions of the study area and collect relevant data on tectonic evolution, familiarizing yourself with the tectonic movements and stress directions experienced by the study area. The main purpose of this step is to gain an understanding of the basic geological framework of the study area, ensuring that the subsequent tectonic interpretation is consistent with the geological conditions of the region, and making the tectonic interpretation more reasonable.

[0055] Step 102: Perform coherent enhancement processing on the seismic data to obtain the coherent enhanced data volume.

[0056] 1. Preprocessing seismic data

[0057] This embodiment targets pre-stack or post-stack seismic data with poor signal-to-noise ratios. First, it performs filtering and noise reduction processing, which can be carried out using structurally guided filtering or dip-guided filtering techniques. The purpose of the filtering processing is to improve the quality of the seismic data, thereby highlighting the seismic profile characteristics of the faults. Figure 2 and Figure 3 This is a comparison profile of seismic activity before and after the construction of guided filtering in this embodiment. Figure 2 and Figure 3 It can be seen that the seismic data before filtering has a lot of noise, poor strata continuity within the fault zone, and unclear fault features. After filtering, the random noise in the seismic data is reduced, the signal-to-noise ratio is significantly improved, the fault points on the profile are clear, and the fault features are more obvious, which is conducive to the identification and interpretation of faults.

[0058] 2. Perform coherent enhancement processing on seismic data.

[0059] The coherence volume is calculated using filtered seismic data. In this embodiment, the coherence volume calculation uses a third-generation coherence algorithm, but other coherence algorithms commonly used in the field can also be employed. Based on the changes in tectonic amplitude and fault development in the study area, calculation parameters are tested to obtain a coherence-enhanced data volume, which highlights fault characteristics.

[0060] Step 103: Interpret the top layer of the work area.

[0061] Based on well logging data from drilled wells, the top interface is interpreted according to the conventional stratigraphic interpretation process. This step is the foundation of this invention. First, filtered seismic data is used to locate marker layers with strong seismic reflection amplitudes and relatively continuous phase axes near the target stratigraphic interval, based on well-seismic calibration results. Point-to-point or dip-based automatic tracking technology is used to interpret the stratigraphic level on the marker layers. For areas with local tracking anomalies, manual interactive modifications are performed. Structural interpretation of the top stratigraphic level is then performed using interpretation software such as SMI and GeoEast. Specifically, well logging curves such as AC (acoustic transit time) and DEN (density) are imported into the Syntool module of the seismic interpretation software to synthesize seismic records, calibrate the stratigraphy, and interpret the stratigraphic level. During the stratigraphic interpretation process, it is necessary to check whether the interpreted stratigraphic level is closed.

[0062] Step 104: Generate a 3D model for layer interpretation

[0063] Import the interpreted top-level data into Petrel, and use the make surface module in Petrel to generate a 3D model. Figure 4 This is a top-level model diagram of an embodiment of this disclosure.

[0064] Step 105: Interpret the fault to generate a fault model.

[0065] The generated top-level stratigraphic model is displayed in three-dimensional space, and the region of a specific fault to be interpreted is magnified. Next, a normal section perpendicular to the fault's strike is opened to clearly display its cross-sectional features. The fault is then interpreted at intervals along its strike. These intervals are not fixed; when the fault is relatively smooth and its strike variation is small, the interpretation interval can be larger. When the fault strike variation is significant, the interpretation interval can be finer near the point of strike change to better control the fault's morphology. After interpretation, linear interpolation techniques are used to generate the fault model.

[0066] Step 106: Based on the fault model, extract two sets of fault bars: the main survey line and the connecting line. Fit the control points on the fault bars to obtain spline fault bars, and perform spline mesh interpolation on the two sets of spline fault bars to generate a three-dimensional fault surface model. Extract the two sets of fault bars (main survey line and connecting line) from the generated fault model, and use B-spline curves to fit the control points on the fault bars to obtain spline fault bars. The interval of these spline fault bars can be set equally according to interpretation requirements. Then, perform spline mesh interpolation on the two sets of spline fault bars to generate a three-dimensional fault surface model. Figure 5 This is a three-dimensional tomographic surface model diagram of an embodiment of this disclosure.

[0067] The equation of the B-spline curve used in this example is:

[0068]

[0069] Among them, P i (i = 0, 1…n) are the control vertices (coordinates), B i,k (i = 0, 1, ..., n) are k-order normalized B-spline basis functions; the basis function is a sequence U of non-decreasing parameter u: u0 ≤ u1 ≤ ... ≤ u n+k+1 The determined k-th degree piecewise polynomial.

[0070] Step 107: Based on the three-dimensional fault surface model, select any direction of the main survey line and the connecting line to adjust the spline fault bars until the entire fault space is automatically closed for interpretation. On the three-dimensional surface model, select spline fault bars in any direction of the main survey line and the connecting line for adjustment. After adjustment, regenerate the spline network surface model based on the spline fault bars in that direction. The entire spline network surface model will be updated in real time. This means that fault interpretation is performed by adjusting the three-dimensional surface instead of adjusting the conventional two-dimensional section fault bars. This ensures that while the fault in this direction is closed for interpretation, the spline fault bars in the other direction are also adjusted in real time. The projection points of the same fault on the same section and the spline fault bars are strictly closed.

[0071] Step 108: After interpreting the stratigraphy and faults, the structural diagram needs to be edited using drawing software. The main task of editing the structural diagram is to edit the faults, which involves classifying and grading them. Faults are classified based on whether they have a displacement at the top interface, and then graded according to the size of the displacement and the fault length. Different line types and line numbers are then used to represent the faults. A specific example is provided. Figure 6 As shown, the process has ended.

[0072] Example 3

[0073] This invention aims to provide a mapping method for fault top surface structural maps, primarily addressing the shortcomings of existing technologies in the analysis of structural strata, especially normal faults. Normal faults are often formed due to crustal extension, and their fault planes and top surfaces are often difficult to represent accurately. This mapping method can more clearly depict the geometry of faults and their relationship with surrounding strata, especially in complex geological environments, effectively improving the accuracy of fault top surface structural maps.

[0074] The objective of this invention can be achieved through the following technical measures: a method for representing the top surface structure of a fault, comprising: Step 1, preprocessing seismic data to improve its quality; Step 2, interpreting the top surface stratigraphic level based on well logging data from drilled wells; Step 3, importing the interpreted stratigraphic level data into Petrel software, and using the make surface module in Petrel to generate stratigraphic levels; Step 4, identifying the fault in the aforementioned three-dimensional model and inserting fault rods to form a fault model. Based on the fault model, extracting two sets of fault rods—the main survey line and the connecting line—fitting control points on the fault rods to obtain spline fault rods, and performing spline mesh interpolation on the two sets of spline fault rods to generate a three-dimensional fault surface model; Step 5, based on the three-dimensional fault surface model, selecting any direction profile of the main survey line and the connecting line to adjust the spline fault rods until the entire fault space is automatically closed and interpreted.

[0075] The objective of this invention can also be achieved through the following technical measures:

[0076] The method for representing the top surface structure of a fault also includes, before step 1, obtaining basic geological information of the study area, collecting relevant data on tectonic evolution, and understanding the tectonic movements and stress directions experienced by the study area. In step 2, seismic records are synthesized based on well logging curves such as sonic transit time and density to calibrate the strata for top surface interpretation. During the strata interpretation process, it is necessary to check whether the interpreted strata are closed. In step 4, the generated strata model is displayed in three-dimensional space, and the area of ​​a fault to be interpreted is magnified; the normal profile perpendicular to the strike of the fault to be interpreted is opened, and the initial fault is interpreted at intervals along the strike of the fault to be interpreted. After the interpretation is completed, linear interpolation technology is used to generate the fault model. In step 5, two sets of fault bars, the main survey line and the connecting line, are extracted from the fault model; B-spline curves are used to fit the control points on the fault bars to obtain spline fault bars, and spline mesh interpolation is performed on the two sets of spline fault bars to generate the three-dimensional surface model.

[0077] The equation of the B-spline curve is:

[0078]

[0079] Among them, P i (i = 0, 1…n) are the control vertices (coordinates), B i,k (i = 0, 1, ..., n) are k-order normalized B-spline basis functions; the basis function is a sequence U of non-decreasing parameter u: u0 ≤ u1 ≤ ... ≤ u n+k+1 The determined k-th degree piecewise polynomial.

[0080] The above description is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any equivalent substitutions or modifications made by those skilled in the art within the scope of the technology disclosed in the present invention, based on the technical solution and concept of the present invention, should be covered within the scope of protection of the present invention.

Claims

1. A novel method for establishing a top surface structural map based on well-seismic integration, characterized in that, Includes the following steps: Step 1: Understand the basic geological conditions of the study area; Step 2: Preprocess the seismic data; Step 3: Perform top-level stratigraphic interpretation; Step 4: Generate a 3D model for layer interpretation; Step 5: Identify faults in the above three-dimensional model and insert fault rods to form a fault model, generating a three-dimensional fault surface model; Step 6: Based on the three-dimensional fault surface model, perform fault interpretation; Step 7: After completing the interpretation of the stratigraphy and faults, edit the structural diagram using drawing software.

2. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 1, characterized in that, In step 1, we learn about the basic geological conditions of the study area, collect relevant data on tectonic evolution, and become familiar with the tectonic movements and stress directions experienced by the study area.

3. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 2, characterized in that, In step 2, for pre-stack or post-stack seismic data with poor signal-to-noise ratio, filtering and denoising processing is first performed on them, using structural guided filtering or dip guided filtering methods.

4. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 3, characterized in that, In step 2, the seismic data is coherently enhanced to obtain a coherent enhanced data volume. The coherent volume is calculated using the filtered seismic data. Based on the changes in tectonic amplitude and the degree of fault development in the study area, the calculation parameters are tested to obtain the coherent enhanced data volume.

5. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 4, characterized in that, In step 3, based on the well logging data from the drilled wells, the top interface is interpreted according to the conventional stratigraphic interpretation process. First, filtered seismic data is applied to find marker layers with strong seismic reflection amplitudes and relatively continuous phase axes near the target stratigraphic segment for structural interpretation, based on the well-seismic calibration results. Point-to-point or dip-based automatic tracking technology is used to interpret the stratigraphic position on the marker layer. For areas with local tracking anomalies, manual interactive modification is performed. SMI and GeoEast interpretation software are used to perform structural interpretation of the top surface of the stratigraphic layer. Specifically, the sonic transit time (AC) and density (DEN) logging curves are imported into the seismic interpretation software—Syntool module—to synthesize seismic records, calibrate the stratigraphic layer, and interpret the stratigraphic position.

6. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 5, characterized in that, In step 4, the interpreted top surface layer data is imported into Petrel, and the make surface module in Petrel is used to generate a three-dimensional model.

7. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 6, characterized in that, In step 5, the fault is interpreted to generate a fault model. The generated top-level stratigraphic model is displayed in three-dimensional space, and the area of ​​a certain fault to be interpreted is magnified. Next, the normal section perpendicular to the fault strike is opened to most clearly show the cross-sectional features of the fault. The fault is interpreted at intervals along the fault strike. After the interpretation is completed, linear interpolation technology is used to generate the fault model.

8. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 7, characterized in that, In step 6, two sets of fault bars, namely the main survey line and the connecting line, are extracted based on the fault model. The control points on the fault bars are fitted to obtain spline fault bars, and spline mesh interpolation is performed on the two sets of spline fault bars to generate a three-dimensional fault surface model. The two sets of fault bars, namely the main survey line and the connecting line, are extracted from the generated fault model, and the control points on the fault bars are fitted using B-spline curves to obtain spline fault bars. The interval of the spline fault bars can be set equally according to the interpretation requirements. Then, spline mesh interpolation is performed on the two sets of spline fault bars to generate a three-dimensional fault surface model. The equation of the B-spline curve used is: Among them, P i (i = 0, 1…n) are the control vertices (coordinates), B i,k (i = 0, 1, ..., n) are k-order normalized B-spline basis functions; the basis function is a sequence U of non-decreasing parameter u: u0 ≤ u1 ≤ ... ≤ u n+k+1 The determined k-th degree piecewise polynomial.

9. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 8, characterized in that, In step 6, based on the three-dimensional fault surface model, the spline fault bars are adjusted along any direction of the main survey line and the connecting line until the entire fault space is automatically closed for interpretation. The spline fault bars in any direction of the main survey line and the connecting line are selected and adjusted on the three-dimensional surface model. After adjustment, a new spline network surface model is generated based on the spline fault bars in that direction. The entire spline network surface model is updated in real time. This means that fault interpretation is performed by adjusting the three-dimensional surface instead of adjusting the conventional two-dimensional section fault bars, ensuring that the fault in this direction is closed for interpretation while also ensuring that the spline fault bars in the other direction are adjusted in real time. The projection points of the same fault on the same section and the spline fault bars are strictly closed.

10. The novel method for establishing a top surface structural map based on well-seismic integration as described in claim 9, characterized in that, In step 7, the faults are classified according to whether there is a fault displacement at the top interface. Then, the faults are graded according to the size of the fault displacement and the fault length, and the faults are represented by different line types and line numbers.

Images