Multi-data fusion complex tectonic zone seismic interpretation and geological modeling method

Abstract

The invention belongs to the technical field of geologic structure interpretation and modeling, and particularly relates to a multi-data fusion seismic interpretation and geologic modeling method for a complex structural zone. The invention discloses a multi-data fusion seismic interpretation and geologic modeling method for a complex tectonic zone. The method comprises the following steps: observing and judging formation and fault development characteristics in a work area through an earth surface outcrop fault; interpreting a stratigraphic dip angle by adopting imaging logging to identify a fault; using 3D-VSP data to enhance seismic imaging and identify faults; a low-angle fault is recognized through conventional logging interpretation; identifying a velocity abrupt change interface by adopting a seismic velocity spectrum to determine a fault; through fault comprehensive judgment identified by different data and methods, the fault position is determined, and a geological structure mode of a research work area is constructed. According to the method, through fusion interpretation, fault position identification, seismic interpretation and geologic structure modeling of various data such as outcrop, logging data, seismic data and velocity field, the accuracy of geologic structure modeling can be improved, the success rate of oil-gas exploration is improved, and meanwhile, the cost is saved.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of geological structure interpretation and modeling, and particularly relates to a multi-data fusion complex structure zone seismic interpretation and geological modeling method. BACKGROUND

[0002] The geological structure model is an important basis for oil and gas exploration. In order to determine the geological structure characteristics, evaluate the oil and gas enrichment regularity, find the favorable factors of oil and gas enrichment, improve the success rate of exploration, reduce the risk, and obtain and establish an accurate underground structure model, it is crucial. The conventional geological structure modeling method basically depends on seismic data, mainly because the seismic data has wide coverage and large depth, and can well display the obvious stratum and fault interface, and generally seismic data acquisition is completed in the early exploration stage.

[0003] The existing conventional method is generally more suitable for medium and large basins with stable tectonic movement, long duration, and strong regularity of sedimentary process. Small basins are more easily affected by surrounding geological activities and sudden geological events due to their small lake basin range, and faults are more developed, which makes the stratum, structure, and lithology change faster and the characteristics more complex, thereby making it difficult to accurately interpret the seismic structure. At the same time, if the surface conditions are complex and the quality of the seismic acquisition data is not high, the structural interpretation will be more difficult. Even if new technologies are used for acquisition in the later exploration stage to improve the accuracy of structural interpretation, it will not only increase the cost, but also be difficult to improve the seismic quality under complex surface conditions.

[0004] At present, there are few methods for using seismic data to fuse multiple data to assist in geological structure interpretation and modeling. Patent CN201410213093.6 discloses a complex mountain front structure zone geological structure comprehensive modeling method, which uses field outcrop, gravity, magnetic method, electrical method, seismic and drilling multiple data to establish a preliminary model, including identifying and depicting faults, dividing the structure system according to the large fault, etc. Finally, the final model is obtained by continuously adjusting and modifying the forward simulation, effectively improving the rationality of modeling, but the use of logging data is not sufficient, only used to judge and interpret the stratum, obtain the velocity, and the method is more conventional, and the 3D-VSP data in the well is not considered. Patent CN202011055399.5 proposes a method of converting a seismic profile into a depth profile and marking the corresponding surface geological measured profile to obtain a preliminary structural modeling interpretation scheme, and finally modifying the initial scheme based on the structural deformation mode to obtain the final interpretation model. This method does not consider more data such as logging and velocity field; patent CN201711363156.6 proposes a method of establishing a three-dimensional space geological model grid according to geological sampling data, but this method is more inclined to optimize the three-dimensional operation and mapping method of the computer. SUMMARY

[0005] In view of the above problems, the purpose of the present application is to provide a multi-data fusion complex structure zone seismic interpretation and geological modeling method, which can improve the accuracy and rationality of geological structure modeling, reduce the error caused by using only single data and method, improve the accuracy of structure interpretation and modeling in oil and gas exploration work, and thus improve the success rate of oil and gas exploration while saving cost.

