CN-121995480-A - Clastic rock deposition landform restoration method, device, medium and computing equipment

Abstract

The invention discloses a clastic rock deposition landform restoration method, a device, a medium and computing equipment. The method comprises the steps of establishing a clastic rock sedimentary landform interpretation research work area, obtaining an earthquake stratum interface of a target recovery interface, obtaining an interface earthquake layering of an overlying stratum of the target stratum, obtaining a minimum isochronous earthquake stratum interface of the overlying stratum of the target stratum, analyzing based on an earthquake on-phase axis contact relation, obtaining the number of times of superpoints of the target recovery interface, obtaining micro-floating quantity of a micro-floating surface, obtaining an attitude angle and a point inclination angle of discrete data points of a three-dimensional earthquake target range, obtaining the point inclination angle and the azimuth angle of the target recovery interface and the micro-floating surface, calculating the impression thickness of the target recovery interface, and recovering the clastic rock sedimentary landform. The method can calculate the thickness of the impression of the target recovery interface more accurately, so as to recover the sedimentary landform of the clastic rock, lay a scientific foundation for well position deployment and the like, and improve the economic benefit of oil field exploration and development.

Inventors

• FU LEI
• WEI PU
• Ji Baoqiang
• CHEN JING
• LI JUNFEI
• XU GANG
• ZHANG CHUNGUANG
• WANG XUEHUA
• ZHONG YANG

Assignees

• 中国石油天然气股份有限公司

Dates

Publication Date

20260508

Application Date

20241107

Claims (14)

