Stratum pinching point and pinching line identification method and system

Through the forward simulation and hierarchical interpretation of the stratigraphic tip point recognition method, combined with seismic data and logging curve, the problem of large error in the stratigraphic tip point recognition is solved, the accuracy and efficiency of identification are improved, and more accurate is achieved identification of stratigraphic tip extinction lines.

CN119105106APending Publication Date: 2024-12-10CHINA PETROLEUM & CHEMICAL CORP +1
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202310608070.4
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2023-05-26
Publication Date
2024-12-10

AI Technical Summary

Technical Problem

The existing stratigraphic tip point recognition methods have large errors when predicting stratigraphic tip line, especially in areas where the stratigraphic inclination changes rapidly, which are greatly affected by factors such as seismic data resolution and signal-to-noise ratio.

Method used

By conducting main frequency analysis of seismic data in the identified area, combining the logging curve and formation inclination of the drilled wells, a geological forward model is established, and a forward simulation is used for different phase wavelets, and a phase wavelet with the smallest error is selected for sharp extinction points. Identify and perform hierarchical interpretation and extrapolation on the optimal phase seismic data body to reduce identification errors.

Benefits of technology

The accuracy of the recognition of stratigraphic tip point is improved, the recognition error is reduced, and the reliability of the recognition results is improved through the control of simulation efficiency and the acquisition of reflection coefficients.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119105106A_ABST
    Figure CN119105106A_ABST
Patent Text Reader

Abstract

The invention belongs to the technical field of oil-gas exploration, and particularly relates to a stratum pinching point and pinching line identification method and system. According to the method, a corresponding geological forward modeling model is established through geophysical parameters of a to-be-identified area, a stratigraphic dip angle and a set stratigraphic vanishing point position, forward modeling of the geological condition of the to-be-identified area is simulated, and the vanishing point identification error of different-phase wavelets with the same frequency is judged; horizon interpretation is carried out on the basis of the phase seismic data volume corresponding to the phase wavelet with the minimum error, so that the error existing in pinout point identification can be reduced as far as possible; moreover, horizon interpretation is carried out on the basis of the phase seismic data volume corresponding to the phase wavelet with the minimum error, and according to the error of the phase wavelet with the minimum error obtained by forward modeling, the identified position of the vanishing point after horizon interpretation is extrapolated so as to compensate the error, so that the error is further eliminated, and the accuracy of the method is improved. And the accuracy of vanishing point identification is improved.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The invention belongs to the technical field of oil and gas exploration, and in particular relates to a formation pinch-out point and pinch-out line identification method and system. Background Art

[0002] Stratigraphic traps are relatively common in sandstone oil and gas reservoirs, but how to accurately depict the stratigraphic pinch-out line is the main difficulty in identifying the stratigraphic pinch-out boundary. The technologies in existing patents and non-patent literature mainly use methods such as well-seismic combination and seismic phase classification to identify stratigraphic pinch-out points; for example, the Chinese invention patent with application publication number CN110794477A discloses a method for identifying pinch-out lines of overlapping stratigraphic formations by combining well-seismic combination. This method uses the structural maps of the two marker layers above and below the layer of the identified pinch-out line to make a difference, and sets a thickness threshold based on the thickness of the drilled formation. This method is effective for formations with stable formation dips and gradual pinch-outs, but is not applicable to areas where the formation dips vary greatly within a small range and the formation thickness changes rapidly; the Chinese invention patent with application publication number CN110794477A discloses a method for identifying pinch-out lines of overlapping stratigraphic formations by combining well-seismic combination. The method uses the structural maps of the two marker layers above and below the layer of the identified pinch-out line to make a difference, and sets a thickness threshold based on the thickness of the drilled formation. This method is effective for formations with stable formation dips and gradual pinch-outs, but is not applicable to areas where the formation dips vary greatly within a small range and the formation thickness changes rapidly; Chinese invention patent number CN113805227A discloses a method for identifying hidden formation pinch-out points, which mainly identifies formation pinch-out lines by using instantaneous frequency attributes along the layer, according to the mutation points of the attributes, combined with drilling and waveform differences; Chinese invention patent number CN112394398A discloses a method for identifying formation pinch-out points based on lithofacies and seismic facies, which mainly uses a geological forward model to compare the seismic profile obtained by seismic simulation, analyzes the changes in seismic amplitude and waveform caused by formation thinning and pinch-out, and artificially identifies formation pinch-out lines in interpretation;

