CN110263763A - Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy - Google Patents
Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy Download PDFInfo
- Publication number
- CN110263763A CN110263763A CN201910615287.1A CN201910615287A CN110263763A CN 110263763 A CN110263763 A CN 110263763A CN 201910615287 A CN201910615287 A CN 201910615287A CN 110263763 A CN110263763 A CN 110263763A
- Authority
- CN
- China
- Prior art keywords
- lithology
- cyclostratigraphy
- resolution
- midpoint
- method based
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to sedimentary formation fields, are related to a kind of method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy.The method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy in the present invention demarcates the middle part point of different lithology unit on High-resolution Cyclostratigraphy section, using point in the middle part of different lithology unit, to indicate the cyclic fluctuation feature of primary deposit signal；The method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy in the present invention, the border issue of different lithology unit is evaded, to substance transfer and compaction of the undebatable lithologic boundary in diagenetic process, and the problem of can effectively solve the problem that a large amount of interference signals caused by the rectangle wave train based on lithologic boundary model foundation, the target astronomy signal in frequency domain has been highlighted, the simulation match result to primary deposit and astronomical signal is improved.
Description
Technical field
The invention belongs to sedimentary formation fields, are related to a kind of based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy
Method.
Background technique
Lithological cycle is lithology units alternating deposit different in sedimentary formation, recurrent phenomenon；Lithological cycle is
Cyclical stratigraphy important geologic climate substitution in avatar the most intuitive and cyclical stratigraphy research on geologic section
Index；Cycle compared to other geologic climate indexs (such as magnetic susceptibility, oxygen isotope) in stratigraphic section changes, lithological cycle
Variation more preferably can record and show palaeoenvironment, palaeoclimatic cyclic fluctuation；By to different lithology unit in High-resolution Cyclostratigraphy
Thickness change establishes the bed succession of variation of lithological as geologic climate index, is a kind of very convenient, effective High-resolution Cyclostratigraphy point
Analysis method；Compared with other geologic climate indexs (such as magnetic susceptibility, oxygen isotope), the foundation of lithology thickness sequence it is only necessary to
Field observation lithology, lithofacies variation and its variation on depth of stratum, directly acquisition lithology thickness data, do not need acquisition sample
Product and progress later period test, realize a kind of research process that is quick, pollution-free and not destroying ecological environment.
Existing lithology thickness sequence, is often based upon lithologic boundary model, assigns different numerical value to different lithology unit, than
Such as calcium carbonate content, high limestone is assigned a value of 1, and the low mud stone of calcium carbonate content is assigned a value of 0；It is demarcated on stratigraphic section different
The boundary point of lithology units, obtains the bed succession of formation thickness and lithologic boundary, is then existed with the value that the lithology is assigned
Homogenization interpolation (being all 1 or 0) is carried out inside lithology units, establishes the rectangle wave train in the domain of stratum based on lithologic boundary model
Column.Existing lithological cycle analysis method is carried out in spectrum analysis, dynamic Fu on the basis of acquiring rectangle wave train
The research methods such as leaf analysis, astronomical tune analysis, so that astronomical chronicle scale is established, the time of Accurate Calibration earth evolution history
Scale；But there is two o'clock in the boundary model of existing lithological cycle: one, spectrum analysis during, rectangle wave train
More Interference Peaks spectrum can be generated in frequency domain, increase the degree of difficulty of target astronomy signal in identification stratum；Two, it deposits
The palaeoenvironment paleo-climate change in period is continuous undulating changing pattern, and the rectangle wave train of existing flat-Catastrophe Model cannot
It is preferable to simulate primary deposit signal fluctuation characteristic at any time.
Summary of the invention
It is a kind of based on lithology midpoint Sine Interpolation model in order to solve the problems in the prior art, the purpose of the present invention is to providing
The method for analyzing High-resolution Cyclostratigraphy.
The specific technical solution of the present invention is as follows:
Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy, comprising the following steps:
Selection lithological cycle stratigraphic section simultaneously establishes initial layers position, carries out lithology units division, mark to lithological cycle stratigraphic section
The middle part point of fixed each lithology units and measure middle part point to the distance of initial layers position and to different types of lithology units into
Row assignment establishes two-dimensional coordinate point to the middle part point distance value and assignment of each lithology units, carries out according to two-dimensional coordinate point uniform
Sine Interpolation establishes lithology intermediate value sine wave sequence in the domain of stratum, carries out spectrum analysis and Fourier transformation to sine wave sequence
Analysis carries out astronomical tune analysis to sine wave sequence, draws, establishes astronomical chronicle scale.
The invention has the following advantages:
By being found after being analyzed practical High-resolution Cyclostratigraphy and theoretical model, the existing rectangle based on lithologic boundary model
In wave train, there are the interference signals within the scope of stronger high and low frequency, and the distribution form of target astronomy signal frequency is in square
It is not highlighted out in shape wave train, needs experience abundant and a variety of basic datas during actual analysis comprehensive
Close the frequency distribution range for judging echo signal.
The method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy in the present invention demarcates High-resolution Cyclostratigraphy section
The middle part point of upper different lithology unit, using point in the middle part of different lithology unit, to indicate the periodicity of primary deposit signal
Fluctuation characteristic；The method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy in the present invention, has evaded different lithology
The border issue of unit, thus substance transfer and compaction of the undebatable lithologic boundary in diagenetic process, and can
It effectively solves the problems, such as a large amount of interference signals caused by the rectangle wave train based on lithologic boundary model foundation, has highlighted frequency
Target astronomy signal in domain improves the simulation match result to primary deposit and astronomical signal.
Detailed description of the invention
Fig. 1 is the rectangle wave train of comparative example 1 of the present invention；
Fig. 2 is the sine wave sequence of the embodiment of the present invention 1；
Fig. 3 is the rectangular wave sequence spectrum analysis chart of comparative example 1 of the present invention；
Fig. 4 is the sine wave sequence spectrum analysis figure of the embodiment of the present invention 1；
Fig. 5 is the dynamic spectrum analysis figure of comparative example 1 of the present invention；
Fig. 6 is the dynamic spectrum analysis figure of the embodiment of the present invention 1；
Fig. 7 is the rectangle wave train of comparative example 2 of the present invention；
Fig. 8 is the sine wave sequence of the embodiment of the present invention 2；
Fig. 9 is the rectangular wave sequence spectrum analysis chart of comparative example 2 of the present invention；
Figure 10 is the sine wave sequence spectrum analysis figure of the embodiment of the present invention 2；
Figure 11 is the dynamic spectrum analysis figure of comparative example 2 of the present invention；
Figure 12 is the dynamic spectrum analysis figure of the embodiment of the present invention 2.
Specific embodiment
Comparative example 1
Selection has confirmed 1300cm thickness silicalite-mud stone High-resolution Cyclostratigraphy comprising slope and precessional cycle astronomy signal, will
SiO2The higher silicalite unit of content is assigned a value of 1, by SiO2The lower mud stone unit of content is assigned a value of -1, with the High-resolution Cyclostratigraphy
The bottom successively chooses 11 boundary points as initial layers position upwards, and distance using each boundary point to initial layers position is as each side
The X-coordinate of boundary's point establishes the seat of each boundary point using the lithology units assignment under each boundary point as the Y-coordinate of each boundary point
Punctuate (10.2,1), (18, -1), (26,1), (35, -1), (46.2,1), (47.8, -1), (50.2,1), (60.2, -1),
(65.8,1), (68.6, -1), (84,1)；Uniformity interpolation is carried out between above-mentioned each coordinate points, such as coordinate points (Xn, 1),
X value successively increases 0.1 when the interpolation of (Xm, -1), and Y value is unified for -1, (Xn, 1), (Xn+0.1, -1), (Xn+0.2, -1),
(Xn+0.3, -1) ... (Xm-0.1, -1), (Xm, -1) obtain rectangle wave train as shown in Figure 1, to the rectangle wave train
Spectrum analysis and Fourier transformation analysis are carried out, rectangular wave sequence spectrum analysis chart as shown in Figure 3 and as shown in Figure 5 is obtained
Dynamic spectrum analysis figure.
Embodiment 1
Selection has confirmed silicalite-mud stone High-resolution Cyclostratigraphy of the 1300cm thickness comprising slope and precessional cycle astronomy signal, will
SiO2The higher silicalite unit of content is assigned a value of 1, by SiO2The lower mud stone unit of content is assigned a value of -1, with the High-resolution Cyclostratigraphy
The bottom is successively chosen point in the middle part of 11 lithology units upwards, is put with each middle part to the distance of initial layers position as initial layers position
As the X-coordinate of each middle part point, lithology units assignment establishes each middle part point as the Y-coordinate of each middle part point where using middle part point
Coordinate points (5.1,1), (14.1, -1), (22,1), (30.5, -1), (40.6,1), (47, -1), (49,1), (55.2, -1),
(63,1), (67.2, -1), (76.3,1)；Uniformity interpolation is carried out between above-mentioned each coordinate points, such as coordinate points (Xn, 1),
X value successively increases 0.1 when the interpolation of (Xm, -1), amounts to insertion N=(Xm-Xn)/0.1 numerical value, Y value is from 1 to -1 according to just
String formulaF(t)=Acos ω tIt is gradually reduced, (Xn, 1), (Xn+0.1, (1-(cos (π * 0.1/0.1)/N+ π)+1))), (Xn+
0.2, (1-(cos (π * 0.2/0.1)/N+ π)+1))), (Xn+0.3, (1-(cos (π * 0.3/0.1)/N+ π)+1))) ...,
(Xm, -1) obtains sine wave sequence as shown in Figure 2, carries out spectrum analysis and Fourier transformation point to the sine wave sequence
Analysis, obtains sine wave sequence spectrum analysis figure and dynamic spectrum analysis figure as shown in FIG. 6 as shown in Figure 4.
It is compared by Fig. 3 and Fig. 4, Fig. 5 and Fig. 6 it is found that rectangle wave train shows its except 2 echo signal peak spectrums
His interference signal (disturbed one -2), interference signal is obvious in low frequency and high-frequency range；Sine wave sequence shows 2 target letters
Number peak spectrum, and apparent compacting (not interfering with 1,2 energy of interference are very low) has been carried out to the signal of low frequency and high-frequency range.
Comparative example 2
Selection synthesizes slope-precession of the equinoxes theoretical model sequence comprising slope and precessional cycle astronomy signal.The theoretical model sequence packet
Containing two column data of XY, the first column data X value is 0 Dao the 3000kyr time using 1kyr as interval, and the second column data Y value is every
The slope of 1kyr-precession of the equinoxes changing value, slope-precession of the equinoxes value fluctuating change between -4 to 5；The theoretical model sequence nucleotide sequence is converted
For ternary sequence of cycles, step includes: that the second column data Y≤0 is all assigned a value of 1, and 0<Y<2 item are assigned a value of 2, Y>=2 and are assigned a value of 3,
The cycle variation stratum for establishing unit 1, unit 2 and unit 3, simulates the High-resolution Cyclostratigraphy of three kinds of lithology units, with theory rotation
The stratum bottom is returned as initial layers position, successively chooses 11 boundary points upwards, is made with the distance of each boundary point to initial layers position
Each boundary is established using the lithology units assignment under each boundary point as the Y-coordinate of each boundary point for the X-coordinate of each boundary point
The coordinate points (1,2) of point, (8,3), (23,2), (37,1), (41,2), (54,3), (58,2), (83,1), (87,2), (99,
3), (101,2)；Uniformity interpolation is carried out between above-mentioned each coordinate points, such as coordinate points (Xn, 2), when the interpolation of (Xm, 3)
X value successively increases 0.1, and Y value is unified for 3, (Xn, 2), (Xn+0.1,3), (Xn+0.2,3), (Xn+0.3,3) ... (Xm-0.1,
3), (Xm, 3) obtains rectangle wave train as shown in Figure 7, carries out spectrum analysis and Fourier transformation to the rectangle wave train
Analysis, obtains rectangular wave sequence spectrum analysis chart as shown in Figure 9 and dynamic spectrum analysis figure as shown in figure 11.
Embodiment 2
Selection synthesizes slope-precession of the equinoxes theoretical model sequence comprising slope and precessional cycle astronomy signal.The theoretical model sequence packet
Containing two column data of XY, the first column data X value is 0 Dao the 3000kyr time using 1kyr as interval, and the second column data Y value is every
The slope of 1kyr-precession of the equinoxes changing value, slope-precession of the equinoxes value fluctuating change between -4 to 5；The theoretical model sequence nucleotide sequence is converted
For ternary sequence of cycles, step includes: that the second column data Y≤0 is all assigned a value of 1, and 0<Y<2 item are assigned a value of 2, Y>=2 and are assigned a value of 3,
The cycle variation stratum for establishing unit 1, unit 2 and unit 3, simulates the High-resolution Cyclostratigraphy of three kinds of lithology units, with theory rotation
The stratum bottom is returned as initial layers position, successively chooses 11 boundary points upwards, is put with each middle part to the distance work of initial layers position
For the X-coordinate of each middle part point, lithology units assignment establishes each middle part point as the Y-coordinate of each middle part point where using middle part point
Coordinate points (0.5,2), (4.5,3), (15.5,2), (30,1), (39,2), (47.5,3), (56,2), (70.5,1), (85,2),
(93,3), (100,2)；Uniformity interpolation is carried out between above-mentioned each coordinate points, such as coordinate points (Xn, 2), (Xm's, 3) is interleave
X value successively increases 0.1 when value, amounts to insertion N=(Xm-Xn)/0.1 numerical value, Y value is from 2 to 3 according to sine formulaF(t)=Acos ωtIt is gradually increased, (Xn, 2), (Xn+0.1, (2+(cos (π * 0.1/0.1)/N+ π)+1))/2), (Xn+0.2, (2+(cos (π *
0.2/0.1)/N+ π)+1))/2), (Xn+0.3,2+(cos (π * 0.3/0.1)/N+ π)+1))/2) ..., (Xm, 3), obtain as
Sine wave sequence shown in Fig. 8 carries out spectrum analysis to the sine wave sequence and Fourier transformation is analyzed, obtains such as Figure 10 institute
The sine wave sequence spectrum analysis figure and dynamic spectrum analysis figure as shown in figure 12 shown.
It is compared by Fig. 9 and Figure 10, Figure 11 and Figure 12 it is found that rectangle wave train is shown except 6 echo signal peak spectrums
Other interference signals (disturbed one, 2,3), interference signal is obvious in low frequency and high-frequency range；Sine wave sequence shows 6
Echo signal peak spectrum, and apparent compacting has been carried out to the signal of low frequency and high-frequency range and (has not interfered with 1,3, interferes 2 energy
It is very low).
Claims (1)
1. the method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy, which comprises the following steps: selection
Lithological cycle stratigraphic section simultaneously establishes initial layers position, carries out lithology units division to lithological cycle stratigraphic section, demarcates each lithology
The middle part point of unit and the distance for measuring middle part point to initial layers position carry out assignment to different types of lithology units, right
The middle part point distance value and assignment of each lithology units establish two-dimensional coordinate point, carry out homogeneous sinusoidal interpolation according to two-dimensional coordinate point,
Lithology intermediate value sine wave sequence in the domain of stratum is established, spectrum analysis is carried out to sine wave sequence and Fourier transformation is analyzed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910615287.1A CN110263763B (en) | 2019-07-09 | 2019-07-09 | Method for analyzing gyrus stratum based on lithology midpoint sine interpolation model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910615287.1A CN110263763B (en) | 2019-07-09 | 2019-07-09 | Method for analyzing gyrus stratum based on lithology midpoint sine interpolation model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110263763A true CN110263763A (en) | 2019-09-20 |
| CN110263763B CN110263763B (en) | 2021-06-29 |
Family
ID=67925187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910615287.1A Active CN110263763B (en) | 2019-07-09 | 2019-07-09 | Method for analyzing gyrus stratum based on lithology midpoint sine interpolation model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110263763B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024036819A1 (en) * | 2022-08-19 | 2024-02-22 | 中国石油天然气集团有限公司 | Chronostratigraphic-domain stratigraphic sedimentary cycle analysis method and apparatus |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1168476A (en) * | 1997-04-16 | 1997-12-24 | 北京大学 | Extra long electromagnetic wave remote detecting method and equipment |
| CN101349764A (en) * | 2007-07-20 | 2009-01-21 | 中国石油大学(北京) | Earthquake gyration analysis method |
| US20130063433A1 (en) * | 2010-05-27 | 2013-03-14 | Landmark Graphics Corporation | Method and system of rendering well log values |
| CN104237945A (en) * | 2014-09-18 | 2014-12-24 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic data self-adaptive high-resolution processing method |
| CN104765065A (en) * | 2015-04-17 | 2015-07-08 | 中国地质大学（北京） | Method and device for determining sedimentary stratum ancient landform |
| CN106526679A (en) * | 2017-01-10 | 2017-03-22 | 中国石油大学(华东) | Method for recovering non-uniform erosion thickness through utilizing cycle analysis method |
| CN109782343A (en) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | Stratigraphic cycles analysis method and device |
-
2019
- 2019-07-09 CN CN201910615287.1A patent/CN110263763B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1168476A (en) * | 1997-04-16 | 1997-12-24 | 北京大学 | Extra long electromagnetic wave remote detecting method and equipment |
| CN101349764A (en) * | 2007-07-20 | 2009-01-21 | 中国石油大学(北京) | Earthquake gyration analysis method |
| US20130063433A1 (en) * | 2010-05-27 | 2013-03-14 | Landmark Graphics Corporation | Method and system of rendering well log values |
| CN104237945A (en) * | 2014-09-18 | 2014-12-24 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic data self-adaptive high-resolution processing method |
| CN104765065A (en) * | 2015-04-17 | 2015-07-08 | 中国地质大学（北京） | Method and device for determining sedimentary stratum ancient landform |
| CN106526679A (en) * | 2017-01-10 | 2017-03-22 | 中国石油大学(华东) | Method for recovering non-uniform erosion thickness through utilizing cycle analysis method |
| CN109782343A (en) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | Stratigraphic cycles analysis method and device |
Non-Patent Citations (5)
| Title |
|---|
| GRAHAM P. WEEDON 等: "Cyclostratigraphy, stratigraphic gaps and the duration of the Hettangian Stage (Jurassic): insights from the Blue Lias Formation of southern Britain", 《GEOLOGICAL MAGAZINE》 * |
| 伊海生: "沉积旋回叠置形式的波形分析及旋回层序划分方法", 《沉积学报》 * |
| 姚旭: "东古特提斯洋大陆边缘二叠纪硅质岩成因研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 王贵文 等: "测井曲线谱分析方法及其在沉积旋回研究中的应用", 《石油勘探与开发》 * |
| 田世峰: "中、新生代旋回地层学研究及其油气地质意义-以惠州凹陷和煤山剖面为例", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024036819A1 (en) * | 2022-08-19 | 2024-02-22 | 中国石油天然气集团有限公司 | Chronostratigraphic-domain stratigraphic sedimentary cycle analysis method and apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110263763B (en) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hinnov et al. | Spectral analysis of the Middle Triassic Latemar limestone | |
| CN101923177A (en) | Velocity analysis method and device of multi-wave seismic data | |
| Li et al. | Cretaceous anomalous subsidence and its response to dynamic topography in the Songliao Basin, Northeast China | |
| CN105954802A (en) | Conversion method and apparatus of lithologic data volumes | |
| Bádenas et al. | Sequence stratigraphy and bedding rhythms of an outer ramp limestone succession (Late Kimmeridgian, Northeast Spain) | |
| Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
| CN109633593A (en) | A kind of Ground Penetrating Radar Signal quantitative analysis method and system | |
| Fischer et al. | Cyclostratigraphic approach to Earth’s history: an introduction | |
| CN110263763A (en) | Method based on lithology midpoint Sine Interpolation model analysis High-resolution Cyclostratigraphy | |
| CN106338764A (en) | Stratigraphic trap oil and gas reservoir overlap-denudation line identification method through bio-inspired computing | |
| Bond et al. | Evidence of astronomical forcing of the Earth's climate in Cretaceous and Cambrian times | |
| CN107507095A (en) | From fast conversion method from depth of stratum section to stratum time section and system | |
| Maurer et al. | Statistical time-series analysis and sedimentological tuning of bedding rhythms in a Triassic basinal succession (Southern Alps, Italy) | |
| CN106908838A (en) | The method for building target area stratum inelastic attenuation quality factor three-dimensional model | |
| CN103901475B (en) | One attribute isogram method for drafting and device | |
| Herkat | Application of correspondence analysis to palaeobathymetric reconstruction of Cenomanian and Turonian (Cretaceous) rocks of Eastern Algeria | |
| Wolak et al. | Reconstruction of depocenter evolution through time using relative stratigraphic thickness | |
| Yu et al. | Continuous wavelet transform and dynamic time warping-based fine division and correlation of glutenite sedimentary cycles | |
| CN110568149B (en) | Fine and rapid quantitative simulation method for hydrocarbon generation and discharge history of sedimentary basin hydrocarbon source rock | |
| CN106291690B (en) | Accumulate pattern the Forward Modeling in the stratum of orbital period tuning | |
| CN104062680A (en) | Method for calculating wave resistance and backstepping gradient of objective function | |
| CN106646608B (en) | A kind of layer model method for building up for observation system parameter demonstration | |
| Koyan et al. | Analyzing 3D multi-frequency ground-penetrating radar (GPR) data simulated across a realistic sedimentary model | |
| Olóriz et al. | The marl-limestone rhythmites from the Lower Kimmeridgian (Platynota Zone) of the central Prebetic and their relationship with variations in orbital parameters | |
| CN117195511B (en) | Quantitative calculation method for initial crust thickness and expansion coefficient |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |