Clastic rock sedimentary landform recovery method and device, medium and computing equipment

Abstract

The invention discloses a clastic rock sedimentary landform recovery method and device, a medium and computing equipment. The method comprises the following steps: establishing a clastic rock sedimentary landform interpretation research work area; acquiring a seismic stratigraphic interface of the target recovery interface; obtaining an interface seismic layering of an overlying stratum of the target stratum; obtaining a minimum isochronous seismic stratum interface of an overlying stratum of the target stratum; performing analysis based on the seismic event contact relationship, and obtaining the number of upper exceeding points of the target recovery interface; obtaining the micro-drift amount of the micro-drift surface; acquiring a spatial azimuth angle and a point inclination angle of discrete data points in a three-dimensional seismic target range; obtaining a point inclination angle and an azimuth angle of the target recovery interface and the micro floating surface; calculating the impression thickness of the target recovery interface; and recovering the sedimentary landform of the clastic rock. According to the method, the impression thickness of the target recovery interface can be calculated more accurately, then the clastic rock sedimentary landform is recovered, a scientific foundation is laid for well position deployment and the like, and the economic benefits of oil field exploration and development are improved.

Description

Technical Field

[0001] This invention relates to the field of oil and gas exploration technology, and in particular to a method, apparatus, medium and computing equipment for restoring clastic sedimentary landforms. Background Technology

[0002] The study of clastic sedimentary reservoirs is quite complex. Macroscopically speaking, tectonics controls sedimentation, and sediments reflect the state of paleotectonic structures. To clarify this relationship, it is essential to identify the intermediate connecting element – ​​paleogeomorphology. The spatial matching and configuration of paleogeomorphology and source area distribution controls the sedimentary sequence and filling pattern. For oil and gas researchers, geomorphological description, especially precise geomorphological characterization, can directly indicate the paleowater depth of a region. In the study of clastic oil and gas reservoirs, paleogeomorphological characterization, combined with vertical lithological combinations from drilling logging, allows for comprehensive analysis of the development locations of sandstone reservoirs, calcareous tight sandstones, and source rocks. Therefore, accurate paleogeomorphological mapping of clastic rocks is of great significance for locating oil and gas reservoirs and deploying well sites for oil and gas exploration and development.

[0003] Currently, there are numerous methods for paleogeographic reconstruction, including the imprinting method, the fill-in method, sedimentological analysis, and sequence stratigraphy. The imprinting method is the most commonly used. The basic principle of the imprinting method is the sedimentary compensation principle. The method involves first selecting the overlying strata as a reference surface, then calculating the thickness of the original strata between the reference surface and the target surface, and finally characterizing the paleogeographic morphology by using the calculated sedimentary thickness of the overlying strata.

[0004] However, the above method cannot find an overlying adjacent isochronous surface as a reference surface, which leads to the drawback of not being able to accurately track the isochronous reference surface and thus cannot fully restore the clastic sedimentary landform. Summary of the Invention

[0005] The main objective of this invention is to provide a method, apparatus, medium, and computing device for restoring clastic sedimentary landforms, so as to solve the technical problem that existing technologies cannot fully restore clastic sedimentary landforms.

[0006] To achieve the above objectives, the present invention provides a method for restoring clastic sedimentary landforms, comprising the following steps: S10, establishing a study area for interpreting clastic sedimentary landforms; S20, obtaining the seismic stratigraphic interface of the target restoration interface; S30, obtaining the seismic stratification of the interface of the overlying strata of the target strata; S40, obtaining the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata; S50, analyzing the contact relationship based on seismic phase axes to obtain the number of overshoot points of the target restoration interface; S60, obtaining the micro-float amount of the micro-float surface; S70, obtaining the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range; S80, obtaining the dip angles and azimuth angles of the target restoration interface and the micro-float surface; S90, calculating the impression thickness of the target restoration interface; S100, restoring the clastic sedimentary landform.

[0007] In some embodiments, step S10 includes the following steps: S110, based on seismic data of the study area, obtaining the seismic trace number and plane geodetic coordinates, obtaining the three-point coordinates of the seismic data, and obtaining the correspondence between the trace number and the geodetic coordinates, thereby establishing a clastic sedimentary geomorphological interpretation study area and loading seismic data; S120, based on the seismic profile, analyzing the continuity variation of the seismic reflection wave phase axis, selecting continuous data as research data, and loading it into the clastic sedimentary geomorphological interpretation study area; S130, obtaining wells that encounter target formation interfaces within the three-dimensional study area, selecting well data that are uniformly distributed within the three-dimensional area and have complete logging curves, and loading them into the clastic sedimentary geomorphological interpretation study area as another research data.

[0008] In some embodiments, step S20 includes the following steps: S210, using the acoustic and density curves of the well to obtain a synthetic seismic record; S220, based on the well-seismic calibration method, analyzing the seismic reflection characteristics of the target recovery interface; S230, based on the seismic reflection characteristics of the target recovery interface, performing seismic stratigraphic tracking and comparison to obtain the seismic stratigraphic interface of the target recovery interface.

[0009] In some embodiments, step S30 includes the following steps: S310, analyzing the correspondence between seismic strata and geological interfaces of the internal phase axis; S320, tracing the seismic strata interface of the internal phase axis close to the target recovery interface to obtain the interface seismic stratification of the overlying strata of the target strata.