[0003] However, all of the above-mentioned methods for identifying formation pinch-out points are based on post-stack pure wave seismic data. Due to many seismic factors such as the resolution and signal-to-noise ratio of seismic data and differences in sedimentary diversity in different regions, there is a large error between the predicted formation pinch-out line and the actual one. The author Yue Xiwei mentioned a method for identifying formation pinch-out lines using amplitude attributes in the "Method for Identifying Formation Pinch-Out Points Using Tuned Body Attributes" published in "Special Oil and Gas Reservoirs" in February 2022. Although this method can reduce errors compared to the identification method based on original seismic data, its accuracy is affected by the effective frequency band range and the main frequency, and there may still be a large error between the depicted pinch-out line and the actual pinch-out line. Summary of the invention

[0004] The object of the present invention is to provide a formation pinch-out point and pinch-out line identification method and system, which are used to solve the problem that there is a large error between the identification result of the existing formation pinch-out point identification method and the actual pinch-out line.

[0005] In order to achieve the above object, the present invention provides a formation pinch-out point identification method, comprising the following steps:

[0006] 1) Perform a main frequency analysis on the seismic data of the area to be identified to determine the main frequency of earthquakes in the area to be identified;

[0007] 2) Obtaining the well logging curves of the wells that have been drilled corresponding to the area to be identified, and determining the geophysical parameters of the area to be identified;

[0008] 3) Determine the formation dip of the area to be identified;

[0009] 4) setting the position of the stratigraphic pinch-out point of the geological forward model, and establishing the geological forward model by combining the geophysical parameters and stratigraphic dip of the area to be identified obtained in steps 2) and 3);

[0010] 5) According to the earthquake main frequency of the area to be identified obtained in step 1) and the geological forward model established in step 4), forward simulation is performed using different phase wavelets within a set phase angle range to obtain the error between the simulated formation pinch-out point corresponding to each different phase wavelet and the set formation pinch-out point of the geological forward model;

[0011] 6) Selecting the phase wavelet with the smallest error as the optimal phase wavelet, obtaining the optimal phase seismic data volume according to the optimal phase wavelet, and identifying the formation pinch-out point according to the optimal phase seismic data volume.

[0012] The beneficial effect of the above technical solution is: by simulating the forward modeling of the geological conditions of the area to be identified, the error in the identification of the pinch-out point caused by different phase wavelets of the same frequency is determined, and the stratigraphic interpretation is performed based on the phase seismic data body corresponding to the optimal phase wavelet with the smallest error, thereby minimizing the error in the identification of the pinch-out point.

[0013] Furthermore, the method for identifying the formation pinch-out point according to the optimal phase seismic data volume is:

[0014] A stratigraphic interpretation is performed on the optimal phase seismic data volume, and the position of the pinch-out point identified after the stratigraphic interpretation is extrapolated, the extrapolated value being the error value between the simulated stratigraphic pinch-out point corresponding to the optimal phase wavelet and the set stratigraphic pinch-out point of the geological forward model, and the new position after extrapolation is used as the identification result of the stratigraphic pinch-out point in the area to be identified.

[0015] The beneficial effects of the above technical solution are: not only is the stratigraphic interpretation performed on the basis of the phase seismic data volume corresponding to the phase wavelet with the smallest error, but also the error in the identification of the vanishing point by the phase wavelet with the smallest error obtained by forward simulation is used to extrapolate the position of the vanishing point identified after the stratigraphic interpretation to compensate for the error, thereby further eliminating the error and improving the accuracy of the identification of the vanishing point.