[0006] The technical scheme of the present application is as follows: a multi-data fusion complex structure zone seismic interpretation and geological modeling method, comprising the following steps: S1: observing the surface outcrop faults in the research area and collecting related stratum data, and judging the stratum and fault development characteristics in the research area; S2: interpreting the stratum dip angle by using imaging logging, and identifying the faults by using the stratum dip angle mutation points; S3: enhancing seismic imaging by using 3D-VSP data, and identifying the faults by using the seismic reflection characteristics in seismic imaging; S4: identifying the stratum repetition phenomenon by using conventional logging interpretation, and determining the low-angle faults; S5: identifying the velocity mutation interface by using seismic velocity spectrum, and determining the faults by using the velocity interface; S6: determining the fault position by comprehensively judging the faults identified by different data and methods in steps S1-S5, and constructing the geological structure model of the research area in combination with conventional seismic structure interpretation.

[0007] In the step S1, the surface outcrop faults in the research area are observed, and the related stratum data are collected, and the stratum and fault development characteristics in the research area are judged, and the specific process is as follows: S11: collecting the field surface outcrop data according to the basin distribution characteristics and geological background of the research area, judging the stratum distribution and the outcropped faults in the research area by using the surface outcrop characteristics, and determining the position, occurrence and footwall and hanging wall characteristics of the faults; S12: based on the stratum and fault characteristics of the surface outcrop, the corresponding position of the faults in the seismic profile is determined by extending to the underground based on the contrast of the seismic data.

[0008] In the step S2, the stratum dip angle is interpreted by using imaging logging, and the faults are identified by using the stratum dip angle mutation points, and the specific process is as follows: S21: interpreting the stratum dip angle and occurrence by using the imaging logging data of the drilled well, and obtaining the well depth-stratum dip angle variation graph; S22: identifying and counting the depth points of the stratum dip angle mutation in the well depth-stratum dip angle variation graph, and determining as the fault breakpoints; S23: using the well-seismic calibration method commonly used in seismic interpretation to calibrate the position of the breakpoint in the seismic profile, and finally determining the corresponding position of the fault in the seismic profile.

[0009] The 3D-VSP data is used in step S3 to enhance seismic imaging, and the fault is identified by using the seismic reflection characteristics in seismic imaging, and the specific process is as follows: S31: using 3D-VSP data to stack in the range partially coinciding with the seismic data of the research area, and improving the imaging effect of the data; S32: using the basic method and characteristics of fault identification in seismic data structural interpretation, including seismic reflection axis bending, faulting, phase conversion and amplitude change, to determine the position and depth of the fault and the breakpoint, and to determine the corresponding position of the fault in the seismic profile.

[0010] In step S4, one or more of the conventional logging interpretation including acoustic logging, resistivity logging, gamma logging, neutron logging or density logging is used, and the stratigraphic repetition is identified according to the characteristics of the conventional logging interpretation curve, including the similarity of the smooth, zoned and gradually changing rules of the logging curve, and the similarity of the curve shape. The specific method is to divide the zones according to the increasing and decreasing rules of the smooth logging curve, starting from the top of the curve to the first curve rule change point to divide a section, and then starting again to the next curve rule change point, and so on until the end of the curve. In all adjacent zones, the interval section with the same change rule is determined as the stratigraphic repetition section, and finally the depth point of the repetition position is determined to identify the low-angle fault.

[0011] In step S5, during the use of seismic data processing, the seismic velocity field is picked up from the seismic velocity spectrum, and the sudden interface of the velocity value in the velocity field is identified, and the seismic reflection axis of the sudden interface is taken as the fault to enhance the reliability of the identification of the position and shape of the fault in the region.

[0012] In step S6, the faults identified by different data and methods in steps S1-S5 are comprehensively judged to determine the position of the fault, and the conventional seismic structural interpretation is combined to build the geological structure model of the research area, which specifically includes: S61: performing logging stratigraphic interpretation and well-seismic calibration to obtain the seismic stratigraphic interpretation result; S62: using the basic method and characteristics of fault identification in seismic data, including one or more of seismic reflection axis bending, faulting, phase conversion and amplitude change, to identify the fault again on the basis of the fault interpretation result obtained in steps S1-S5, and to improve the interpretation result; S63: finally, the geological structure modeling is performed according to the seismic interpretation result to establish the geological structure model.

[0013] The technical effect of the present application is that the present application uses outcrops, logging data, seismic data, velocity field and other data to identify the position of the fault from multiple aspects, and uses 3D-VSP data to stack seismic data to locally improve the seismic imaging effect, which is suitable for the case that the surface condition is complex, the natural environment is harsh, the originally collected seismic data is old and of low quality, and it is difficult to collect again, the cost is high, and the quality is difficult to improve, starting from other data which is easy to collect and low in cost, and then performing seismic interpretation and geological structure modeling, so as to improve the accuracy and rationality of geological structure modeling to a certain extent, reduce the error caused by using single data and method, finally improve the exploration success rate, reduce the risk, and save the cost.

