AU2022399667A1 - Method and apparatus for inversion of crustal structure of passive continental margin based on subsidence in stratigraphic record - Google Patents
Method and apparatus for inversion of crustal structure of passive continental margin based on subsidence in stratigraphic record Download PDFInfo
- Publication number
- AU2022399667A1 AU2022399667A1 AU2022399667A AU2022399667A AU2022399667A1 AU 2022399667 A1 AU2022399667 A1 AU 2022399667A1 AU 2022399667 A AU2022399667 A AU 2022399667A AU 2022399667 A AU2022399667 A AU 2022399667A AU 2022399667 A1 AU2022399667 A1 AU 2022399667A1
- Authority
- AU
- Australia
- Prior art keywords
- crustal
- basement
- stratigraphic
- subsidence
- density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 241001191378 Moho Species 0.000 claims abstract description 33
- 238000005553 drilling Methods 0.000 claims abstract description 18
- 238000004590 computer program Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 abstract description 10
- 230000015572 biosynthetic process Effects 0.000 abstract description 8
- 230000007246 mechanism Effects 0.000 abstract description 5
- 101100457415 Rattus norvegicus Mmp23 gene Proteins 0.000 abstract 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 abstract 1
- 239000003795 chemical substances by application Substances 0.000 abstract 1
- 239000004148 curcumin Substances 0.000 abstract 1
- 230000006870 function Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000007704 transition Effects 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
J 9hW(1-NiJPM tIN J z EP li l
(1)19|)i
() Pd VT, WO 2023/098441 A1
2023 4 6 )P 8 (08.06.2023) W IPO I PCT
(51) M pj9 41J)}% A : ( 72 ) L : 5& (ZHAO, Zhongxian); F 3
GO1VJ/30 (2006.01) fil 4$ ~ M 6 X *A Ml1119 ,Guangdong 511458
()%¶fl119t,Guangdong 511458 (CN)o
(22) M pj: $1i ( : 2022 411 ) 10 H (10.11.2022) (7){ii1i tu]gogi1458 flNRQ
(74) JtTT A : f- 1,)tl f§4 51 -! 1]i M~ $f it TT phR q]
(25) $+ (GUANGZHOU KEYUE I.P. LAW OFFICE); +I
(26) Qt : [ 3T j F r A A E100 rx23-1
(30)Vt -t616 Guangdong 510070 (CN)o
202210951025.4 202248)A9H (09.08.2022) CN (81)1"(R MMHA,
gtr): AE, AG, AL, AM, AO, AT, AU, AZ, BA, BB, BG,
(71) A 44p~n;t t4 iTi-M R Fi(SOUTH CHI- BH, BN, BR, BW, BY, BZ, CA, CH, CL, CN, CO, CR, CU,
NA SEA INSTITUTE OF OCEANOLOGY, CHINESE CV, CZ, DE, DJ, DK, DM, DO, DZ, EC, EE, EG, ES, FI,
ACADEMY OF SCIENCES) [CN/CN]; + GB, GD, GE, GH, GM, GT, HN, HR, HU, ID, IL, IN, IQ, )'I T1 )+1 M 61 X * 1119 , Guangdong IR, IS, IT, JM, JO, JP, KE, KG, KH, KN, KP, KR, KW, KZ,
511458 (CN)° LA, LC, LK, LR, LS, LU, LY, MA, MD, MG, MK, MN,
(54) Title: METHOD AND APPARATUS FOR INVERTING CRUSTAL STRUCTURE OF PASSIVE CONTINENTAL MARGIN
ON BASIS OF SUBSIDENCE IN STRATIGRAPHIC RECORD
AA BID cc
DD
AA Deep reflection seism
1E 313lhtS~h~d]B Oil andogas driling
CC Ocean drilling
EE FF DD Shallow stratigraphic reflection interface, stratigraphic
framework and stratigraphic depth profile; sedimentary
basement reflection and Moho reflection
- fEE Actually measured total tectonic subsidence amount
of abasement
r FF Crustal round-trip reflection time
GG III h GG Assign the value of an initial crustal average density
HH According to a speed-density empirical formula,
calculate a crustal average speed and thickness and the
HH depth of a Moho surface
|| Increase the crustal average density by 1 kg/mA3
JJ Theoretical total tectonic subsidence amount of the
basement
1J KK Is theoretical subsidence consistent with actually
I
kL measured subsidence?
ILL No
KKMM Yes
KK NN Crustal density, speed and thickness, and depth of
ii. the Moho
NN MM
(57) Abstract: A method and apparatus for inverting a crustal structure of a passive continental margin on the basis of subsidence in
ena stratigraphic record. The method comprises seven steps. In comparison with the prior art, without the need to additionally provide
gravity or refraction seismic data and simply by means of reflection information of a multichannel seism to a stratum, a sedimentary
basement and the Moho, the crustal average density, speed and thickness can be quickly and efficiently calculated within a large range
W O 2023/098441 A 1 |||||111I||III||I|||||I||||||I|||||||||||||i|iD |||||||||||||||||||||||||||
MW, MX, MY, MZ, NA, NG, NI, NO, NZ, OM, PA, PE,
PG, PH, PL, PT, QA, RO, RS, RU, RW, SA, SC, SD, SE,
SG, SK, SL, ST, SV, SY, TH, TJ, TM, TN, TR, TT, TZ,
UA, UG, US, UZ, VC, VN, WS, ZA, ZM, ZWc
(8 4 ) - -M(j hrH R P ) fAln -¶E
M': ARIPO (BW, CV, GH, GM, KE, LR, LS, MW,
MZ, NA, RW, SD, SL, ST, SZ, TZ, UG, ZM, ZW), RlE
(AM, AZ, BY, KG, KZ, RU, TJ, TM), [III (AL, AT, BE,
BG, CH, CY, CZ, DE, DK, EE, ES, Fl, FR, GB, GR, HR,
HU, IE, IS, IT, LT, LU, LV, MC, ME, MK, MT, NL, NO,
PL, PT, RO, RS, SE, SI, SK, SM, TR), OAPI (BF, BJ, CF,
CG, CI, CM, GA, GN, GQ, GW, KM, ML, MR, NE, SN,
TD, TG).
kJBAiR Mifr e( WJ48. 2 (h))
- tfiEAQ1t$, i #21& (2) (a) fTfjl
to establish a crustal attribute and a depth structure. An effective means is provided for obtaining space variation information of a crustal
structure of a passive continental margin that lacks data, and an important basis is thus provided for understanding the formation, basin
formation and reservoir formation mechanisms of the passive continental margin.
(57) rjR: -fg ATh ji-AAt h ft 'j7ht,-W JkA tHM
~~AttS~~~al'~ tAtAtT-' LL, TTiAhTh
Description
[0001] The present invention relates to the field of crustal structures, and in particular to a method and apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record.
[0002] Passive continental margins are stable transition zones between continents and oceans, extending about 105,000 km around the world, about 2.6 times the circumference of the earth. They are the most important areas for sedimentation and mineral resource (such as oil and gas) accumulation on the earth's surface, providing an important guarantee for the sustainable economic and social development of mankind. The passive continental margins are formed through extension, thinning and rupture of the continental lithosphere, showing various structural types with different magmatic, metamorphic and sedimentary processes, which significantly affect the evolution of continental margin basins and the formation of oil and gas resources. Therefore, the crustal structures of the passive continental margins are of great significance for understanding the continental margin extension, thinning process, continent-ocean transition, basin formation mechanism and oil and gas evaluation.
[0003] The crustal structure of the continental margin is explored mainly on basis of drilling, submarine seism and multi-channel reflection seism. However, offshore drilling and submarine seism are expensive and little data are acquired, which limits the understanding of the spatial characteristics of the crustal structure. The multi-channel reflection seism is relatively low in detection cost and readily available for high-density and high-accuracy measurements. The seismic reflection data can provide clear imaging for shallow strata, intracrustal reflection and deep Moho, which is an effective method for understanding the deep crustal reflection structure. However, it is still difficult to determine information such as crustal attributes and thickness based on seismic reflection alone, and it has to combine other data for further analysis. For example, under the constraints of multi-channel seism on the shallow strata, gravity data is often combined to carry out gravity-seismic simulation of the density attribute and thickness of the crust, or refraction seismic travel time information is combined to carry out ray tracing for simulation of the velocity attribute and thickness of the crust. Obviously, these methods not only require additional data and increase the cost, but also show simplicity of simulated crustal attribute and possible multiplicity of solution of simulation results.
[0004] Chinese Patent No. CN113740915A discloses a method for joint inversion of parameters of a crustal structure through gravity and a receiving function of a spherical coordinate system. In this method, the gravity data and the receiving function data are fitted at the same time in an inversion process, and the complementation effect of the gravity and the receiving function is achieved by a joint inversion algorithm, so that the multiplicity of solution of the single data inversion is reduced. Moreover, the influence of the earth curvature is considered in the joint inversion, and a forward modeling method of a Tesseroid unit body in the spherical coordinate system is introduced. In this method, the high resolution of the gravity in the transverse direction and the high resolution of the receiving function in the depth direction near a station are considered, so that the more accurate crustal structure parameters are obtained. However, this method still requires additional data.
[0005] The present invention provides a rapid and efficient method and apparatus for inversion of crustal density, velocity structure and thickness of a passive continental margin based on the multi-channel reflection seismic interpretations of the strata, sedimentary basement and crust mantle boundary (Moho) as well as the tectonic subsidence in a stratigraphic record in the absence of other data source such as gravity and refraction seismic data.
[0006] To achieve the object above, the technical solutions of the present invention are as follows.
[0007] In a first aspect, the present invention provides a method for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record. The method includes:
[0008] step 1: carrying out stratigraphic correlation and division as well as interpretations of a sedimentary basement and Moho according to multi-channel seismic data under constraints of oil and gas drilling data or ocean drilling data, establishing a shallow stratigraphic framework, and performing time-depth conversion to obtain a stratigraphic depth profile;
[0009] step 2: carrying out backstripping analysis according to the principle of Airy isostasy under constraints of well lithology to backstrip all strata and water above the sedimentary basement, and then performing eustasy correction to obtain an observed total basement tectonic subsidence amount;
[0010] step 3: calculating a crustal two-way reflection travel time between the sedimentary basement and the Moho according to the interpretations of the sedimentary basement and the crust-mantle boundary Moho along the multi-channel seismic profile;
[0011] step 4: assigning an initial value to a mean crustal density;
[0012] step 5: calculating a mean crustal velocity according to an empirical formula of a velocity-density relation, and then calculating a crustal thickness and a depth of the Moho surface using the mean crustal velocity and the two-way travel time;
[0013] step 6: calculating a corresponding theoretical total basement tectonic subsidence amount according to the crustal thickness and the crustal density; and
[0014] step 7: comparing the calculated total basement tectonic subsidence amount and the observed total basement tectonic subsidence amount, if the two are identical, considering that the assigned crustal density is reasonable, and if the two are non-identical, increasing the mean crustal density by a set value each time and repeating steps 5-7 until a result meeting an accuracy requirement is obtained.
[0015] In a second aspect, the present invention provides an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record. The apparatus includes:
[0016] a stratigraphic depth profile module configured to perform analysis, carry out stratigraphic correlation and division as well as interpretations of a sedimentary basement and Moho according to multi-channel seismic data under constraints of oil and gas drilling data or ocean drilling data, establish a shallow stratigraphic framework, and perform time-depth conversion to obtain a stratigraphic depth profile;
[0017] an observed total basement tectonic subsidence amount module configured to carry out backstripping analysis according to the principle of Airy isostasy under constraints of well lithology to backstrip all strata and water above a sedimentary basement, and then perform eustasy correction to obtain an observed total basement tectonic subsidence amount;
[0018] a crustal two-way reflection travel time module configured to calculate a crustal two-way reflection travel time between the sedimentary basement and the Moho according to the interpretations of the sedimentary basement and the Moho along the multi-channel seismic profile;
[0019] an initial value module configured to assign an initial value to a mean crustal density;
[0020] a calculation module configured to calculate a mean crustal velocity according to an empirical formula of a velocity-density relation, and calculate a crustal thickness and a depth of the Moho surface using the mean crustal velocity and the two-way travel time;
[0021] a theoretical total basement tectonic subsidence amount module configured to calculate a corresponding theoretical total basement tectonic subsidence amount according to the crustal thickness and the crustal density; and
[0022] a comparison module configured to compare the calculated total basement tectonic subsidence amount and the observed total basement tectonic subsidence amount; if the two are identical, consider that the assigned crustal density is reasonable; and if the two are non-identical, increase the mean crustal density by a set value each time and execute repeatedly the calculation module, the theoretical total basement tectonic subsidence amount module and the comparison module until a result meeting an accuracy requirement is obtained.
[0023] In a third aspect, the present invention provides an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in stratigraphic record. The apparatus includes a memory, a processor, and a computer program stored in the memory and runnable on the processor, wherein the processor, when executing the computer program, implements the steps of the method defined as above.
[0024] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the method described as above.
[0025] Compared with the prior art, the present invention has the following beneficial effects.
[0026] The present invention can rapidly and efficiently establish crustal attributes and depth structures in a wide range through calculating the mean crustal density, velocity and thickness only according to seismic reflection data, without additional provision of gravity or refraction seismic data. The present invention provides an effective means for acquiring information on the spatial variation of the crustal structure of passive continental margins where lack data, and also provides an important basis for understanding the formation mechanisms of the passive continental margins as well as the sedimentary basins and hydrocarbon resources.
[0027] FIG. 1 shows a flowchart of a method for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record according to Embodiment 1;
[0028] FIG. 2 shows a schematic diagram of establishing a crustal structure depth profile by inversion of a crustal density, velocity attribute and thickness based on tectonic subsidence in a stratigraphic record according to the interpretations of shallow strata, a sedimentary basement and Moho along the multi-channel seismic profile;
[0029] FIG. 3 shows a schematic diagram of the composition of an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record according to Embodiment 2, and
[0030] FIG. 4 shows a schematic diagram of the composition of an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record according to Embodiment 3.
[0031] The technical solutions of the present invention will be further explained below in conjunction with the attached drawings and the embodiments.
[0032] Embodiment 1
[0033] The continent-ocean transition process of the passive continental margin in the South China Sea is the focus and frontier of current ocean drilling and international geoscientific researches. In the ocean-continent transition zone of the northern continental margin of the South China Sea, multi-channel seism clearly images the reflection in shallow strata, a sedimentary basement and deep crust-mantle boundary. Here, a specific process of applying a method for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record is illustrated by taking this profile as an embodiment. As shown in FIG. 1 and FIG. 2, the method specifically includes the following steps.
[0034] In step 1, stratigraphic correlation and division as well as reflection interpretations of a sedimentary basement and a crust-mantle boundary (Moho) (FIG. 2a) are first carried out according to multi-channel seismic reflection data and in combination with existing rock core data of oil and gas drilling and ocean drilling; a shallow stratigraphic framework is established; and time-depth conversion is performed to obtain a stratigraphic depth profile (FIG. 2f).
[0035] In step 2, backstripping analysis is carried out on shallow strata according to the principle of Airy isostasy by utilizing stratigraphic lithology information revealed through drilling or drilling to backstrip all strata and water above the sedimentary basement; then a difference value between an initial rifting and present eustasy is corrected; and an observed total basement tectonic subsidence amount is calculated (FIG. 2e).
[0036] In step 3, a crustal two-way reflection travel time between the sedimentary basement and the Moho is calculated according to the interpretations of the sedimentary basement and the Moho along the multi-channel seismic profile (FIG. 2a).
[0037] In step 4, an initial value of 1500 kg/m3 is assigned to a mean crustal density.
[0038] In step 5, a mean crustal velocity (km/s) is then calculated according to an empirical formula 1 of a velocity (V)-density (D) relation, and a crustal thickness and a depth of the Moho surface are calculated using the mean crustal velocity and the two-way travel time:
[00391 V=0.002831-D-1.593994 (1).
[0040] In step 6, a corresponding theoretical total basement tectonic subsidence amount is calculated according to the crustal thickness and the crustal density.
[0041] In step 7, the calculated total basement tectonic subsidence amount and the observed total basement tectonic subsidence amount are compared; if the two are identical (FIG. 2e), the assigned crustal density is considered reasonable and thus corresponding crustal structure information such as the crustal density (FIG. 2b), velocity (FIG. 2c), thickness (FIG. 2d) and Moho depth (FIG. 2f) are obtained; and if the two are non-identical, the mean crustal density is increased by 1 kg/m3 each time, and steps 5-7 are repeated until a result meeting an accuracy requirement is obtained.
[0042] It thus can be seen that, compared with the prior art, the method of the present invention can rapidly and efficiently establish crustal attributes and depth structures in a wide range by calculating the mean crustal density, velocity and thickness only according to seismic reflection data, without additional provision of gravity or refraction seismic data. The present invention provides an effective means for acquiring information on the spatial variation of the crustal structure of the passive continental margin that lacks data, and then provides an important basis for understanding the formation mechanisms of the passive continental margin as well as the sedimentary basins and hydrocarbon resources.
[0043] Embodiment 2
[0044] Referring to FIG. 3, this embodiment provides an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record. The apparatus includes:
[0045] a stratigraphic depth profile module configured to perform analysis, carry out stratigraphic correlation and division as well as interpretations of a sedimentary basement and Moho according to multi-channel seismic data under constraints of oil and gas drilling data or ocean drilling data, establish a shallow stratigraphic framework, and perform time-depth conversion to obtain a stratigraphic depth profile;
[0046] an observed total basement tectonic subsidence amount module configured to carry out backstripping analysis according to the principle of Airy isostasy under constraints of well lithology to backstrip all strata and water above the sedimentary basement, and then perform eustasy correction to obtain an observed total basement tectonic subsidence amount;
[0047] a crustal two-way reflection travel time module configured to calculate a crustal two-way reflection travel time between the sedimentary basement and the Moho according to the interpretations of the sedimentary basement and the Moho along the multi-channel seismic profile;
[0048] an initial value module configured to assign an initial value to a mean crustal density, which is 1500 kg/m3 in this embodiment;
[0049] a calculation module configured to calculate a mean crustal velocity according to an empirical formula of a velocity-density relation, and calculate a crustal thickness and a depth of the Moho surface using the mean crustal velocity and the two-way travel time;
[0050] a theoretical total basement tectonic subsidence amount module configured to calculate a corresponding theoretical total basement tectonic subsidence amount according to the crustal thickness and the crustal density; and
[0051] a comparison module configured to compare the calculated total basement tectonic subsidence amount and the observed total basement tectonic subsidence amount; if the two are identical, consider the assigned crustal density reasonable; and if the two are non-identical, increase the mean crustal density by 1 kg/m3 each time and execute repeatedly the calculation module, the theoretical total basement tectonic subsidence amount module and the comparison module until a result meeting an accuracy requirement is obtained.
[0052] The present invention can rapidly and efficiently establish crustal attributes and depth structures in a wide range by calculating mean crustal density, velocity and thickness only according to seismic reflection information of multi-channel seism to strata, the sedimentary basement and the Moho, without additional provision of gravity or refraction seismic data. The present invention provides an effective means for acquiring information on the spatial variation of the crustal structure of the passive continental margin that lacks data, and then provides an important basis for understanding the formation mechanisms of the passive continental margin as well as the sedimentary basins and hydrocarbon resources.
[0053] Embodiment 3:
[0054] Referring to FIG. 4, this embodiment provides an apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record. The apparatus includes a processor, a memory, and a computer program stored in the memory and runnable on the processor. For example, the computer program is a processing program for the method for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record. The processor, when executing the computer program, implements the steps of Embodiment 1 above, for example, as shown in FIG. 1. Alternatively, the processor, when executing the computer program, achieves the function of various modules in Embodiment 2 above.
[0055] Exemplarily, the computer program may be divided into one or more modules/units, which are stored in the memory and executed by the processor to achieve the present invention. The one or more modules/units may be a series of computer program instruction segments capable of completing a specific function, and the instruction segments are configured to describe the execution process of the computer program in the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record.
[0056] The apparatus for inversion of the crustal structure of the passive continental margin based on the subsidence in the stratigraphic record may be a computing device such as a desktop computer, a notebook, a handheld computer, and a cloud server. The apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record may include, but not limited to, a processor and a memory. It can be understood by those skilled in the art that FIG. 4 only shows an example of the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record, and is not intended to limit the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record, and the apparatus may include more or fewer components than those as illustrated, or a combination of some components, or components different from those as illustrated. For example, the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record may further include an input/output device, a network access device, a bus, or the like.
[0057] The processor may be a central processing unit (CPU) or otherwise, a general-purpose processor, a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, or discrete hardware components or the like. The general-purpose processor may be a microprocessor, or any conventional processor or the like.
[0058] The memory may be an internal storage unit, for example a hard disk or an internal memory, of the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record. The memory may also be an external storage device, for example an installed plug-in hard disk, smart media card (SMC), a secure digital (SD) card, a flash card or the like, of the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record. Further, the memory may further include both the internal storage unit and the external storage device of the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record. The memory is configured to store the computer program and other programs and data required by the apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record. The memory may also be configured to temporarily store data that have been output or are to be output.
[0059] Embodiment 4
[0060] This embodiment provides a computer-readable storage medium storing a computer program. The computer program, when executed by a processor, implements the steps of the method described in Embodiment 1.
[0061] The computer-readable medium as shown may be any apparatus that may contain, store, communicate, propagate, or transmit programs for use in or in conjunction with an instruction execution system, apparatus, or device. More specific examples (non-exhaustive listing) of the computer-readable medium include: an electrical connection (electronic apparatus) with one or more arranged wires, a portable computer cassette (magnetic apparatus), a random-access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber apparatus, and a portable compact disc read-only memory (CD-ROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program may be printed. For example, the paper or other medium is optically scanned, and then, editing, interpreting or, if necessary, other appropriate processing are performed to obtain the program electronically, which is then stored in a computer memory.
[0062] The above embodiments are merely for illustrating the technical concepts and features of the present invention, and are intended to allow those of ordinary skills in the art to understand and thereby implement the content of the present invention, and the protection scope of the present invention cannot be limited thereto. Any equivalent changes or modifications made according to the spirit of the present invention should be construed as falling within the protection scope of the present invention.
Claims (10)
- CLAIMS What is claimed is: 1. A method for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record, characterized by comprising: step 1: carrying out stratigraphic correlation and division as well as interpretations of a sedimentary basement and a crust-mantle boundary Moho according to multi-channel seismic data under constraints of oil and gas drilling data or ocean drilling data, establishing a shallow stratigraphic framework, and performing time-depth conversion to obtain a stratigraphic depth profile; step 2: carrying out backstripping analysis according to the principle of Airy isostasy under constraints of well lithology to backstrip all strata and water above the sedimentary basement, and then performing eustasy correction to obtain a measured total tectonic subsidence amount of the basement; step 3: calculating a crustal two-way reflection travel time between the sedimentary basement and the Moho according to the interpretations of the sedimentary basement and the Moho along the multi-channel seismic profile; step 4: assigning an initial value to a mean crustal density; step 5: calculating a mean crustal velocity according to an empirical formula of a velocity density relation, and calculating a crustal thickness and a depth of Moho using the mean crustal velocity and a two-way travel time; step 6: calculating a corresponding theoretical total tectonic subsidence amount of the basement according to the crustal thickness and the crustal density; and step 7: comparing the calculated total tectonic subsidence amount of the basement with the measured total tectonic subsidence amount of the basement, if the two are identical, considering that the assigned crustal density is reasonable, and if the two are non-identical, increasing the mean crustal density by a set value each time and repeating steps 5-7 until a result meeting an accuracy requirement is obtained.
- 2. The method for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 1, characterized in that the empirical formula of the velocity V-density D relation is as follows: V = 0.002831- D -1.593994
- 3. The method for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 1, characterized in that in step 7, the set value is 1 kg/m3
- 4. The method for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 1, characterized in that the initial value is 1500 kg/m3
- 5. An apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record, comprising: a stratigraphic depth profile module configured to perform analysis, carry out stratigraphic correlation and division as well as interpretations of a sedimentary basement and Moho according to multi-channel seismic data under constraints of oil and gas drilling data or ocean drilling data, establish a shallow stratigraphic framework, and perform time-depth conversion to obtain a stratigraphic depth profile; a measured total basement tectonic subsidence amount module configured to carry out backstripping analysis according to the principle of Airy isostasy under constraints of well lithology to backstrip all strata and water above the sedimentary basement, and then perform eustasy correction to obtain a measured total basement tectonic subsidence amount; a crustal two-way reflection travel time module configured to calculate a crustal two-way reflection travel time between the sedimentary basement and the Moho according to the interpretations of the sedimentary basement and the Moho along the multi-channel seismic profile; an initial value module configured to assign an initial value to a mean crustal density; a calculation module configured to calculate a mean crustal velocity according to an empirical formula of a velocity-density relation, and calculate a crustal thickness and depth of Moho using the mean crustal velocity and the two-way travel time; a theoretical total basement tectonic subsidence amount module configured to calculate a corresponding theoretical total basement tectonic subsidence amount according to the crustal thickness and the crustal density; and a comparison module configured to compare the calculated total basement tectonic subsidence amount and the measured total basement tectonic subsidence amount; if the two are identical, consider that the assigned crustal density is reasonable; and if the two are non-identical, increase the mean crustal density by a set value each time and execute repeatedly the calculation module, the theoretical total basement tectonic subsidence amount module and the comparison module until a result meeting an accuracy requirement is obtained.
- 6. The apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 5, characterized in that the empirical formula of the velocity V-density D relation is as follows: V = 0.002831-D -1.593994.
- 7. The apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 5, characterized in that in the comparison module, the set value is 1 kg/m3 .
- 8. The apparatus for inversion of the crustal structure of the passive continental margin based on subsidence in the stratigraphic record according to claim 5, characterized in that the initial value is 1500 kg/m3 .
- 9. An apparatus for inversion of a crustal structure of a passive continental margin based on subsidence in a stratigraphic record, comprising a memory, a processor, and a computer program stored in the memory and runnable on the processor, characterized in that the processor, when executing the computer program, implements the steps of the method according to anyone of claims 1 to 4.
- 10. A computer-readable storage medium storing a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method according to anyone of claims 1 to 4.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210951025.4A CN115373024B (en) | 2022-08-09 | 2022-08-09 | Method and device for inverting passive land edge crustal structure based on stratum recording settlement |
CN202210951025.4 | 2022-08-09 | ||
PCT/CN2022/131035 WO2023098441A1 (en) | 2022-08-09 | 2022-11-10 | Method and apparatus for inverting crustal structure of passive continental margin on basis of subsidence in stratigraphic record |
Publications (1)
Publication Number | Publication Date |
---|---|
AU2022399667A1 true AU2022399667A1 (en) | 2024-02-22 |
Family
ID=84064344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
AU2022399667A Pending AU2022399667A1 (en) | 2022-08-09 | 2022-11-10 | Method and apparatus for inversion of crustal structure of passive continental margin based on subsidence in stratigraphic record |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN115373024B (en) |
AU (1) | AU2022399667A1 (en) |
WO (1) | WO2023098441A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117195511B (en) * | 2023-08-23 | 2024-04-30 | 中国科学院南海海洋研究所 | Quantitative calculation method for initial crust thickness and expansion coefficient |
CN117555025B (en) * | 2024-01-11 | 2024-04-02 | 应急管理部国家自然灾害防治研究院 | Multi-layer crust structure inversion method based on gravity data |
CN117572530B (en) * | 2024-01-17 | 2024-04-05 | 自然资源部第二海洋研究所 | Method for jointly determining ocean land boundary by gravity inversion Moholo surface and submarine earthquake |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2497576B (en) * | 2011-12-15 | 2017-06-14 | Statoil Petroleum As | ASEP+D Method:identifying anomalous areas of the earth's lower crust |
CN104459795A (en) * | 2014-12-08 | 2015-03-25 | 中国科学院南海海洋研究所 | Depth-varying-to-density earth crust extension coefficient thermal calibration gravity anomaly retrieval method |
CN106094019A (en) * | 2016-04-12 | 2016-11-09 | 中国石油化工股份有限公司 | The Depth Domain earth formation inversion method mapped based on geological information |
CN107589448B (en) * | 2017-07-13 | 2019-02-05 | 西安交通大学 | A kind of multitrace seismogram reflection coefficient sequence Simultaneous Inversion method |
CN110244352B (en) * | 2019-06-11 | 2020-12-25 | 西安石油大学 | Variable density-based crust thickness gravity inversion method |
CN110221344B (en) * | 2019-06-17 | 2020-08-28 | 中国地质大学(北京) | Seismic full-waveform and gravity joint inversion method for three-dimensional density structure of earth crust |
CN111337993A (en) * | 2020-03-30 | 2020-06-26 | 中国地质科学院 | Variable density and variable depth constraint-based gravity density interface inversion method |
CN111610561B (en) * | 2020-06-05 | 2021-05-07 | 中国地质大学(北京) | Calculation method for quantitatively establishing sea level change and sea advance-sea retreat curve |
CN112285771A (en) * | 2020-09-03 | 2021-01-29 | 田晓峰 | Hidden fracture imaging method in urban active fault detection |
CN113514900B (en) * | 2021-07-12 | 2022-05-17 | 吉林大学 | Density constraint-based spherical coordinate system gravity and gravity gradient joint inversion method |
-
2022
- 2022-08-09 CN CN202210951025.4A patent/CN115373024B/en active Active
- 2022-11-10 AU AU2022399667A patent/AU2022399667A1/en active Pending
- 2022-11-10 WO PCT/CN2022/131035 patent/WO2023098441A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN115373024A (en) | 2022-11-22 |
WO2023098441A1 (en) | 2023-06-08 |
CN115373024B (en) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2022399667A1 (en) | Method and apparatus for inversion of crustal structure of passive continental margin based on subsidence in stratigraphic record | |
CN109709603B (en) | Seismic horizon identification and tracking method and system | |
CN109800863B (en) | Logging phase identification method based on fuzzy theory and neural network | |
CA2890187C (en) | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks | |
CN105353412A (en) | Calculating method and system of well-to-seismic integration average speed field | |
US20140350905A1 (en) | Asep+d method: identifying anomalous areas of the earth's lower crust | |
US11169287B2 (en) | Method and system for automated velocity model updating using machine learning | |
CN107783185B (en) | A kind of processing method and processing device of tomographic statics | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
CN109188522B (en) | Velocity field construction method and device | |
CN107991705B (en) | Log curve correction method and apparatus based on Two-dimensional Statistical feature | |
CN108491482B (en) | Geological map dynamic synthesis method and system considering proximity relation | |
CN108646288B (en) | Method and device for establishing near-surface model | |
CN111638555B (en) | Method and device for explaining micro-logging | |
Xu et al. | Seismic singularity attribute and its applications in sub-seismic faults detection | |
CN113341461B (en) | Earthquake velocity prediction method, earthquake velocity prediction device and server | |
CN117094152B (en) | Basin sedimentation history simulation method and device for coupling sedimentation | |
Chen et al. | Seismic ahead-prospecting based on deep learning of retrieving seismic wavefield | |
CN116449422A (en) | Method, device, equipment and medium for establishing three-dimensional speed model of disc on positive fault | |
Guo et al. | GeoPDNN: A Semisupervised Deep Learning Neural Network Using Pseudolabels for Three-dimensional Urban Geological Modelling and Uncertainty Analysis from Borehole Data | |
CN118033736A (en) | Method, device, equipment, storage medium and product for constructing speed field model | |
CN114706126A (en) | Quantitative carving method and system for seismic geologic body under geological cognition constraint | |
Pei | Reservoir Data Modeling Technology Based on Well Seismic Information Fusion | |
Su et al. | Seismic prediction of porosity in tight reservoirs based on transformer | |
Wrona et al. | Complex fault system revealed by 3-D seismic reflection data with deep learning and fault network analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
DA3 | Amendments made section 104 |
Free format text: THE NATURE OF THE AMENDMENT IS: AMEND THE INVENTION TITLE TO READ METHOD AND APPARATUS FOR INVERSION OF CRUSTAL STRUCTURE OF PASSIVE CONTINENTAL MARGIN BASED ON SUBSIDENCE IN STRATIGRAPHIC RECORD |