[0010] In some embodiments, step S40 includes the following steps: S410, obtaining most wells with gamma value abrupt change points, connecting the wells together as the maximum reference surface; S420, starting from the maximum reference interface, sequentially searching downwards for more continuous phase axes in the three-dimensional region until the target recovery interface is the maximum interface, and taking the continuous phase axis closest to the target recovery interface as the minimum isochronous seismic stratigraphic interface of the overlying strata of the target formation.

[0011] In some embodiments, step S50 includes the following steps: S510, subtracting the minimum isochronous seismic stratigraphic interface of the overlying strata from the seismic stratigraphic interface of the target recovery interface to obtain the time-depth difference, and recording the time-depth difference with the largest value as the maximum thickness Thickmax; S520, starting from the maximum thickness Thickmax, analyzing the stratigraphic over-contact relationship between the seismic phase axis and the target recovery bottom interface, marking the over-contact points of the main seismic layers and the target recovery interface, and projecting them onto a horizontal plane; S530, starting from the maximum thickness Thickmax, obtaining the number of over-contact points O of the target recovery interface. i , where i is the data point number starting from the target recovery interface.

[0012] In some embodiments, step S60 includes the following steps: S610, based on the location point of the maximum thickness Thickmax, a profile is drawn, three points uniformly distributed in the seismic profile are selected, and the time span of a peak at the seismic target recovery interface is read to obtain the average value h. 平 S620, shift the target restoration boundary upwards by h 平 The calculation line for the sedimentary bottom surface restoration parameters is compared with the minimum isochronous interface line of the overlying strata to check whether the calculated parameter line is below the minimum isochronous interface line of the overlying strata. If it exceeds the minimum isochronous interface line of the overlying strata, the parameter interface is shifted downward by a certain amount of value. This process is iterated and traversed through all profiles to ensure that all adjusted calculated parameter lines are below the minimum isochronous interface line of the overlying strata. The surface formed by connecting all isochronous lines is denoted as the micro-drift surface. S630, the time values ​​of all points on the micro-drift surface are subtracted from the target restoration interface to obtain the drift amount of the micro-drift surface, denoted as ΔH. i .

[0013] In some embodiments, step S70 includes the following steps: S710, taking the minimum depth point of the minimum isochronous interface of the overlying strata as the top, the maximum time point of the target reconstructed paleogeographic strata as the bottom, and the three-dimensional work area line as the horizontal lateral boundary, a three-dimensional cube window is formed; S720, within the three-dimensional cube window, a seismic data point P0 is selected, and n points above and below the seismic data point P0 are selected respectively to form the amplitude of n+1 points, and the correlation of the n+1 points in the adjacent traces is obtained; S730, the correlation of the adjacent traces is iteratively calculated from top to bottom, and when the correlation is the maximum, the point of the adjacent trace is recorded as the correlation similarity starting point P1, and the depth and planar coordinates of the point are obtained; S740, using the spatial coordinate information of P0 and P1, a line is drawn, and the spatial azimuth and dip angle of each seismic data point are obtained when the angle between the line and the horizontal plane projection and the due north direction is calculated; S750, within the three-dimensional solid window, all points of the seismic data are obtained respectively, thereby obtaining the spatial azimuth and dip angle of the discrete data points of the three-dimensional seismic target range.

[0014] In some embodiments, step S80 includes the following steps: S810, using the coordinates and time points of the corresponding three-dimensional seismic event on the target recovery interface (L... i ,T j Z k ), obtain the tilt angle and azimuth angle of the target recovery interface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For the j-th channel number, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. S820 represents the azimuth angle of the spatial location ijk; using the micro-floating surface to correspond to the coordinates and time points of the three-dimensional seismic event (L... i ,T j Z k ), obtain the point tilt angle and azimuth angle of the micro-floating surface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For the j-th channel number, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. Let be the azimuth angle of the spatial position ijk.

[0015] In some embodiments, step S90 includes the following step: S910, obtaining a sedimentary stratum calculation reference surface based on the following formula:

[0016] In the formula: The thickness of the imprint at the target interface is defined by i, where i is the line direction number, j is the trace direction number, and k is the longitudinal depth sampling point number. Subtract the time values ​​of all points on the micro-floating surface from the target interface at the spatial ijk position. To recover the point tilt angle of the interface data for the target at spatial ijk location, To recover the azimuth angle of the interface data for the target at the spatial point ijk location, Let be the point tilt angle of the micro-drift surface data at location ijk in space. The azimuth angle of the micro-float data at spatial position ijk; S920, starting from the maximum thickness, accumulate along the line direction, using the super-point as the boundary control point, and calculate the cumulative thickness using the following formula. :

[0017] In the formula: i represents the data point from the target recovery interface. and These represent the number of times the target recovery interface at points i-1 and i has exceeded the specified threshold. and Let be the imprint thicknesses of the target recovery interface at points i-1 and i, respectively. and The cumulative thicknesses at points i-1 and i are respectively; S930, median filtering is used to eliminate... The abrupt jump point is used to obtain the smoothed data of the formation thickness, which is the imprint thickness of the target recovery interface.

[0018] In some embodiments, step S100 includes the following steps: S1100, obtaining a contour map of the impression thickness by using the inverse distance weighted interpolation method based on the discrete points of the impression thickness; S1200, dividing the sedimentary paleomorphology according to the thickness variation characteristics and the variation depth law of the clastic rock water body, and restoring the clastic rock sedimentary landform.