[0016] Furthermore, different phase wavelets within the set phase angle range are used for forward simulation, and the error between the simulated formation pinch-out point corresponding to each different phase wavelet and the set formation pinch-out point of the geological forward model is obtained in the following manner:

[0017] Within the set phase angle range, at each interval of the set phase angle, the corresponding phase wavelet is taken once for forward simulation to obtain the simulated formation pinch-out points corresponding to each different phase wavelet, and the differences between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points are respectively used as the errors between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points of the geological forward model; the frequencies of the different phase wavelets are all the main frequencies of earthquakes in the area to be identified.

[0018] The beneficial effects of the above technical solution are: while ensuring the comprehensiveness of the search for the optimal phase wavelet within the set phase angle range, the number of forward simulations required to search for the optimal phase wavelet is controlled, thereby improving the efficiency of determining the optimal phase wavelet.

[0019] Furthermore, the method of obtaining the optimal phase seismic data volume according to the optimal phase wavelet is: obtaining the optimal phase seismic data volume according to the convolution of the optimal phase wavelet and the reflection coefficient.

[0020] Furthermore, the reflection coefficient is extracted from the post-stack pure-wave seismic data volume.

[0021] The beneficial effects of the above technical solution are: the process of obtaining the reflection coefficient is simple, and the obtained reflection coefficient is more consistent with the actual situation.

[0022] Furthermore, the logging curves of the drilled wells mainly include acoustic wave curves, density curves, natural gamma curves and neutron curves.

[0023] Furthermore, there are at least two drilled wells, and the method for determining the formation inclination of the area to be identified is: obtaining the burial depth and actual horizontal well spacing of the drilled wells corresponding to the area to be identified, and determining the geological layer passed by the drilled wells according to the burial depth of the drilled wells; according to the actual horizontal well spacing of each drilled well and the altitude difference of each drilled well corresponding to the same geological layer, calculate the formation inclination of the area to be identified.

[0024] The present invention also provides a formation pinch-out point identification system, comprising a processor, wherein the processor is used to execute program instructions to implement the formation pinch-out point identification method.

[0025] The formation pinch-out point identification system can achieve the same beneficial effects as the formation pinch-out point identification method described above.

[0026] The present invention also provides a formation pinch-out line identification method, which determines the plane position of the formation pinch-out line according to the formation pinch-out point identification result, and the formation pinch-out point identification adopts the formation pinch-out point identification method as described above.

[0027] The formation pinch-out line identification method can achieve the same beneficial effects as the formation pinch-out point identification method described above.

[0028] The present invention also provides a formation pinch-out line identification system, comprising a processor, wherein the processor is used to execute program instructions to implement the formation pinch-out line identification method mentioned above.

[0029] The formation pinch-out line identification system can achieve the same beneficial effects as the formation pinch-out line identification method described above. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 It is a flowchart of a formation pinch-out point identification method in an embodiment of the formation pinch-out point identification method of the present invention;

[0031] Figure 2 It is a schematic diagram of the error comparison between the simulated formation pinch-out points and the actual formation pinch-out points corresponding to different phase wavelets in the embodiment of the formation pinch-out point identification method of the present invention;

[0032] Figure 3 A schematic diagram of a formation pinch-out point determined by a negative 90-degree phase seismic data volume in an embodiment of a formation pinch-out point identification method of the present invention;

[0033] Figure 4 It is a schematic diagram of the plane position of the formation pinch-out line determined by superimposing the recognition results of multiple areas to be recognized in the embodiment of the formation pinch-out point recognition method of the present invention. DETAILED DESCRIPTION

[0034] In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention is further described in detail below with reference to the accompanying drawings and embodiments.