[0014] Further description will be made below in combination with the drawings. BRIEF DESCRIPTION OF DRAWINGS

[0015] Figure 1 A flow chart of a complex structure zone seismic interpretation and geological modeling method of multiple data fusion of the embodiment of the present application.

[0016] Figure 2 A surface fault outcrop map of a western faulted basin foothill zone of the embodiment of the present application.

[0017] Figure 3 A dip logging interpretation map of L114 well of the embodiment of the present application.

[0018] Figure 4 A dip logging interpretation map of L8 well of the embodiment of the present application.

[0019] Figure 5 A fault distribution form preliminary interpretation map according to the surface outcrop and the logging stratum dip of the embodiment of the present application.

[0020] Figure 6 A 3D-VSP data stacking enhanced seismic fault interpretation map of the embodiment of the present application.

[0021] Figure 7 A well-to-well stratum correlation map of the embodiment of the present application.

[0022] Figure 8 A fault interpretation result map of repeated strata according to well-to-well correlation of the embodiment of the present application.

[0023] Figure 9 A seismic stacking velocity field map of the embodiment of the present application.

[0024] Figure 10 A final seismic structure interpretation result map of the embodiment of the present application.

[0025] Figure 11 A final seismic structure interpretation result map of the embodiment of the present application. Figure 12 This is a geological structure model diagram of an embodiment of the present invention. Detailed Implementation Example 1

[0026] like Figure 1 As shown, a method for seismic interpretation and geological modeling of complex tectonic zones based on multi-data fusion includes the following steps: S1: Observe the exposed faults in the research area and collect relevant stratigraphic data to determine the strata and fault development characteristics in the research area; S2: Use imaging logging to interpret formation dip angles and use formation dip angle abrupt change points to identify faults; S3: Use 3D-VSP data to enhance seismic imaging and identify faults using seismic reflection features in seismic imaging; S4: Use conventional well logging interpretation to identify formation repetition phenomena and determine low-angle faults; S5: Use seismic velocity spectrum to identify velocity abrupt transition interfaces and use velocity interfaces to determine faults; S6: By comprehensively judging the faults identified by different data and methods in steps S1 to S5, the fault locations are determined, and a geological structure model of the study area is constructed in combination with conventional seismic tectonic interpretation.

[0027] In step S1, the exposed faults in the research area are observed and relevant stratigraphic data are collected to determine the strata and fault development characteristics of the research area. The specific process is as follows: S11: Based on the basin distribution characteristics and geological background of the research area, collect field surface outcrop data, use the surface outcrop characteristics to determine the stratigraphic distribution and faults exposed on the surface in the research area, and determine the location, attitude, and hanging wall and footwall characteristics of the faults. S12: Based on the stratigraphic and fault characteristics of the surface outcrops, and by referring to the seismic data, extend underground to determine the corresponding location of the fault in the seismic profile.

[0028] In step S2, imaging logging is used to interpret formation dip angles, and abrupt changes in formation dip angles are used to identify faults. The specific process is as follows: S21: Using drilling imaging logging data, interpret formation dip and attitude to obtain a well depth-formation dip variation diagram; S22: Identify and statistically analyze the depth points where the formation dip angle changes abruptly in the well depth-formation dip angle variation diagram, and determine them as fault breakpoints; S23: Use well-seismic calibration methods commonly used in seismic interpretation to calibrate the position of the fault in the seismic profile, and finally determine the corresponding position of the fault in the seismic profile.

[0029] The step S3 uses 3D-VSP data to enhance the seismic imaging, and uses the seismic reflection characteristics in the seismic imaging to identify the faults, and the specific process is as follows: S31: using 3D-VSP data, stacking in the range of local coincidence with the seismic data of the research area, improving the imaging effect of the data; S32: using the basic method and characteristics of the fault identification in the seismic data structure interpretation, including the seismic reflection axis bending, faulting, phase conversion and amplitude change, determining the fault and the breakpoint position and depth, and determining the corresponding position of the fault in the seismic profile.