[0019] Furthermore, to achieve the above objectives, this application also provides a clastic sedimentary landform restoration device, which includes: a work area establishment module for establishing a clastic sedimentary landform interpretation and research work area; a first interface acquisition module for acquiring the seismic stratigraphic interface of the target restoration interface; a seismic stratification acquisition module for acquiring the interface seismic stratification of the overlying strata of the target strata; a second interface acquisition module for acquiring the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata; an overshoot point number acquisition module for analyzing the contact relationship of seismic phase axes and acquiring the overshoot point number of the target restoration interface; a micro-drift acquisition module for acquiring the micro-drift amount of the micro-drift surface; a first azimuth and point dip acquisition module for acquiring the spatial azimuth and point dip of discrete data points within the three-dimensional seismic target range; a second point dip and azimuth acquisition module for acquiring the point dip and azimuth of the target restoration interface and the micro-drift surface; an impression thickness calculation module for calculating the impression thickness of the target restoration interface; and a clastic sedimentary landform restoration module for restoring the clastic sedimentary landform.

[0020] In addition, to achieve the above objectives, embodiments of this application also provide a computer-readable storage medium including instructions that, when run on a computer, cause the computer to perform the clastic sedimentary landform restoration method of any embodiment of this application.

[0021] In addition, to achieve the above objectives, embodiments of this application also provide a computing device, which includes at least one processor, a memory, and an input / output unit; wherein the memory is used to store a computer program, and the processor is used to call the computer program stored in the memory to execute the clastic rock sedimentary landform restoration method of any embodiment of this application.

[0022] Compared with the prior art, the beneficial effects of the present invention are as follows: The clastic sedimentary landform restoration method provided in this application analyzes the contact relationship of seismic phase axes to obtain the number of overshoot points on the target restoration interface, thereby acquiring the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range. Using these spatial azimuth and dip angles, the dip and azimuth angles of the target restoration surface and the micro-float surface are calculated, allowing for a more accurate calculation of the impression thickness of the target restoration interface and thus restoring the clastic sedimentary landform. This method avoids the drawback of existing technologies that require finding an overlying adjacent isochronous surface as a reference surface, which leads to the inability to accurately track isochronous reference surfaces. This lays a scientific foundation for well location deployment and other aspects, thereby improving the economic efficiency of oilfield exploration and development. Attached Figure Description

[0023] Figure 1 This is a flowchart of a method for restoring clastic sedimentary landforms according to an embodiment of this application; Figure 2This is a structural block diagram of a clastic rock sedimentary landform restoration device provided in an embodiment of this application; Figure 3 This is a seismic profile provided according to an embodiment of the present invention; Figure 4 This is a point tilt angle cross-sectional view provided according to an embodiment of the present invention; Figure 5 This is an azimuth cross-sectional view provided according to an embodiment of the present invention; Figure 6 This is a schematic diagram of the mold thickness according to an embodiment of the present invention; Figure 7 This is a restored clastic sedimentary landform map provided according to an embodiment of the present invention; Figure 8 This is a seismic stratigraphic tracing profile provided according to an embodiment of the present invention; Figure 9 This is a schematic diagram of the structure of a medium provided in an embodiment of this application; Figure 10 This is a schematic diagram of the structure of a computing device provided in an embodiment of this application.

[0024] The realization of the purpose, functional features and advantages of this application will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0025] It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of the application. Rather, these embodiments are provided to make the disclosure more thorough and complete, and to fully convey the scope of the disclosure to those skilled in the art.

[0026] To address the aforementioned technical problems, embodiments of this application provide a method for restoring clastic sedimentary landforms. This method can be executed by a computer, such as... Figure 1 As shown, the method may include the following steps: S10, establish a study area for interpreting clastic sedimentary landforms.

[0027] In an exemplary embodiment, step S10 may include the following steps: S110, based on seismic data of the study area, obtains the seismic trace number and plane geodetic coordinates, obtains the three-point coordinates of the seismic trace number and the correspondence between the geodetic coordinates, thereby establishing the study area for the interpretation of clastic sedimentary landforms and loading seismic data; S120. Based on the seismic profile, the continuity of the seismic reflection wave phase axis is analyzed, and continuous data is selected as the research data and loaded into the clastic sedimentary geomorphological interpretation study area. S130: Obtain wells that encounter target formation interfaces within the three-dimensional work area, select well data that are uniformly distributed within the three-dimensional area and have complete logging curves, and load them into the clastic sedimentary geomorphological interpretation study area as another set of research data.

[0028] In this embodiment, the seismic data for the study area includes seismic stacking migration data and drilling data. A seismic profile, also called a seismic record profile, is a seismic data map indicating a specific seismic line, as shown below. Figure 3 As shown.

[0029] Furthermore, there are currently two basic formats for field seismic data: one is a multi-channel transmission record arranged in chronological order of sampling time, called a time-series record; the other is a record ordered by seismic trace, called a trace-series record. Generally, the time-series field data is converted into a trace-series record through decompilation, and then the seismic records are stored in the order of shot and trace. Each seismic trace consists of two parts: the trace head and the data, with the trace number stored therein to distinguish the location of the receiver.

[0030] S20, Obtain the seismic stratigraphic interface of the target recovery interface.