[0035] Embodiment of the method for identifying formation pinch-out points

[0036] This embodiment provides a formation pinch-out point identification method, referring to Figure 1 , including the following steps:

[0037] 1) Performing a main frequency analysis on the seismic data of the area to be identified to determine the main frequency of the seismic in the area to be identified; the main frequency of the seismic in the area to be identified is the main frequency of the target layer of the seismic data in the area to be identified;

[0038] This example takes the Upper Paleozoic A group in a certain area as an example. There are 46 completed wells (drilled wells), all of which penetrate the Paleozoic strata. The actual drilling shows that there is a missing stratum in the A group of the study area (i.e., the area to be identified), and there is a pinch-out line. The study area covers an area of ​​1000km 2 , the entire area is fully covered by three-dimensional seismic data.

[0039] The main frequency analysis of the seismic data in the study area showed that the main frequency of the seismic target layer in the study area is 30HZ.

[0040] 2) Obtaining the well logging curves of the wells that have been drilled corresponding to the area to be identified, and determining the geophysical parameters of the area to be identified;

[0041] In this embodiment, the well logging curves of the drilled wells corresponding to the area to be identified mainly include acoustic wave curves, density curves, natural gamma curves and neutron curves, that is, the acoustic wave, density and neutron logging curves of 46 drilled wells are sorted out to obtain geophysical parameters such as velocity and wave impedance at each well point; since the drilled wells have been logged, the inverse of the acoustic wave time difference obtained during acoustic logging (i.e., acoustic wave time difference logging) is the velocity at the well point, which is actually the velocity of the longitudinal wave propagation of the formation from the surface to the bottom of the well at the drilled well; the wave impedance is the product of the velocity and the density.

[0042] 3) Determine the formation dip of the area to be identified;

[0043] Specifically, there are at least two drilled wells, and the method for determining the formation dip of the area to be identified is: obtain the burial depth and actual horizontal well spacing of the drilled wells corresponding to the area to be identified, and determine the geological layer passed by the drilled wells according to the burial depth of the drilled wells; calculate the formation dip of the area to be identified according to the actual horizontal well spacing of each drilled well and the altitude difference of each drilled well corresponding to the same geological layer; in this embodiment, the formation dip is determined to be 4 degrees according to the burial depth of the bottom of group A of 46 drilled wells and the actual well spacing.

[0044] 4) Setting the position of the stratigraphic pinch-out point of the geological forward model, and combining the geophysical parameters and stratigraphic dip of the area to be identified obtained in steps 2) and 3), establishing a geological forward model; in this embodiment, the geological modeling work is carried out using Tesserial Pro forward simulation software, and the established geological forward model is a wedge-shaped geological forward model; inputting different phase sub-waves and main frequencies into the forward simulation software, using a self-excited and self-collected simulation method, running the geological model, and then forward simulation can be performed.

[0045] 5) According to the earthquake main frequency of the area to be identified obtained in step 1) and the geological forward model established in step 4), forward simulation is performed using different phase wavelets within a set phase angle range to obtain the error between the simulated formation pinch-out point corresponding to each different phase wavelet and the set formation pinch-out point of the geological forward model;

[0046] Specifically, within the set phase angle range, at each interval of the set phase angle, a corresponding phase wavelet is taken for forward simulation to obtain the simulated formation pinch-out points corresponding to each different phase wavelet, and the difference between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points is used as the error between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points of the geological forward model; wherein the frequencies of the different phase wavelets are all the main frequencies of earthquakes in the area to be identified;