[0030] The step S4 uses one or more of the conventional logging interpretation including acoustic logging, resistivity logging, gamma logging, neutron logging or density logging, and identifies the stratigraphic repetition according to the conventional logging interpretation curve characteristics, the conventional logging interpretation curve characteristics including the similarity of the smooth logging curve, the similarity of the zoning after the gradual change rule, the specific method is that the smooth logging curve is divided into a section from the top of the curve to the first curve rule change point, and then starts again to the next curve rule change point, and the same is true to the end of the curve, the interval section with the same change rule in all adjacent zones is determined as the stratigraphic repetition section, and finally the depth point of the repetition position is determined, and then the low-angle fault is identified.

[0031] The step S5 uses the seismic velocity field picked up from the seismic velocity spectrum in the process of using the seismic data processing, identifies the sudden interface of the velocity value in the velocity field, and uses the seismic reflection axis of the sudden interface as the fault, to enhance the reliability of the fault position and form identification in the region.

[0032] The step S6 determines the fault position by comprehensively judging the faults identified by different data and methods in steps S1-S5, and combines the conventional seismic structure interpretation to build the geological structure model of the research area, and the specific steps are as follows: S61: performing logging stratigraphic interpretation and well-seismic calibration to obtain the seismic stratigraphic interpretation result; S62: using the basic method and characteristics of the fault identification in the seismic data, including one or more of the seismic reflection axis bending, faulting, phase conversion and amplitude change, identifying the fault again on the basis of the fault interpretation result obtained in steps S1-S5, and perfecting the interpretation result; S63: finally, the geological structure modeling is performed according to the seismic interpretation result to establish the geological structure model.

[0033] The present application uses outcrops, well logging data, seismic data, velocity field and other various data to identify the fault position in multiple aspects, and simultaneously uses 3D-VSP data to locally improve the seismic imaging effect by stacking seismic data, which is suitable for the case that the surface condition is complex, the natural environment is harsh, the originally collected seismic data is old and of low quality, and it is difficult to collect seismic data again due to high cost and low quality. Starting from other various data which is easy to collect and of low cost, the present application fuses multiple data together, and then performs seismic interpretation and geological structure modeling, so as to improve the accuracy and rationality of the geological structure modeling to a certain extent, reduce the error caused by using single data and method, finally improve the exploration success rate and reduce the risk, and save the cost. Embodiment 2

[0034] The complex structure zone seismic interpretation and geological modeling method of multiple data fusion described in embodiment 1 is used to perform seismic interpretation and geological modeling on a complex structure zone in a western faulted basin piedmont, and the specific process is as follows: S1: The outcrop faults in the complex structure zone in the piedmont of a western faulted basin are observed and relevant stratum data are collected to determine the stratum and fault development characteristics in the working area; The surface condition of the piedmont zone in the block is complex and harsh, the altitude is high, the seismic data collection is difficult, the seismic data quality is not high due to the influence of the surface loose gravel layer, and the underground structure is complex and the large and small faults are very developed due to the influence of the orogenic movement. Through field geological investigation, the outcropped fault shape is collected, as shown in Figure 2 It can be found that the fault is nearly vertical with high angle, which indicates that there is one or more large-scale high-angle faults extending to the surface, which can roughly infer the development characteristics and shape of the fault in the block, as shown in Figure 5 ; S2: The stratum dip angle is interpreted by imaging logging, and the fault is identified by using the stratum dip angle mutation point; The stratum dip angle is identified by the imaging logging interpretation method, and the well depth-stratum dip angle change graph is obtained, as shown in Figure 3 , Figure 4 The dip angle change graph of L114 and L8 in Figure 3 It can be seen from that the stratum above 3670 meters in the L114 well is mainly east-dipping, the stratum dip angle below the stratum starts to change southward, and the stratum after 3990 meters is basically south-dipping. The L8 well has the same situation, and the stratum from 3920 meters to 3970 meters also changes from east-dipping to south-dipping. Then, the fault breakpoints of multiple wells in the adjacent area are identified, and two faults are basically determined. Then, the well-seismic calibration method commonly used in seismic interpretation is used to determine the corresponding position of the fault in the seismic profile, and the statistics are shown in Table 1. Thus, in combination with the nearly vertical and steep shape characteristics of the fault determined by the outcrop observation in step S1, the distribution shape of the two faults in the seismic profile can be roughly determined, as shown in Figure 5; Table 1 Breakpoint identification statistics of inclination information (partial wells)