[0031] In an exemplary embodiment, step S20 may include the following steps: S210 uses the acoustic and density curves of the well to obtain synthetic seismic records; S220, based on well-seismic calibration method, analyzes the seismic reflection characteristics of the target recovery interface; S230, based on the seismic reflection characteristics of the seismic reflection interface of the target recovery interface, the seismic strata are tracked and compared to obtain the seismic stratigraphic interface of the target recovery interface.

[0032] Specifically, well-seismic calibration serves as a bridge connecting well logging and seismic information. The accuracy of the calibration results directly determines the accuracy of seismic inversion. Therefore, calibration should aim for the highest possible correlation between well logging and seismic information, especially ensuring precise calibration of the target layer, and the calibration results should be self-checked and verified. Typically, the influence of wavelet waves must be directly considered during calibration.

[0033] S30, obtain the seismic layering of the interface of the overlying strata of the target stratum.

[0034] In an exemplary embodiment, step S30 may specifically include the following steps: S310, Analyze the correspondence between seismic strata and geological interfaces of the internal phase axis; S320, track the seismic stratigraphic interface of the internal phase axis near the target recovery interface to obtain the interface seismic stratification of the overlying strata of the target strata.

[0035] Specifically, when analyzing the correspondence between seismic strata and geological interfaces of internal phase axes, the seismic strata interface of the target recovery interface should be taken as the bottom boundary, and the minimum isochronous traceable interface of the overlying strata should be taken as the top boundary.

[0036] S40, obtain the minimum isochronous seismic stratigraphic interface of the overlying strata of the target stratum.

[0037] In an exemplary embodiment, step S40 may specifically include the following steps: S410, Obtain the gamma value abrupt change points in most wells, connect the wells together as the maximum reference surface; S420, starting from the maximum reference interface, search downwards for relatively continuous phase axes in the three-dimensional region until the target recovery interface is the maximum interface, and take the continuous phase axis closest to the target recovery interface as the minimum isochronous seismic stratigraphic interface of the overlying strata of the target stratum.

[0038] In step S410, the drilling spontaneous potential curve or spontaneous gamma curve is used, with the target recovery interface as the bottom surface, and the sedimentation cycle is divided from bottom to top to obtain the gamma value mutation point in most wells.

[0039] Furthermore, the gamma value is a physical property of the display; it is fixed, unchanging, and uncorrectable. In this exemplary embodiment, the gamma value refers to the specific numerical value of the distortion of the input signal by the output image of the display.

[0040] S50 analyzes the contact relationship of seismic phase axes to obtain the number of overshoot points on the target recovery interface.

[0041] In an exemplary embodiment, step S50 may specifically include the following steps: S510, subtract the minimum isochronous seismic stratigraphic interface of the overlying strata from the seismic stratigraphic interface of the target recovery interface to obtain the time-depth difference, and record the time-depth difference with the largest value as the maximum thickness Thickmax. S520, starting from the maximum thickness Thickmax, analyzes the stratigraphic over-contact relationship between the seismic phase axis and the target reconstructed bottom interface, marking the over-contact points of the main seismic layers and the target reconstructed interface, projecting them onto the horizontal plane, which can be used as a reference. Figure 8 In the diagram, Upward 1, Upward 2, and Upward 3 represent the three upward points indicated by the arrows, H. 微 H represents a slightly floating surface. 目 H represents the target layer. 最大 H represents the maximum reference plane. 最小 The smallest isochronous surface is indicated; the H appearing in the other figures... 微 H 目 H 最大 H最小 All of the above content; S530, starting from the maximum thickness Thickmax, obtain the number of times the target recovery interface exceeds the upper limit point O. i , where i is the data point number starting from the target recovery interface.

[0042] Specifically, in step S520, a profile can be taken along the seismic line direction, and the stratigraphic over-contact relationship between the seismic phase axis and the target recovery bottom interface can be analyzed between the seismic stratigraphic interface and the maximum reference surface of the target recovery interface, thereby marking the over-contact points between the main seismic layers and the target recovery interface.

[0043] S60, obtain the amount of micro-float.

[0044] In an exemplary embodiment, step S60 may specifically include the following steps: S610, based on the location of the maximum thickness Thickmax, pulls out the profile, selects three points uniformly distributed in the seismic profile, and reads the time span of a peak at the seismic target recovery interface to obtain the average value h. 平 ; S620, shift the target restoration boundary upwards by h 平 As the calculation line for the restoration parameters of the sedimentary bottom surface, compared with the minimum isochronous interface line of the overlying strata, all the adjusted parameter calculation lines are below the minimum isochronous line of the overlying strata, and the surface formed by connecting all the isochronous lines is called the micro-floating surface. S630, subtract all the point time values ​​on the micro-drift surface from the target recovery interface to obtain the drift amount of the micro-drift surface, denoted as △H. i .

[0045] A wave crest is the maximum amplitude of a wave within a certain wavelength range, while the opposite minimum amplitude is called a wave trough. Taking a transverse wave as an example, the highest point of the bulge is the wave crest, and the lowest point of the trough is the wave trough.

[0046] Further, in step S620, after comparing the calculated parameters of the sedimentary bottom surface with the minimum isochronous interface line of the overlying strata, the parameter calculation line below the minimum isochronous line of the overlying strata is obtained.

[0047] In addition, if there are parameter calculation lines that exceed the minimum isochronous line of the overlying strata, the parameter calculation lines are shifted downward by a certain amount of value, and all profiles are iterated and examined in turn to ensure that all adjusted parameter calculation lines are below the minimum isochronous line of the overlying strata.