1. 1. A clastic rock deposition topography restoration method, which is characterized by comprising the following steps: s10, establishing a clastic rock deposition landform interpretation research work area; s20, acquiring a seismic stratum interface of a target recovery interface; S30, obtaining the interfacial seismic stratification of the stratum overlying the target stratum; s40, acquiring a minimum isochronous earthquake stratum interface of the stratum covered by the target stratum; S50, analyzing based on the earthquake phase axis contact relation to obtain the number of times of the upper super point of the target recovery interface; s60, obtaining micro-floating quantity of the micro-floating surface; s70, acquiring the attitude angle and the dot inclination angle of discrete data points of the three-dimensional seismic target range; S80, acquiring the point inclination angle and the azimuth angle of the target recovery interface and the micro-floating surface; s90, calculating the thickness of the impression of the target recovery interface; S100, recovering the sedimentary landform of the clastic rock.
2. 2. The clastic rock deposition profile restoration method according to claim 1, wherein step S10 includes the steps of: S110, based on the seismic data of the research area, acquiring the channel number and the plane geodetic coordinates of the earthquake, acquiring the three-point coordinates of the earthquake, and acquiring the corresponding relation between the channel number and the geodetic coordinates, thereby establishing a clastic rock sedimentary landform interpretation research work area, and loading the seismic data; s120, analyzing the continuous change condition of the same phase axis of the seismic reflection wave based on the seismic section, selecting the continuous data as research data, and loading the research data into the clastic rock sedimentary landform interpretation research work area; S130, obtaining wells which are drilled in the three-dimensional work area and meet the target to recover the stratum interface, selecting well data which are uniformly distributed in the three-dimensional area and have complete logging curves, and loading the well data as another research data to the clastic rock deposit landform interpretation research work area.
3. 3. The clastic rock deposition profile restoration method according to claim 1, wherein step S20 includes the steps of: S210, acquiring an artificial synthetic seismic record by using a sound wave and density curve of the well drilling; s220, analyzing the seismic reflection characteristics of the seismic reflection interface of the target recovery interface based on a well earthquake calibration method; s230, tracking and comparing the seismic stratum based on the seismic reflection characteristics of the seismic reflection interface of the target recovery interface, and obtaining the seismic stratum interface of the target recovery interface.
4. 4. The clastic rock deposition profile restoration method according to claim 1, wherein step S30 includes the steps of: S310, analyzing the corresponding relation between the earthquake stratum of the internal phase shaft and the geological interface; S320, tracking the seismic stratum interface of which the internal phase axis is close to the target recovery interface so as to acquire the interfacial seismic stratification of the stratum covered by the target stratum.
5. 5. The clastic rock deposition profile restoration method according to claim 1, wherein step S40 includes the steps of: S410, acquiring gamma value mutation points of most wells, and connecting the wells to be used as the maximum reference surface; S420, sequentially searching continuous homophase axes of the three-dimensional region downwards by taking the maximum reference interface as an initial interface until the target recovery interface is the maximum interface, and taking the continuous homophase axes closest to the target recovery interface as the minimum isochronous earthquake stratum interface of the stratum overlying the target stratum.
6. 6. The clastic rock deposition profile restoration method according to claim 1, wherein step S50 includes the steps of: S510, subtracting the minimum equivalent seismic stratum interface of the overlying stratum from the seismic stratum interface of the target recovery interface, obtaining a time depth difference, and recording the time depth difference with the maximum value as a maximum thickness Thickmax; S520, starting from the position of the maximum thickness Thickmax, analyzing the upper super-contact relation between the stratum of the same phase axis of the earthquake and the target recovery bottom interface, marking the upper super-point of the main earthquake layering and the target recovery interface, and projecting the upper super-point on a horizontal plane; S530, respectively obtaining the number O i of the superpoints of the target recovery interface by taking the position of the maximum thickness Thickmax as the beginning, wherein i is the number of data points from the target recovery interface.
7. 7. The clastic rock deposition profile restoration method according to claim 1, wherein step S60 includes the steps of: S610, pulling a section based on the position point of the maximum thickness Thickmax, selecting three points uniformly distributed in the seismic section, and reading the time span of a wave crest at a seismic target recovery interface to obtain an average value h Flat plate ; S620, upwards translating a target recovery limit by h Flat plate to serve as a deposited bottom surface recovery parameter calculation line, comparing the deposited bottom surface recovery parameter calculation line with an overlying stratum minimum isochrone, checking whether the parameter calculation line is below the overlying stratum minimum isochrone, if the parameter calculation line exceeds the overlying stratum minimum isochrone, downwards translating the parameter interface by a certain amount of numerical values, sequentially iterating and checking all sections, ensuring that all the adjusted parameter calculation lines are below the overlying stratum minimum isochrone, and recording the surface formed by connecting all the isochrones as a micro-floating surface; S630, subtracting all point time values on the micro-drift surface from a target recovery interface to obtain the micro-drift surface drift amount, and recording the drift amount as delta H i .
8. 8. The clastic rock deposition profile restoration method according to claim 1, wherein step S70 includes the steps of: S710, taking the minimum depth point of the interface when the overlying stratum is minimum as the top, taking the maximum time point of the target recovery paleomorphic stratum as the bottom, and taking the three-dimensional work area line as the horizontal and transverse boundary to form a three-dimensional cube window; S720, in a three-dimensional cube window, selecting a seismic data point P 0 , respectively selecting upper and lower n points at the seismic data point P 0 to form the amplitude of n+1 points, and obtaining the correlation of the n+1 points in adjacent channels; S730, iteratively solving the correlation from top to bottom of the adjacent channels, and when the correlation is maximum, marking the point of the adjacent channels as a correlation similarity starting point P 1 and obtaining the depth and plane coordinates of the point; S740, using the space coordinate information of the P 0 and the P 1 as connecting lines, and obtaining the space azimuth angle and the point inclination angle of each seismic data point when the connecting lines form an included angle with the horizontal plane projection and the north direction; S750, respectively acquiring all points of the seismic data in the three-dimensional stereo window, thereby acquiring the attitude angle and the point inclination angle of the discrete data points of the three-dimensional seismic target range.
9. 9. The clastic rock deposition profile restoration method according to claim 1, wherein step S80 includes the steps of: S810, obtaining the point inclination angle and azimuth angle of the target restoration interface by using the coordinates and time points (L i ,T j ,Z k ) of the corresponding three-dimensional earthquake of the target restoration interface, and recording as (L i ,T j ,Z k , , ) Where i is a line direction number, j is a track direction number, k is a longitudinal depth sampling point number, L i is an i-th line number, T j is a j-th track number, Z k is a depth value, For the dot tilt angle of the spatial ijk position, Azimuth for spatial ijk position; S820, the coordinate and time point (L i ,T j ,Z k ) of the three-dimensional earthquake are corresponding to the micro-drift plane, the point inclination angle and azimuth angle of the micro-drift plane are obtained and recorded as (L i ,T j ,Z k , , ) Where i is a line direction number, j is a track direction number, k is a longitudinal depth sampling point number, L i is an i-th line number, T j is a j-th track number, Z k is a depth value, For the dot tilt angle of the spatial ijk position, Is the azimuth of the spatial ijk position.
10. 10. The clastic rock deposition profile restoration method according to claim 1, wherein step S90 includes the steps of: s910, acquiring a sedimentary stratum calculation reference surface based on the following formula: Wherein: I is a line direction number, j is a track direction number, k is a longitudinal depth sampling point number, Subtracting all point time values on the micro-flyover for the spatial ijk position target restoration interface, The dot tilt angle of the interface data is restored for the spatial ijk position object, Restoring the azimuth of the interface data for the spatial point ijk position target, For the dot inclination angle of the spatial ijk position micro-drift plane data, Azimuth angle of the micro-drift surface data of the spatial ijk position; S920, starting from the maximum thickness, accumulating along the line direction by taking the upper super point as the boundary control point, and obtaining the accumulated thickness by using the following formula : Where i is the data point of the recovery interface from the target, And The number of superpoints experienced by the labeled target recovery interface for points i-1 and i respectively, And The stamp thickness of the target recovery interface for i-1 and i-point respectively, And The accumulated thickness of i-1 and i points respectively; S930, using median filtering to eliminate And (3) obtaining data after the stratum thickness is smoothed, wherein the data after the stratum thickness is smoothed is the impression thickness of the target recovery interface.
11. 11. The clastic rock deposition profile restoration method according to claim 1, wherein step S100 includes the steps of: s1100, obtaining a contour map of the thickness of the stamp by discrete points of the thickness of the stamp based on an inverse distance weighted interpolation method; s1200, dividing the sedimentary palace according to the thickness change characteristics and the depth change rule of the clastic rock water body, and recovering the clastic rock sedimentary palace.
12. 12. A clastic rock deposition morphology restoration device, comprising: The work area establishing module is used for establishing a clastic rock deposition landform interpretation and research work area; the first interface acquisition module is used for acquiring the seismic stratum interface of the target recovery interface; The earthquake layering acquisition module is used for acquiring the interfacial earthquake layering of the stratum covered by the target stratum; The second interface acquisition module is used for acquiring the minimum isochronous earthquake stratum interface of the stratum covered by the target stratum; the super point number acquisition module is used for analyzing based on the same-phase-axis contact relation of the earthquake to acquire the super point number of the target recovery interface; the micro-drift amount acquisition module is used for acquiring the micro-drift amount of the micro-drift surface; The first azimuth angle and point inclination angle acquisition module is used for acquiring the attitude angle and the point inclination angle of discrete data points of the three-dimensional seismic target range; The second dot inclination angle and azimuth angle acquisition module is used for acquiring the dot inclination angle and azimuth angle of the target recovery interface and the micro-drift surface; The die thickness calculation module is used for calculating the die thickness of the target recovery interface; And the clastic rock deposition topography recovery module is used for recovering the clastic rock deposition topography.
13. 13. A computer readable storage medium comprising instructions which, when run on a computer, cause the computer to perform the clastic rock deposition profile restoration method of any one of claims 1-11.
14. 14. A computing device, the computing device comprising: At least one processor, memory, and input output unit; The memory is used for storing a computer program, and the processor is used for calling the computer program stored in the memory to execute the clastic rock deposit topography restoration method according to any one of claims 1-11.