[0035] S3: 3D-VSP data is used to enhance seismic imaging, and faults are identified by seismic reflection characteristics in seismic imaging; The 3D-VSP data in the area and the original seismic data are used for stacking, and the imaging effect is locally improved, as shown in Figure 6 It can be found that the stacked data is more easily identified for faults 1 and 2, and the faults are roughly determined in step two Figure 5 The rough fault distribution pattern is determined, and the faults are further described according to the basic method and characteristics of identifying faults in seismic data, so that the position and shape of the faults are more accurate; S4: Conventional well logging interpretation is used to identify stratigraphic repetition and determine low-angle faults; The conventional well logging data of adjacent wells are compared and analyzed, and the well logging curves are smoothed and processed in zones, and the stratigraphic repetition is identified by the stratigraphic change characteristics, and the low-angle faults are determined, as shown in Figure 7 Taking well L8 as an example, the resistivity curve is smoothed and divided into several increasing and decreasing regular change zones, Figure 8 It can be found that the change characteristics in the 3820m~3960m section and the 3960m~4100m section are basically the same, both showing a change characteristic of first rapid decrease and then slow increase, and the lengths of each change section are similar, which can be judged as stratigraphic repetition corresponding to Figure 7 The 2-3 layer is repeated, and according to the same curve rhythm characteristics, it can be judged that wells L114, Q2-23 and Q2-39 also have stratigraphic repetition, but the repeated stratigraphic thickness of wells Q2-23 and Q2-39 is smaller than that of wells L114 and L8. According to the stratigraphic repetition here, the positions of the fault breakpoints of wells L114, L8, Q2-23 and Q2-39 are determined at 3990m, 3970m, 3980m and 3950m. Using the basic method of well-seismic calibration in seismic interpretation, the fault can be described in the seismic section as a low-angle and flat fault, as shown in Figure 9 ; S5: Seismic velocity spectrum is used to identify velocity mutation interface, and the velocity interface is used to determine the fault; During the stacking process of seismic data, the seismic velocity field picked up from the velocity spectrum is used Figure 10 For example, the velocity value mutation interface is identified as the fault position, Figure 10 The position where the velocity field changes from light gray to dark gray and then to light gray can be used as two velocity mutation interfaces, and the determined fault position and shape are basically the same as the high-angle fault described in step three. This step can again enhance the reliability of the fault position and shape in this area; S6: The fault position is determined by the comprehensive judgment of the faults identified by different data and methods in steps S1-S5, and the conventional seismic structure interpretation is combined to build the geological structure model of the study area, and the specific process is as follows: The faults determined by different data in steps S1-S5 are superimposed, and the fault position is determined again by comprehensive interpretation, and the basic structure interpretation of seismic data is carried out, see Figure 11 , the dashed line in the figure is the stratum interpretation result determined by well logging stratum interpretation, well-seismic calibration with seismic data, and stratum interpretation combined with seismic reflection characteristics, the middle solid line is the main fault structure and distribution position of the region determined by steps S1-S5, and the upper and lower solid lines are the faults identified by seismic reflection characteristics, seismic reflection axis bending, faulting, phase conversion and amplitude change, the structure characteristics of which are consistent with the red faults, including high steep faults and low angle flat faults, finally, the geological structure modeling is carried out based on the interpretation result of Figure 10 , and the geological structure model diagram is determined, see Figure 11 .

[0036] According to the exploration practice of a western faulted basin, it is proved that the accurate and reasonable geological structure model of the piedmont zone area is basically established by using the present application, Figure 12 , the high angle nappe fault, fault 1 and fault 2, and the low angle slip fault, fault 3 and fault 4, two sets of fault systems are clearly defined, and the piedmont thrust zone lithology-structure reservoir is discovered, and the cumulative reported predicted oil geological reserves are ten million tons, and the small oilfield with an average annual output of nearly 100,000 tons is built.

[0037] The above is only the preferred specific embodiment of the present application, but the protection scope of the present application is not limited thereto, any change or replacement easily thought of by those skilled in the art within the technical range disclosed by the present application should be covered in the protection scope of the present application.

Claims