[0048] S70 acquires the spatial azimuth and dip angles of discrete data points within the 3D seismic target area.

[0049] In an exemplary embodiment, step S70 may specifically include the following steps: S710, with the minimum depth point of the minimum isochronous interface of the overlying strata as the top, the maximum time point of the target restored paleogeographic strata as the bottom, and the three-dimensional work area line as the horizontal lateral boundary, forms a three-dimensional cube window. S720, within the three-dimensional cube window, select seismic data point P0, and select n points above and below the seismic data point P0 respectively to form the amplitude of n+1 points, and obtain the correlation of the n+1 points in adjacent traces; S730, the correlation between adjacent channels is iteratively calculated from top to bottom. When the correlation is the maximum, the point of the adjacent channel is recorded as the correlation similarity starting point P1, and the depth and planar coordinates of the point are obtained. S740, using the spatial coordinate information of P0 and P1 to draw a line, when the angle between the line and the horizontal plane projection and the due north direction is obtained, the spatial azimuth and dip angle of each seismic data point are obtained; S750: Within the three-dimensional window, all points of the seismic data are acquired, thereby obtaining the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range.

[0050] Amplitude refers to the maximum value that a vibrating physical quantity can reach, usually denoted by A, representing the range and intensity of the vibration. Amplitude is a scalar quantity, having only magnitude and no direction, and is expressed in units of length such as meters or centimeters. The larger the amplitude, the greater the deviation of the object from its equilibrium position during vibration, and the greater the energy of the vibration.

[0051] Furthermore, in a system, the movement of one point causes changes in the movement of other points; this is known as point correlation.

[0052] S80, obtain the tilt angle and azimuth angle of the target recovery interface and the micro-floating surface.

[0053] In an exemplary embodiment, reference is made to Figure 4 and Figure 5 Step S80 may specifically include the following steps: S810, using the coordinates and time points of the corresponding three-dimensional earthquake on the target recovery interface (L i ,T j Z k ), obtain the tilt angle and azimuth angle of the target recovery interface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For the j-th channel number, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. Let be the azimuth angle of the spatial position ijk; S820 utilizes micro-floating surfaces to correspond to the coordinates and time points of a 3D seismic event (L...). i ,T j Z k ), obtain the point tilt angle and azimuth angle of the micro-floating surface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For the j-th channel number, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. Let be the azimuth angle of the spatial position ijk.

[0054] Among them, the 3D seismic survey is conducted in a linear pattern, with surface element reception, and is a seismic survey conducted in both longitudinal and transverse directions. The processed results can be cut into cross-sections in any direction.

[0055] S90, calculate the imprint thickness of the target recovery interface.

[0056] In an exemplary embodiment, step S90 may specifically include the following steps: S910, the reference surface for sedimentary strata calculation is obtained based on the following formula: (1) In the formula: The thickness of the imprint at the target interface is defined by i, where i is the line direction number, j is the trace direction number, and k is the longitudinal depth sampling point number. Subtract the time values ​​of all points on the micro-floating surface from the target interface at the spatial ijk position. To recover the point tilt angle of the interface data for the target at spatial ijk location, To recover the azimuth angle of the interface data for the target at the spatial point ijk location, Let be the point tilt angle of the micro-drift surface data at location ijk in space. The azimuth angle of the spatial ijk position micro-drift surface data; S920, starting from the maximum thickness, accumulate along the track direction, using the super-point as the boundary control point, and calculate the cumulative thickness using the following formula. : (2) In the formula: i represents the data point from the target recovery interface. and These represent the number of times the target recovery interface at points i-1 and i has exceeded the specified threshold. and Let be the imprint thicknesses of the target recovery interface at points i-1 and i, respectively. and These are the cumulative thicknesses at points i-1 and i, respectively; S930, using median filtering to eliminate The abrupt jump point is identified, and the smoothed formation thickness data is obtained. This smoothed formation thickness data represents the imprint thickness of the target recovery interface. A schematic diagram of the imprint thickness is shown below. Figure 6 As stated above.

[0057] Among them, sedimentary strata refer to strata on the periphery of quartz sand grains where inert films are deposited. The sedimentary layer contains fine-grained materials, dead clay, and ash, which are firmly aggregated on the sand grains to form an uneven outer layer, reducing the flowability of the sand grains, decreasing the energy transfer efficiency of the molding sand, and causing uneven compaction of the sand mold.

[0058] Furthermore, median filtering is a nonlinear signal processing technique based on sorting statistics theory that effectively suppresses noise. The basic principle of median filtering is to replace the value of a point in a digital image or sequence with the median value of all points in its neighborhood, making the surrounding pixel values ​​closer to the true value, thereby eliminating isolated noise points. The method involves using a two-dimensional sliding template of a certain structure to sort the pixels within the template according to their pixel values, generating a monotonically increasing (or decreasing) sequence as a two-dimensional data sequence.

[0059] S100, restoring the clastic sedimentary landform.

[0060] In an exemplary embodiment, step S100 may specifically include the following steps: S1100: Discretize the imprint thickness points and obtain the imprint thickness contour map based on the inverse distance weighted interpolation method; Based on the thickness variation characteristics and the depth variation patterns of clastic rock water bodies, sedimentary paleomorphology was divided in S1200, and the clastic sedimentary landform was reconstructed. The reconstructed clastic sedimentary landform map is shown below. Figure 7 As shown.

[0061] Specifically, inverse distance weighted interpolation, or IDW (Inverse Distance Weight), also known as the inverse distance multiplication method, is a weighted average interpolation method that can perform exact or smooth interpolation. The power parameter controls how the weight coefficients decrease as the distance from a grid node increases. For a larger power, closer data points are given a higher weight share, while for a smaller power, the weights are distributed more evenly among the data points.

[0062] The clastic sedimentary landform restoration method provided in this application analyzes the contact relationship of seismic phase axes to obtain the number of overshoot points on the target restoration interface, thereby acquiring the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range. Using these spatial azimuth and dip angles, the dip and azimuth angles of the target restoration surface and the micro-float surface are calculated, allowing for more accurate calculation of the impression thickness of the target restoration interface and thus restoring the clastic sedimentary landform. This method avoids the drawback of existing technologies that require finding an overlying adjacent isochronous surface as a reference surface, which leads to the inability to accurately track isochronous reference surfaces. This lays a scientific foundation for well location deployment and other aspects, thereby improving the economic benefits of oilfield exploration and development.

[0063] Based on the above embodiments, refer to Figure 2 Another embodiment of this application also provides a clastic rock sedimentary landform restoration device, which may include the following modules: Module 210 is used to establish a comprehensive work area for interpreting clastic sedimentary landforms. The first interface acquisition module 220 is used to acquire the seismic stratigraphic interface of the target recovery interface; Seismic layering acquisition module 230 is used to acquire the interface seismic layering of the overlying strata of the target strata. The second interface acquisition module 240 is used to acquire the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata. The module 250 for obtaining the number of times the upper cross points are used to analyze the contact relationship of the seismic phase axis and obtain the number of times the upper cross points are obtained on the target recovery interface. The micro-drift quantity acquisition module 260 is used to acquire the micro-drift quantity of the micro-drift surface. The first azimuth and dip angle acquisition module 270 is used to acquire the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range; The second point tilt angle and azimuth angle acquisition module 280 is used to acquire the point tilt angle and azimuth angle of the target recovery interface and the micro-floating surface; The impression thickness calculation module 290 is used to calculate the impression thickness of the target recovery interface; Clastic sedimentary landform restoration module 2100 is used to restore clastic sedimentary landforms.

[0064] Based on the above embodiments, this application also provides a computer-readable storage medium, see reference. Figure 9The computer-readable storage medium shown is an optical disc 50, on which a computer program (i.e., a program product) is stored. When the computer program is run by a processor, it will implement the steps described in the above-described method implementation, such as: establishing a study area for interpreting clastic sedimentary landforms; obtaining the seismic stratigraphic interface of the target restoration interface; obtaining the interface seismic stratification of the overlying strata of the target strata; obtaining the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata; analyzing the contact relationship based on seismic phase axes to obtain the number of overshoot points of the target restoration interface; obtaining the micro-float amount; obtaining the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range; obtaining the dip angles and azimuth angles of the target restoration interface and the micro-float; calculating the impression thickness of the target restoration interface; and restoring the clastic sedimentary landforms. The specific implementation methods of each step will not be repeated here.

[0065] It should be noted that examples of computer-readable storage media may also include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other optical and magnetic storage media, which will not be elaborated here.

[0066] Furthermore, based on the above embodiments, this application also provides a computing device. Figure 10 A block diagram is shown of an exemplary computing device 60 suitable for implementing embodiments of the present application. The computing device 60 may be a computer system or a server. Figure 10 The computing device 60 shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this application.

[0067] like Figure 10 As shown, the components of computing device 60 may include, but are not limited to: one or more processors or processing units 601, system memory 602, and bus 603 connecting different system components (including system memory 602 and processing unit 601).

[0068] The computing device 60 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the computing device 60, including volatile and non-volatile media, removable and non-removable media.

[0069] System memory 602 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 6021 and / or cache memory 6022. Computing device 60 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, ROM 6023 may be used to read and write non-removable, non-volatile magnetic media (…). Figure 10 (Not shown in the image, usually referred to as "hard drive"). Although not shown in... Figure 10 The diagram illustrates that a disk drive for reading and writing to removable non-volatile disks (e.g., "floppy disks") and an optical disk drive for reading and writing to removable non-volatile optical disks (e.g., CD-ROMs, DVD-ROMs, or other optical media) can be provided. In these cases, each drive can be connected to a bus 603 connecting different system components via one or more data media interfaces. The system memory 602 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of the embodiments of this application.

[0070] A program / utility 6025 having a set (at least one) of program modules 6024 may be stored, for example, in system memory 602, and such program modules 6024 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment. Program modules 6024 typically perform the functions and / or methods described in the embodiments of this application.

[0071] The computing device 60 can also communicate with one or more external devices 604 (such as a keyboard, pointing device, display, etc.). This communication can be performed via input / output (I / O) interface 605. Furthermore, the computing device 60 can also communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via network adapter 606. Figure 10 As shown, network adapter 606 communicates with other modules of computing device 60 (such as processing unit 601, etc.) via bus 603, which connects different system components. It should be understood that, although... Figure 10 As not shown, it may be used in conjunction with computing device 60 with other hardware and / or software modules.