[0047] In this embodiment, the phase angle range is set to be from negative 90 degrees to positive 90 degrees, and the main frequency of the earthquake in the area to be identified is 30 Hz, that is, phase wavelets with a frequency of 30 Hz from negative 90 degrees to positive 90 degrees are used respectively, with an interval of 5 degrees, to carry out forward simulation, and finally 37 simulated formation pinch-out points are obtained, which are subtracted from the set formation pinch-out point position to calculate the error data between each simulated formation pinch-out point and the actual formation pinch-out point (that is, the set formation pinch-out point position). This is equivalent to first setting an actual formation pinch-out point (that is, setting the formation pinch-out point position) when modeling a geological forward model, and then using different phase wavelets with a frequency of the main frequency of the earthquake for forward modeling based on the geological forward model. Since the geological forward model is established according to the actual seismic data related to the area to be identified, the simulated formation pinch-out points obtained by the forward modeling are compared with the actual formation pinch-out points to obtain the error size of pinch-out identification of each phase wavelet under the geological conditions of the area to be identified.

[0048] In this embodiment, the actual formation pinch-out point (i.e., the set formation pinch-out point position) set during modeling is at 0 meters. Figure 2 , the top part is a simplified topographic map of the area to be identified corresponding to the modeling, from which the shape of the stratum and the pinch-out position can be seen; the lower part shows the forward modeling results of each phase wavelet (due to the large amount of data, only the forward modeling results of the zero-phase wavelet, the 90-degree phase wavelet and the negative 90-degree phase wavelet are shown), where the wave group in the dotted area represents the response of the group A stratum. As the thickness of the group A stratum decreases, the wave group inside the dotted line gradually disappears. In the prior art, the place where the wave group disappears in the three-dimensional seismic data is considered to be the pinch-out point of the stratum; the results show that when the zero-phase wavelet is used, the simulated stratum pinch-out point is at 200 meters, compared with the designed 0 meters, the error is 200 meters, and the error value is the largest; when the 90-degree phase wavelet is used, the simulated stratum pinch-out point is at 150 meters, with an error of 150 meters; when the negative 90-degree phase wavelet is used, the simulated stratum pinch-out point is at 130 meters, with an error of 130 meters, and the error value is the smallest.

[0049] 6) Select the phase wavelet with the smallest error as the optimal phase wavelet, obtain the optimal phase seismic data volume based on the optimal phase wavelet, and identify the formation pinch-out point based on the optimal phase seismic data volume to determine the position of the formation pinch-out line, specifically: perform stratigraphic interpretation on the optimal phase seismic data volume, and extrapolate the position of the pinch-out point identified after the stratigraphic interpretation, the extrapolated value is the error value between the simulated stratigraphic pinch-out point corresponding to the optimal phase wavelet and the set stratigraphic pinch-out point of the geological forward model, and use the new position after extrapolation as the identification result of the stratigraphic pinch-out point in the area to be identified.

[0050] according to Figure 2 , determine that the negative 90-degree phase is the optimal phase, and use the negative 90-degree phase data volume to identify the pinch-out point of the sand body, that is, convolve the negative 90-degree phase wavelet with the reflection coefficient to obtain the negative 90-degree phase seismic data volume, and perform stratigraphic interpretation on the data volume; in this embodiment, since the error value between the simulated stratigraphic pinch-out point corresponding to the negative 90-degree phase wavelet and the set stratigraphic pinch-out point of the geological forward model is 130 meters, the pinch-out point position (such as Figure 3 The area surrounded by the dashed line corresponds to the location of the pinch-out point) and is extrapolated by 130 meters, that is, Figure 3 The pinch-out point on the left side of the center is pushed 130 meters to the right, and the pinch-out point on the right side is pushed 130 meters to the left. The new position of the stratigraphic pinch-out point is determined as the identification result of the stratigraphic pinch-out point in the area to be identified (such as Figure 3 The area surrounded by the solid line in the middle corresponds to the pinch-out point position); the pinch-out point is identified for each profile in the three-dimensional seismic data of the area where the embodiment is located, and the identification results of multiple areas to be identified are superimposed to finally determine the plane position of the stratigraphic pinch-out line. Since the identification method of this embodiment not only interprets the horizon based on the seismic data body corresponding to the phase wavelet with the smallest error to minimize the existing error, but also compensates for the error by extrapolating the pinch-out point position identified after the horizon interpretation, thereby further eliminating the error and improving the accuracy of pinch-out point identification.