1. A method for multi-data fusion based complex structure zone seismic interpretation and geological modeling, characterized in that: The method comprises the following steps: S1: observing the surface outcrop faults in the research area and collecting related stratum data, judging the stratum and fault development characteristics in the research area; S2: interpreting the stratum dip angle by using imaging logging, and identifying the faults by using the stratum dip angle mutation points; S3: enhancing the seismic imaging by using 3D-VSP data, and identifying the faults by using the seismic reflection characteristics in the seismic imaging; S4: identifying the stratum repetition phenomenon by using the conventional logging interpretation, and determining the low-angle faults; S5: identifying the velocity mutation interface by using the seismic velocity spectrum, and determining the faults by using the velocity interface; S6: determining the fault positions by comprehensively judging the faults identified by different data and methods in steps S1-S5, and combining the conventional seismic structure interpretation to build the geological structure mode of the research area.

2. The method according to claim 1, wherein the method is characterized in that: In the step S1, the surface outcrop faults in the research area are observed and the related stratum data are collected to judge the stratum and fault development characteristics in the research area, and the specific process is as follows: S11: according to the basin distribution characteristics and geological background of the research area, collecting the field surface outcrop data, judging the stratum distribution and the outcropped faults in the research area by using the surface outcrop characteristics, and determining the position, occurrence, and upper and lower disc characteristics of the faults; S12: according to the stratum and fault characteristics of the surface outcrop, extending to the underground based on the contrast of the seismic data to determine the corresponding position of the faults in the seismic profile.

3. The method of claim 1, wherein: In the step S2, the stratum dip angle is interpreted by using the imaging logging, and the faults are identified by using the stratum dip angle mutation points, and the specific process is as follows: S21: interpreting the stratum dip angle and occurrence by using the imaging logging data of the drilled well to obtain the well depth-stratum dip angle change graph; S22: identifying and counting the depth points of the stratum dip angle mutation in the well depth-stratum dip angle change graph to determine the fault points; S23: using the well-seismic calibration method commonly used in seismic interpretation to calibrate the position of the points in the seismic profile, and finally determining the corresponding position of the faults in the seismic profile.

4. The method for complex structure zone seismic interpretation and geological modeling of multi-data fusion according to claim 1, characterized in that: In the step S3, the seismic imaging is enhanced by using the 3D-VSP data, and the faults are identified by using the seismic reflection characteristics in the seismic imaging, and the specific process is as follows: S31: using the 3D-VSP data to stack in the range locally coinciding with the seismic data of the research area to improve the imaging effect of the data; S32: using the basic method and characteristics of identifying the faults in the seismic data structure interpretation, including the seismic reflection axis bending, faulting, phase conversion and amplitude change, to determine the position and depth of the faults and points, and to determine the corresponding position of the faults in the seismic profile.

5. The method for complex structure zone seismic interpretation and geological modeling of multi-data fusion according to claim 1, characterized in that: The step S4 uses conventional well logging interpretation including one or more of acoustic wave, resistivity, gamma, neutron or density logging interpretation, and identifies formation repetition according to conventional well logging interpretation curve characteristics, including similarity of smooth, zoned, and regularity of the curve, and the specific method is to divide the curve into zones according to the regularity of the smooth curve, and to divide a zone from the top of the curve to the first curve regularity change point, and then to start again from the next curve regularity change point, and to proceed in the same way to the end of the curve, to determine the interval segment with the same change rule in all adjacent zones as the formation repetition segment, and to determine the depth point of the repetition position, and to identify low-angle faults.

6. The method for complex structure zone seismic interpretation and geological modeling of multi-data fusion according to claim 1, characterized in that: In the step S5, the seismic velocity field picked up from the seismic velocity spectrum in the seismic data processing process is used to identify the sudden interface of the velocity value in the velocity field, and the seismic reflection axis of the sudden interface is used as the fault to enhance the reliability of the identification of the fault position and the form in the region.

7. The method for complex structure zone seismic interpretation and geological modeling of multi-data fusion according to claim 1, characterized in that: In the step S6, the faults identified by different data and methods in the steps S1-S5 are comprehensively judged to determine the fault position, and the conventional seismic structure interpretation is combined to build the geological structure model of the study area, and the specific method includes: S61: well logging formation interpretation and well-to-seismic calibration are performed to obtain the seismic formation interpretation result; S62: the basic method and characteristics of seismic data fault identification are used, including one or more of seismic reflection axis bending, faulting, phase conversion and amplitude change, and the fault interpretation result obtained in the steps S1-S5 is used as the basis to identify the fault again to improve the interpretation result; S63: finally, the geological structure modeling is performed according to the seismic interpretation result to establish the geological structure model.

Images