[0072] The processing unit 601 executes various functional applications and data processing by running programs stored in the system memory 602. For example, it establishes a comprehensive study area for interpreting clastic sedimentary landforms; obtains the seismic stratigraphic interface of the target restoration interface; establishes the study area for interpreting clastic sedimentary landforms; obtains the seismic stratigraphic interface of the target restoration interface; obtains the seismic stratification of the interface of the overlying strata of the target strata; obtains the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata; analyzes the contact relationship based on seismic phase axes to obtain the number of overshoot points of the target restoration interface; obtains the micro-float amount; obtains the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range; obtains the dip angles and azimuth angles of the target restoration interface and the micro-float; calculates the impression thickness of the target restoration interface; and restores the clastic sedimentary landform. The specific implementation methods of each step will not be repeated here. It should be noted that although several units / modules or sub-units / sub-modules of the clastic sedimentary landform restoration device are mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, according to the embodiments of this application, the features and functions of two or more units / modules described above can be embodied in one unit / module. Conversely, the features and functions of one unit / module described above can be further divided into multiple units / modules for embodiment.

[0073] In the description of this application, it should be noted that the terms "first", "second", and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0074] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0075] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative; for example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interface; the indirect coupling or communication connection between apparatuses or units may be electrical, mechanical, or other forms.

[0076] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0077] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0078] If a function is implemented as a software functional unit and sold or used as an independent product, it can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0079] Finally, it should be noted that the above embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The protection scope of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the technical scope disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the protection scope of this application. Therefore, the protection scope of this application should be determined by the protection scope of the claims.

[0080] Furthermore, although the operations of the method of this application are described in a specific order in the accompanying drawings, this does not require or imply that these operations must be performed in that specific order, or that all the operations shown must be performed to achieve the desired result. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

Claims

1. A method for restoring clastic sedimentary landforms, characterized in that, The method includes the following steps: S10, establish a research area for interpreting clastic sedimentary landforms; S20, Obtain the seismic stratigraphic interface of the target recovery interface; S30, obtain the seismic layering of the interface of the overlying strata of the target strata; S40, Obtain the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata; S50, based on the contact relationship of seismic phase axes, analyzes and obtains the number of overshoot points on the target recovery interface; S60, obtain the micro-float amount; S70, acquires the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target area; S80, obtain the tilt angle and azimuth angle of the target recovery interface and the micro-floating surface; S90, calculate the imprint thickness of the target recovery interface; S100, restoring the clastic sedimentary landform.

2. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S10 includes the following steps: S110, based on seismic data of the study area, obtains the seismic trace number and plane geodetic coordinates, obtains the three-point coordinates of the seismic trace number and the correspondence between the geodetic coordinates, thereby establishing the study area for the interpretation of clastic sedimentary landforms and loading seismic data; S120. Based on the seismic profile, the continuity of the seismic reflection wave phase axis is analyzed, and continuous data is selected as the research data and loaded into the clastic sedimentary geomorphological interpretation study area. S130: Obtain wells that encounter target formation interfaces within the three-dimensional work area, select well data that are uniformly distributed within the three-dimensional area and have complete logging curves, and load them into the clastic sedimentary geomorphological interpretation study area as another set of research data.

3. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S20 includes the following steps: S210 uses the acoustic and density curves of the well to obtain synthetic seismic records; S220, based on well-seismic calibration method, analyzes the seismic reflection characteristics of the target recovery interface; S230, based on the seismic reflection characteristics of the seismic reflection interface of the target recovery interface, the seismic strata are tracked and compared to obtain the seismic stratigraphic interface of the target recovery interface.

4. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S30 includes the following steps: S310, Analyze the correspondence between seismic strata and geological interfaces of the internal phase axis; S320, track the seismic stratigraphic interface of the internal phase axis near the target recovery interface to obtain the interface seismic stratification of the overlying strata of the target strata.

5. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S40 includes the following steps: S410, obtain the gamma value abrupt change points in most wells, connect the wells together as the maximum reference surface; S420, starting from the maximum reference interface, search downwards for relatively continuous phase axes in the three-dimensional region until the target recovery interface is the maximum interface, and take the continuous phase axis closest to the target recovery interface as the minimum isochronous seismic stratigraphic interface of the overlying strata of the target stratum.

6. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S50 includes the following steps: S510, subtract the minimum isochronous seismic stratigraphic interface of the overlying strata from the seismic stratigraphic interface of the target recovery interface to obtain the time-depth difference, and record the time-depth difference with the largest value as the maximum thickness Thickmax. S520, starting from the maximum thickness Thickmax, analyze the stratigraphic over-contact relationship between the seismic phase axis and the target recovery interface, mark the over-contact points of the main seismic layers and the target recovery interface, and project them onto the horizontal plane. S530, starting from the maximum thickness Thickmax, obtain the number of times the target recovery interface exceeds the upper limit point O. i , where i is the data point number starting from the target recovery interface.

7. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S60 includes the following steps: S610, based on the location of the maximum thickness Thickmax, pulls out the profile, selects three points uniformly distributed in the seismic profile, and reads the time span of a peak at the seismic target recovery interface to obtain the average value h. 平 ; S620, shift the target restoration boundary upwards by h 平 As the calculation line for the restoration parameters of the sedimentary bottom surface, it is compared with the minimum isochronous interface line of the overlying strata to check whether the parameter calculation line is below the minimum isochronous line of the overlying strata. If it exceeds the minimum isochronous line of the overlying strata, the parameter interface is shifted downward by a certain amount of value. All profiles are iterated and checked in turn to ensure that all parameter calculation lines after adjustment are below the minimum isochronous line of the overlying strata. The surface formed by connecting all isochronous lines is called the micro-floating surface. S630, subtract all the point time values ​​on the micro-drift surface from the target recovery interface to obtain the drift amount of the micro-drift surface, denoted as △H. i .

8. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S70 includes the following steps: S710, with the minimum depth point of the minimum isochronous interface of the overlying strata as the top, the maximum time point of the target to restore the ancient landform strata as the bottom, and the three-dimensional work area line as the horizontal boundary, forms a three-dimensional cube window. S720, within the three-dimensional cube window, select seismic data point P0, and select n points above and below the seismic data point P0 respectively to form the amplitude of n+1 points, and obtain the correlation of the n+1 points in adjacent traces; S730, the correlation between adjacent channels is iteratively calculated from top to bottom. When the correlation is the maximum, the point of the adjacent channel is recorded as the correlation similarity starting point P1, and the depth and planar coordinates of the point are obtained. S740, using the spatial coordinate information of P0 and P1 to draw a line, when the angle between the line and the horizontal plane projection and the due north direction is obtained, the spatial azimuth and dip angle of each seismic data point are obtained; S750: Within the three-dimensional window, all points of the seismic data are acquired, thereby obtaining the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range.

9. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S80 includes the following steps: S810, using the coordinates and time points of the corresponding three-dimensional earthquake on the target recovery interface (L i ,T j Z k ), obtain the tilt angle and azimuth angle of the target recovery interface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For channel j, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. Let be the azimuth angle of the spatial position ijk; S820 utilizes micro-floating surfaces to correspond to the coordinates and time points of a 3D seismic event (L...). i ,T j Z k ), obtain the point tilt angle and azimuth angle of the micro-floating surface, denoted as (L i ,T j Z k , , ), where i is the line direction index, j is the trace direction index, k is the longitudinal depth sampling point number, and L i For line number i, T j For channel j, Z k This is a depth value. Let be the tilt angle of the point at position ijk in space. Let be the azimuth angle of the spatial position ijk.

10. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S90 includes the following steps: S910, the reference surface for sedimentary strata calculation is obtained based on the following formula: In the formula: The thickness of the imprint at the target interface is defined by i, where i is the line direction number, j is the trace direction number, and k is the longitudinal depth sampling point number. Subtract the time values ​​of all points on the micro-floating surface from the target interface at the spatial ijk position. To recover the point tilt angle of the interface data for the target at spatial ijk location, To recover the azimuth angle of the interface data for the target at the spatial point ijk location, Let be the point tilt angle of the micro-drift surface data at location ijk in space. The azimuth angle of the spatial ijk position micro-drift surface data; S920, starting from the maximum thickness, accumulate along the track direction, using the super-point as the boundary control point, and calculate the cumulative thickness using the following formula. : In the formula: i represents the data point from the target recovery interface. and These represent the number of times the target recovery interface at points i-1 and i has exceeded the specified threshold. and Let be the imprint thicknesses of the target recovery interface at points i-1 and i, respectively. and These are the cumulative thicknesses at points i-1 and i, respectively; S930, using median filtering to eliminate The abrupt jump point is used to obtain the smoothed data of the formation thickness, which is the imprint thickness of the target recovery interface.

11. The method for restoring clastic sedimentary landforms according to claim 1, characterized in that, Step S100 includes the following steps: S1100: Discretize the imprint thickness points and obtain the imprint thickness contour map based on the inverse distance weighted interpolation method; S1200, based on the characteristics of thickness variation and the depth variation of clastic rock water bodies, is divided into sedimentary paleomorphology and the clastic rock sedimentary landform is restored.

12. A device for restoring clastic sedimentary landforms, characterized in that, include: The work area establishment module is used to establish work areas for interpreting clastic sedimentary landforms. The first interface acquisition module is used to acquire the seismic stratigraphic interface of the target recovery interface; The seismic layering acquisition module is used to acquire the interface seismic layers of the overlying strata of the target strata; The second interface acquisition module is used to acquire the minimum isochronous seismic stratigraphic interface of the overlying strata of the target strata. The module for obtaining the number of times the upper boundary point is used to analyze the contact relationship of the seismic phase axis and obtain the number of times the upper boundary point is obtained on the target recovery interface. The micro-drift quantity acquisition module is used to acquire the micro-drift quantity of the micro-drift surface; The first azimuth and dip angle acquisition module is used to acquire the spatial azimuth and dip angles of discrete data points within the three-dimensional seismic target range; The second point tilt angle and azimuth angle acquisition module is used to acquire the point tilt angle and azimuth angle of the target recovery interface and the micro-floating surface; The impression thickness calculation module is used to calculate the impression thickness of the target recovery interface; The clastic sedimentary landform restoration module is used to restore clastic sedimentary landforms.

13. A computer-readable storage medium, characterized in that, It includes instructions that, when run on a computer, cause the computer to perform the clastic sedimentary landform restoration method according to any one of claims 1 to 11.

14. A computing device, characterized in that, The computing device includes: At least one processor, memory, and input / output unit; The memory is used to store computer programs, and the processor is used to call the computer programs stored in the memory to execute the clastic rock sedimentary landform restoration method according to any one of claims 1 to 11.

Images