[0051] The plane position of the stratigraphic pinch-out line determined by superimposing the recognition results of multiple areas to be recognized is as follows: Figure 4 As shown in the figure, the white area is the area where the pinch-out occurs. Among them, the reflection coefficient is extracted from the post-stack pure wave seismic data volume. The 3D seismic must first be collected. The principle is to use the speed difference of sound waves propagating in different rock media. Generally, when selecting an area for 3D seismic, the blasting position and the geophone position must be designed first. After the design is completed, construction will be carried out and the explosive source will be excited. According to the records obtained by the set geophone, a series of complex processing processes such as denoising, stacking, offset, and homing are carried out by a special 3D seismic processing department. Finally, a set of superimposed pure wave data is provided to the interpreter, which is the post-stack pure wave data volume, which is the original seismic data volume that the interpreter can obtain.

[0052] The process of obtaining the reflection coefficient by this method is simple, and the result is consistent with the actual situation; the specific method of generating the plane position of the formation pinch-out line according to the position of the formation pinch-out point is an existing technology, so it will not be repeated here.

[0053] Embodiment of formation pinch-out point identification system

[0054] This embodiment provides a formation pinch-out point recognition system, which includes a processor, and the processor is used to execute program instructions to implement the formation pinch-out point recognition method in the above formation pinch-out point recognition method embodiment.

[0055] Since the specific working mode and principle of the formation pinch-out point identification system in this embodiment have been described in detail in the above-mentioned formation pinch-out point identification method embodiment, they will not be repeated here.

[0056] Embodiment of the method for identifying formation pinch-out line

[0057] This embodiment provides a formation pinch-out line identification method, which determines the plane position of the formation pinch-out line based on the result of formation pinch-out point identification. The formation pinch-out point identification adopts the formation pinch-out point identification method in the above-mentioned formation pinch-out point identification method embodiment.

[0058] Since the specific method and principle of the formation pinch-out line identification method in this embodiment have been described in detail in the above-mentioned formation pinch-out point identification method embodiment, they will not be repeated here.

[0059] Embodiment of formation pinch-out line identification system

[0060] This embodiment provides a formation pinch-out line identification system, which includes a processor, and the processor is used to execute program instructions to implement the formation pinch-out line identification method in the above-mentioned formation pinch-out line identification method embodiment.

[0061] Since the specific working mode and principle of the formation pinch-out line identification system in this embodiment have been described in detail in the above-mentioned formation pinch-out point identification method embodiment and formation pinch-out line identification method embodiment, they will not be repeated here.

[0062] The present invention has the following characteristics:

[0063] 1) By simulating the forward modeling of the geological conditions in the area to be identified, the error of different phase wavelets of the same frequency in identifying the pinch-out point is determined, and the horizon interpretation is performed based on the phase seismic data volume corresponding to the phase wavelet with the smallest error, so as to minimize the error in the identification of the pinch-out point;

[0064] 2) Not only is the stratigraphic interpretation performed based on the phase seismic data volume corresponding to the phase wavelet with the smallest error, but the position of the pinch-out point identified after the stratigraphic interpretation is extrapolated to compensate for the error in the identification of the pinch-out point based on the phase wavelet with the smallest error obtained by forward simulation, thereby further eliminating the error and improving the accuracy of pinch-out point identification.

[0065] It should be understood that the above specific embodiments of the present invention are only used to illustrate or explain the principles of the present invention, but do not constitute a limitation to the present invention.

Claims