Description

Clastic rock deposition landform restoration method, device, medium and computing equipment Technical Field The invention relates to the technical field of oil and gas exploration, in particular to a clastic rock sedimentary landform restoration method, a device, a medium and computing equipment. Background The study of clastic rock sedimentary reservoirs is quite complex, macroscopically, the structure controls the sediments to reflect the state of ancient structures, and the matching relationship must be clarified to define the intermediate connection element-ancient landform. The paleo-topography and the space of the source region distribution match the configuration relationship, thereby controlling the deposition sequence and the filling mode of the deposition. For oil and gas workers, the description of the landform, in particular the accurate depiction of the landform, can directly indicate the regional paleo-water depth, and in the research process of clastic rock oil and gas reservoirs, the development positions of a sand reservoir, calcium-containing compact sandstone and hydrocarbon source rock can be comprehensively analyzed by utilizing the lithology longitudinal combination of drilling and logging through the depiction of the paleo-landform. Therefore, the accurate programming of the paleo-topography of the clastic rock has important significance for searching oil and gas reservoirs and oil and gas exploration and development well site deployment. At present, many ancient landform restoration methods exist, and the main working methods include a stamping method, a filling and leveling method, a sedimentology analysis method, a layer sequence stratigraphy method and the like, and the stamping method is the most commonly used method. The basic principle of the impression method is a deposition compensation principle, and the method comprises the steps of firstly selecting an overlying stratum as a reference surface, then calculating the thickness of the reference surface and a target restored original stratum, and representing the paleomorphic form through the calculated deposition thickness of the overlying stratum. However, the method cannot find the overlying adjacent isochronous surface as the reference surface, so that the disadvantage of being unable to accurately track the isochronous reference surface is caused, and the sedimentary topography of the clastic rock cannot be completely recovered. Disclosure of Invention The invention mainly aims to provide a clastic rock deposition landform restoration method, a device, a medium and computing equipment, so as to solve the technical problem that the clastic rock deposition landform cannot be fully restored in the prior art. The invention provides a clastic rock sedimentary topography restoration method, which comprises the following steps of S10, establishing a clastic rock sedimentary topography interpretation research work area, S20, obtaining an earthquake stratum interface of a target restoration interface, S30, obtaining an interface earthquake layering of an overlying stratum of the target stratum, S40, obtaining a minimum isochronous earthquake stratum interface of the overlying stratum of the target stratum, S50, analyzing based on a same-phase-axis contact relation of earthquakes, obtaining the number of times of superpoints of the target restoration interface, S60, obtaining micro-floating surface micro-floating quantity, S70, obtaining an attitude angle and a point inclination angle of discrete data points of a three-dimensional earthquake target range, S80, obtaining the point inclination angle and the azimuth angle of the target restoration interface, S90, calculating the impression thickness of the target restoration interface, and S100, and restoring the clastic rock sedimentary topography. In some embodiments, the step S10 comprises the steps of S110, based on seismic data of a research area, obtaining a channel number and a plane geodetic coordinate of an earthquake, obtaining a three-point coordinate of the earthquake, obtaining a corresponding relation between the channel number and the geodetic coordinate, thereby establishing a clastic rock sedimentary relief interpretation research area, loading seismic data, S120, based on a seismic profile, analyzing continuous variation conditions of a phase axis of a seismic reflection wave, selecting continuous data, loading the continuous data as the research data into the clastic rock sedimentary relief interpretation research area, S130, obtaining wells in the three-dimensional work area, drilling wells meeting a target and recovering stratum interfaces, and selecting well data which are uniformly distributed in the three-dimensional area and have complete logging curves as other research data, and loading the well data into the clastic rock sedimentary relief interpretation research area. In some embodiments, the step S20 comprises the steps of S210, S220, analyzing the seismic