1. A method for identifying formation pinch-out points, characterized in that: The steps include: 1) Perform a main frequency analysis on the seismic data of the area to be identified to determine the main frequency of earthquakes in the area to be identified; 2) Obtaining the well logging curves of the wells that have been drilled corresponding to the area to be identified, and determining the geophysical parameters of the area to be identified; 3) Determine the formation dip of the area to be identified; 4) setting the position of the stratigraphic pinch-out point of the geological forward model, and establishing the geological forward model by combining the geophysical parameters and stratigraphic dip of the area to be identified obtained in steps 2) and 3); 5) According to the earthquake main frequency of the area to be identified obtained in step 1) and the geological forward model established in step 4), forward simulation is performed using different phase wavelets within a set phase angle range to obtain the error between the simulated formation pinch-out point corresponding to each different phase wavelet and the set formation pinch-out point of the geological forward model; 6) Selecting the phase wavelet with the smallest error as the optimal phase wavelet, obtaining the optimal phase seismic data volume according to the optimal phase wavelet, and identifying the formation pinch-out point according to the optimal phase seismic data volume.

2. The method for identifying formation pinch-out points according to claim 1, characterized in that: The method for identifying the formation pinch-out point according to the optimal phase seismic data volume is: A stratigraphic interpretation is performed on the optimal phase seismic data volume, and the position of the pinch-out point identified after the stratigraphic interpretation is extrapolated, the extrapolated value being the error value between the simulated stratigraphic pinch-out point corresponding to the optimal phase wavelet and the set stratigraphic pinch-out point of the geological forward model, and the new position after extrapolation is used as the identification result of the stratigraphic pinch-out point in the area to be identified.

3. The method for identifying formation pinch-out points according to claim 1 or 2, characterized in that: Different phase wavelets within the set phase angle range are used for forward modeling, and the error between the simulated formation pinch-out point corresponding to each different phase wavelet and the set formation pinch-out point of the geological forward model is obtained as follows: Within the set phase angle range, at each interval of the set phase angle, the corresponding phase wavelet is taken once for forward simulation to obtain the simulated formation pinch-out points corresponding to each different phase wavelet, and the differences between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points are respectively used as the errors between the simulated formation pinch-out points corresponding to each different phase wavelet and the set formation pinch-out points of the geological forward model; the frequencies of the different phase wavelets are all the main frequencies of earthquakes in the area to be identified.

4. The method for identifying formation pinch-out points according to claim 1 or 2, characterized in that: The method of obtaining the optimal phase seismic data volume according to the optimal phase wavelet is: obtaining the optimal phase seismic data volume according to the convolution of the optimal phase wavelet and the reflection coefficient.

5. The method for identifying formation pinch-out points according to claim 4, characterized in that: The reflection coefficient is extracted from the post-stack pure-wave seismic data volume.

6. The method for identifying formation pinch-out points according to claim 1 or 2, characterized in that: The logging curves of the drilled wells mainly include sonic curves, density curves, natural gamma curves and neutron curves.

7. The method for identifying formation pinch-out points according to claim 1 or 2, characterized in that: There are at least two drilled wells, and the method for determining the formation dip of the area to be identified is: obtaining the burial depth and actual horizontal well spacing of the drilled wells corresponding to the area to be identified, and determining the geological layer passed by the drilled wells according to the burial depth of the drilled wells; and calculating the formation dip of the area to be identified according to the actual horizontal well spacing of each drilled well and the altitude difference of each drilled well corresponding to the same geological layer.

8. A formation pinch-out point identification system, characterized in that: It includes a processor, and the processor is used to execute program instructions to implement the formation pinch-out point identification method described in any one of claims 1 to 7.

9. A formation pinch-out line identification method, characterized in that: According to the result of formation pinch-out point identification, the plane position of the formation pinch-out line is determined, wherein the formation pinch-out point identification adopts the formation pinch-out point identification method as described in any one of claims 1 to 7 above.

10. A formation pinch-out line identification system, characterized in that: It includes a processor, and the processor is used to execute program instructions to implement the formation pinch-out line identification method described in claim 9.