CN110244352A - A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density - Google Patents
A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density Download PDFInfo
- Publication number
- CN110244352A CN110244352A CN201910502882.4A CN201910502882A CN110244352A CN 110244352 A CN110244352 A CN 110244352A CN 201910502882 A CN201910502882 A CN 201910502882A CN 110244352 A CN110244352 A CN 110244352A
- Authority
- CN
- China
- Prior art keywords
- density
- data
- thickness
- gravity anomaly
- crustal
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000005484 gravity Effects 0.000 claims abstract description 112
- 241001191378 Moho Species 0.000 claims abstract description 75
- 238000011160 research Methods 0.000 claims description 49
- 230000008859 change Effects 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 10
- 238000004513 sizing Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 5
- 238000005553 drilling Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000008021 deposition Effects 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 4
- 238000003475 lamination Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013019 agitation Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003653 coastal water Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012360 testing method 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. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The CRUSTAL THICKNESS gravitational inversion method based on variable density that the invention discloses a kind of obtains 1.0 CRUSTAL THICKNESS of Crustal Model CRUST and density data and density of upper mantle data that spatial resolution covering the whole world is 1 ° × 1 °；The shell curtain density contrast changed with horizontal position is established according to the upper, middle and lower CRUST DENSITY and thickness of Crustal Model CRUST 1.0 and density of upper mantle data；Gravity anomaly caused by Moho relief underneath changes is obtained by bouguer gravity anomaly；Gravity anomaly inversion Moho depth caused by being changed by Moho relief underneath simultaneously obtains CRUSTAL THICKNESS.Compared with prior art, the present invention establishes practical shell curtain density model using based on published data, technically easy to accomplish, and inversion result more meets practical geologic feature, and accuracy is higher.
Description
Technical field
The present invention relates to field of earth sciences, especially a kind of CRUSTAL THICKNESS gravitational inversion method based on variable density.
Background technique
The variation of earth gravitational field is unevenly distributed related, and substance with the material density from adjacent ground surface to Earth
Density Distribution is the reflection of earth system structure and variation, therefore, observation earth gravitational field and the distribution of Study of The Underground material density
Gravity changer (i.e. gravity anomaly) caused by uneven, to understanding earth's layers structure, geodynamics processes, resource and environment
Variation has important scientific meaning.In addition, gravimetric observation is simple, data acquisition is convenient, in contrast, gravimetric prospecting has warp
The advantages that Ji, depth of exploration are big and are quickly obtained information on area, therefore be widely used.
Gravity anomaly inversion has infrastructure research, energy resource investigation and earthquake prediction study etc. important
Value.According to calculating parameter.Difference, gravity anomaly inversion are divided into inversion for physical properties and body inverting.The former inverting underground space
The density of each discrete unit is determined the actual distribution of field source body by Density Distribution；The latter is inverting on the basis of given density
The characteristic point coordinate of underground space subdivision unit, to obtain the rolling shape of geologic body upper and lower interface.In recent years, different using gravity
Normal inverting Moho depth is further to obtain the main research that CRUSTAL THICKNESS has become gravitational inversion.
When using gravity anomaly inversion Moho depth and CRUSTAL THICKNESS, the density contrast of the earth's crust and earth mantle is a crucial ginseng
Number.Currently, scholar is using the progress inverting of normal density mode mostly.Liu Zuhui etc. (nineteen eighty-three) is utilizing the gravity anomaly inversion South Sea
When the Moho depth of sea area, earth's crust averag density takes 2.67g/cm3, outer mantle substance takes 3.27g/cm3, i.e., shell curtain density contrast be-
0.6g/cm3.It is 2.73g/cm that Xu Deqiong and Jiang Jiazhen (1989), which takes earth's crust averag density,3(shell curtain density contrast is -0.54g/
cm3) inverting Mid-northern South China Sea Moho depth, Tsai et al. (2005) take ground in inverting Northern Part of South China Sea Mohorovich discontinuity
The averag density of shell and outer mantle is respectively 2.84g/cm3And 3.28g/cm3, shell curtain density contrast is -0.44g/cm3, Li Shuling etc.
(2012) and Guan et al. (2016) also use the difference inverting of this density value South Sea Moho depth and east
The Moho depth in sea and the South Sea and its adjacent area, and it is -0.30g/cm that Hao Tian Yao etc. (2008), which takes shell curtain density contrast,3Inverting
Flow In Northeastern South China Sea Moho depth.Then, many scholars are in inverting Moho depth, take mostly shell curtain density contrast be-
0.50g/cm3(Liu Jianhua, 1993；Wang et al., 2017；Wu's trick ability etc., 2017；Unnikrishnan et
Al., 2018), difference is only that earth's crust averag density is different with outer mantle averag density value, as Liu Jianhua (1993) takes
The averag density of the earth's crust and outer mantle is respectively 2.82g/cm3And 3.32g/cm3, Unnikrishnan et al. (2018) takes
The averag density of the earth's crust and outer mantle is respectively 2.80g/cm3And 3.30g/cm3.In addition to the above numerical value, Trungetal. (2004
Year) take shell curtain density contrast for -0.40g/cm in inverting South Sea Moho depth3, Han Bo etc. (2014) also chooses density contrast
For -0.40g/cm3Inverting Taiwan and adjacent area Moho depth.Liand Wang (2016) is according to selection shell curtain density contrast
For -0.32g/cm3The inverting Moho depth of east Asia and region of Southeast.
Although above each normal density value solves the anti-of research area's Moho depth and CRUSTAL THICKNESS to a certain extent
Problem is drilled, but when inverting region is larger, it is inaccurate that shell curtain density contrast is portrayed merely with a constant value.For this purpose,
Some scholars simulate density using mathematical function with the variation of depth, and are incorporated in Moho depth inverting.Such as Feng Juan
Chen Yan etc. (2017) is also using referring to using exponential density-depth function inverting North China's Moho depth within (2014)
Number density-depth function inverting Subei Basin Moho depth；Utilize parabolic line density-depth letter within Zhang Enhui etc. (2015)
Count inverting Sichuan and Yunnan area Moho depth.In addition to this 2 kinds of density-depth functions, color density depth function is (such as multinomial
Formula) it also can be used for Moho depth inverting, but have not yet to see correlative study.In addition, in sea area Moho depth inverting,
For the accuracy for improving inverting, thermal agitation gravity anomaly can be corrected generally before inverting (Liand Wang, 2016
Year；Wu recruit just etc., 2017, Liu hold it is identical, 2018), to obtain accurate Mohorovich discontinuity gravity anomaly.
The above variable density function considers shell curtain density contrast along the variation of depth, more accurate than Chang Midu.But the density is only
Vertically change, when studying area and constructing complex (such as subduction zone, mid-oceanic ridge region), the variation of density transversely cannot
Ignore, for this purpose, utilizing the density in transverse direction when utilizing gravity interpretation the middle of the South China Sea Basin crustal structure within Yang Jinyu etc. (2011)
Variation combination has carried out Interpretation of profile.Liu Jun etc. (2009) is adopted using data such as seismic profile, sonobuoys as control point
Gravity profile fitting is carried out with lateral density variation, has obtained Northern Part of South China Sea crustal structure.In addition to sea area, in land-based area crustal structure
Changed when research using lateral density and carries out gravity profile fitting also as main research means.It establishes within Zhang Yongqian etc. (2014)
Middle and lower reaches of the Yangtze River and adjacent areas density structure of the crust；Wang Qianshen etc. (2015) has studied Erdos-Central Qinling-Si Chuandong
The deep-level rockburst in portion；It is had studied by the fitting of Lijing-Kweiyang gravity profile density within Zhang et al. (2018) high
The deep-level rockburst of Meishan Were Large Igneous Provinces.
The Lateral Change Characteristics that shell curtain density contrast is preferably reflected along the density contrast of profile change, with practical geologic feature
It more coincide, however it is only two-dimensional result, and zonal Moho depth and CRUSTAL THICKNESS variation characteristic can not be presented.This
Outside, the process nature of section gravitational inversion is that constantly adjustment model carries out forward fitting, and key is to utilize accurate elder generation
It tests information (such as deep reflection seismic profile) and establishes initial model, this is difficult to realize in 3-d inversion, it is necessary to rely on it
His method establishes the shell curtain density contrast model of variation, to improve the accurate of regional Moho depth and CRUSTAL THICKNESS inversion result
Property, it has important practical significance for infrastructure research, regional structure research etc..
Summary of the invention
The invention aims to solve the deficiencies in the prior art, a kind of CRUSTAL THICKNESS based on variable density is provided
Gravitational inversion method is established on the interface (i.e. Mohorovich discontinuity) of the earth's crust and earth mantle according to the density data that earth's crust different layers position changes
Under with change in location density contrast, using force density Interface Inversion technology inverting Moho depth, and then obtain the thickness of the earth's crust
Degree.
In order to achieve the above objectives, the present invention is implemented according to following technical scheme:
A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density, comprising the following steps:
S1,1.0 CRUSTAL THICKNESS of Crustal Model CRUST and density that spatial resolution covering the whole world is 1 ° × 1 ° are obtained
Data and density of upper mantle data；
S2, it is built according to the upper, middle and lower CRUST DENSITY and thickness and density of upper mantle data of Crustal Model CRUST 1.0
The vertical shell curtain density contrast changed with horizontal position；
S3, gravity anomaly caused by Moho relief underneath changes is obtained by bouguer gravity anomaly；
S4, caused gravity anomaly inversion Moho depth is changed by Moho relief underneath and obtains CRUSTAL THICKNESS.
Further, the specific steps of the S1 are as follows:
S11, spatial resolution is opened for 1 ° × 1 ° global Crustal Model CRUST1.0 database, in Download
CRUST1.0 data file crust1.0.tar.gz is downloaded in Section, includes using in file crust1.0.tar.gz
The code that formula translation is write, to extract CRUSTAL THICKNESS and density data and density of upper mantle data；
S12, getCN1xyz.f file in downloading file is run using Microsoft Visual Studio platform, from number
According to 179.5 ° of W~179.5 ° E covering the whole world are read in file, 8 density layers within the scope of 89.5 ° of S~89.5 ° N longitudes and latitudes
The thickness file of (water layer, ice sheet, upper sedimentary, middle sedimentary, lower sedimentary, upper crust, the middle earth's crust, lower crust) and this 8
The density data file of a density layer and outer mantle；
S13, according to research area specific longitude and latitude range, from the data file of S12 respectively extract research area within the scope of
Each density layer thickness and density data, and it is organized into each density layer thickness and each density layer density both of these documents respectively.
Further, the specific steps of the S2 are as follows:
S21, projective transformation is carried out to each density layer thickness and density data put in order in S13, geographical coordinate is projected
For plane rectangular coordinates；
S22, the shell curtain density contrast that each point at research area's different level position is calculated using following formula:
Wherein i indicates different level location point, Δ ρiFor the shell curtain density contrast at the i-th point, The density value of upper crust, the middle earth's crust, lower crust and outer mantle at respectively i-th point,The thickness of upper crust, the middle earth's crust, lower crust at respectively i-th point；
S23, to Δ ρiData gridding is carried out, the shell with the research same size of area's bouguer gravity anomaly data grids is obtained
Then curtain density contrast grid data is filtered shell curtain density difference data, obtain filtered shell curtain density difference data.
Further, the specific steps of the S3 are as follows:
S31, the deposit thickness number that research area is obtained according to shallow earthquake, drilling well or comprehensive geophysical inversion result
According to, and by U.S.National Oceanic and Atmospheric Administration NOAA acquisition sea area deposit thickness data and by CRUST1.0 data
Deposit thickness data in file obtain land-based area deposit thickness data；The sea area deposit thickness data resolution is 5 '
× 5 ', contain the sea area data in the whole world；
S32, projective transformation is carried out to deposit thickness data, geographical coordinate is projected as plane rectangular coordinates, it is right later
Deposit thickness data carry out gridding, and sizing grid and bouguer gravity anomaly data grids are in the same size；
S33, gravity anomaly caused by deposit thickness changes is calculated, terrain data and deposit thickness is utilized when calculating
Data obtain bottom boundary of sedimentary stratum depth data；The terrain data is obtained by U.S.National Oceanic and Atmospheric Administration NOAA
The terrain elevation data ETOPO1 in area is studied, which is 1 ' × 1 ', and is converted into plane rectangular coordinates,
Gridding is carried out to data, sizing grid and bouguer gravity anomaly data grids are in the same size；Frequency domain density circle is utilized later
Face gravity anomaly Forward Formula calculates gravity anomaly caused by deposit thickness changes:
Wherein G is universal gravitational constant；Indicate two-dimentional wave number field； It is heavy calculating
When gravity anomaly caused by lamination thickness change, the density contrast of sedimentary and the earth's crustIt is established according to CRUST1.0 data；It
Afterwards using the frequency spectrum of gravitational field caused by frequency domain density interface gravity anomaly forward modelling deposit thickness, then pass through Fourier
Inverse transformation obtains GRAVITY ANOMALIES；
S34, gravity anomaly caused by deposit thickness changes, the sedimentary that is eliminated shadow are subtracted from bouguer gravity anomaly
Gravity anomaly after sound is filtered the gravity anomaly, eliminates the gravity anomaly of local density's heterogeneous body in the earth's crust, it
Utilize afterwards the corresponding GRAVITY ANOMALIES of Mohorovich discontinuity known depth point obtained by deep reflection seismic profile data to above-mentioned filtering after
Gravity anomaly carry out background adjustment, obtain Moho relief underneath change caused by gravity anomaly.
Further, the specific steps of the S4 are as follows:
S41, research area's Mohorovich discontinuity is determined according to deep reflection seismic profile, submarine earthquake section OBS, CRUST1.0 data
Mean depth；
S42, caused gravity anomaly, filtered shell curtain density contrast and Mohorovich discontinuity are changed averagely deeply with Moho relief underneath
Degree is basic data, obtains Moho depth using frequency domain density-interface inversion method Parker-Oldenburg inverting；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect CRUSTAL THICKNESS
Variation characteristic, smooth hypsography is subtracted from Moho depth later, obtain research area's CRUSTAL THICKNESS.
In addition, technical solution as a further improvement of the present invention, the specific steps of the S4 can be with are as follows:
S41, research area's Mohorovich discontinuity is determined according to deep reflection seismic profile, submarine earthquake section OBS, CRUST1.0 data
Mean depth；
S42, it establishes after the shell curtain density contrast of cross directional variations, utilizes density-interface inversion inversion method Moho depth；
S421, its average shell curtain density contrast Δ ρ is calculated according to the shell curtain density contrast in entire research area0, calculate shell curtain density
The deviation δ ρ of differencei=Δ ρi-Δρ0；
S422, Moho depth data and Mohorovich discontinuity mean depth and shell curtain density according to CRUST1.0 database
The deviation of difference is drawn using gravity anomaly caused by formula (1) forward modelling shell curtain Density inhomogeneity, and from Moho relief underneath variation
It is eliminated in the gravity anomaly risen；
S423, to eliminate gravity anomaly after shell curtain Density inhomogeneity and average shell curtain density contrast Δ ρ0Based on number
According to obtaining Moho depth using frequency domain density-interface inversion method Parker-Oldenburg inverting；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect CRUSTAL THICKNESS
Variation characteristic, smooth hypsography is subtracted from Moho depth later, obtain research area's CRUSTAL THICKNESS.
Compared with prior art, the present invention establishes the shell curtain density contrast of cross directional variations, and is used for Moho depth and ground
Moho relief underneath and CRUSTAL THICKNESS feature can be accurately presented in thickness of the shell inverting, when survey region range is larger, or construction
More complicated region is not inconsistent using the shell curtain density contrast and practical geologic feature of constant, and only close along the two-dimentional shell curtain of section
Degree difference can not objectively reflect the CRUSTAL THICKNESS variation of whole region again, and technology provided by the invention is made using published data
Based on, practical shell curtain density model is established, it is technically easy to accomplish, and inversion result more meets practical geologic feature,
Accuracy is higher.
Detailed description of the invention
Fig. 1 is the flow chart of the CRUSTAL THICKNESS gravitational inversion method provided by the invention based on variable density.
Fig. 2 is the topographic map of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area.
Fig. 3 is the Bouguer map of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area.
Fig. 4 is the shell curtain density contrast isogram of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area.
Fig. 5 is the deposit thickness isogram of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area.
Fig. 6 is the sedimentary of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area and the density contrast isogram of the earth's crust.
Fig. 7 is gravity contour figure caused by embodiment Indian Ocean southeast mid-oceanic ridge and the deposit thickness of adjacent area.
Fig. 8 is the Mohorovich discontinuity gravity anomaly isogram of embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area.
Fig. 9 is embodiment Indian Ocean southeast mid-oceanic ridge and the Moho depth isogram that adjacent area inverting obtains.
Figure 10 is embodiment Indian Ocean southeast mid-oceanic ridge and the CRUSTAL THICKNESS isogram that adjacent area inverting obtains.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention more comprehensible, with reference to embodiments, to the present invention into
Row is further to be described in detail.Described herein the specific embodiments are only for explaining the present invention, is not used to limit invention.
As shown in Figure 1, a kind of CRUSTAL THICKNESS gravitational inversion method based on variable density provided by the invention, including following step
It is rapid:
S1, Crustal Model CRUST1.0 CRUSTAL THICKNESS and density number that spatial resolution covering the whole world is 1 ° × 1 ° are obtained
Accordingly and density of upper mantle data:
S11, open spatial resolution be 1 ° × 1 ° whole world Crustal Model CRUST1.0 database (network address:
Igppweb.ucsd.edu/~gabi/crust1.html), CRUST1.0 data file is downloaded in Download Section
It include the code write with formula translation in crust1.0.tar.gz, file crust1.0.tar.gz, to extract ground
Thickness of the shell and density data and density of upper mantle data；
S12, getCN1xyz.f file in downloading file is run using Microsoft Visual Studio platform, from number
According to 179.5 ° of W~179.5 ° E covering the whole world are read in file, 8 density layers within the scope of 89.5 ° of S~89.5 ° N longitudes and latitudes
The thickness file of (water layer, ice sheet, upper sedimentary, middle sedimentary, lower sedimentary, upper crust, the middle earth's crust, lower crust) and this 8
The density data file of a density layer and outer mantle；
S13, according to research area specific longitude and latitude range, from the data file of S12 respectively extract research area within the scope of
Each density layer thickness and density data, and it is organized into each density layer thickness and each density layer density both of these documents respectively.
S2, it is built according to the upper, middle and lower CRUST DENSITY and thickness and density of upper mantle data of Crustal Model CRUST 1.0
The vertical shell curtain density contrast changed with horizontal position:
S21, projective transformation is carried out to each density layer thickness and density data put in order in S13, geographical coordinate is projected
For plane rectangular coordinates；
S22, the shell curtain density contrast that each point at research area's different level position is calculated using following formula:
Wherein i indicates different level location point, Δ ρiFor the shell curtain density contrast at the i-th point, The density value of upper crust, the middle earth's crust, lower crust and outer mantle at respectively i-th point,The thickness of upper crust, the middle earth's crust, lower crust at respectively i-th point；
S23, to Δ ρiData gridding is carried out, the shell with the research same size of area's bouguer gravity anomaly data grids is obtained
Then curtain density contrast grid data is filtered shell curtain density difference data, obtain filtered shell curtain density difference data.
S3, gravity anomaly caused by Moho relief underneath changes is obtained by bouguer gravity anomaly:
S31, the deposit thickness number that research area is obtained according to shallow earthquake, drilling well or comprehensive geophysical inversion result
According to, if above data is insufficient, can by U.S.National Oceanic and Atmospheric Administration NOAA (network address:
Www.ngdc.noaa.gov/mgg/sedthick/ sea area deposit thickness data) are obtained and by CRUST1.0 data file
Deposit thickness data obtain land-based area deposit thickness data；The sea area deposit thickness data resolution is 5 ' × 5 ',
Contain the sea area data in the whole world；
S32, projective transformation is carried out to deposit thickness data, geographical coordinate is projected as plane rectangular coordinates, it is right later
Deposit thickness data carry out gridding, and sizing grid and bouguer gravity anomaly data grids are in the same size；
S33, gravity anomaly caused by deposit thickness changes is calculated, terrain data and deposit thickness is utilized when calculating
Data obtain bottom boundary of sedimentary stratum depth data；The terrain data is obtained by U.S.National Oceanic and Atmospheric Administration NOAA
Study the terrain elevation data ETOPO1 (network address: https: //www.ngdc.noaa.gov/mgg/global/ in area
Global.html), which is 1 ' × 1 ', and is converted into plane rectangular coordinates, carries out grid to data
Change, sizing grid and bouguer gravity anomaly data grids are in the same size；The gravity anomaly forward modeling of frequency domain density interface is utilized later
Formula calculates gravity anomaly caused by deposit thickness changes:
Wherein G is universal gravitational constant；Indicate two-dimentional wave number field； It is heavy calculating
When gravity anomaly caused by lamination thickness change, the density contrast of sedimentary and the earth's crustIt is established according to CRUST1.0 data；It
Afterwards using the frequency spectrum of gravitational field caused by frequency domain density interface gravity anomaly forward modelling deposit thickness, then pass through Fourier
Inverse transformation obtains GRAVITY ANOMALIES；
S34, gravity anomaly caused by deposit thickness changes, the sedimentary that is eliminated shadow are subtracted from bouguer gravity anomaly
Gravity anomaly after sound is filtered the gravity anomaly, eliminates the gravity anomaly of local density's heterogeneous body in the earth's crust, it
Utilize afterwards the corresponding GRAVITY ANOMALIES of Mohorovich discontinuity known depth point obtained by deep reflection seismic profile data to above-mentioned filtering after
Gravity anomaly carry out background adjustment, obtain Moho relief underneath change caused by gravity anomaly.
S4, caused gravity anomaly inversion Moho depth is changed by Moho relief underneath and obtains CRUSTAL THICKNESS:
S41, research area's Mohorovich discontinuity is determined according to deep reflection seismic profile, submarine earthquake section OBS, CRUST1.0 data
Mean depth；
S42, caused gravity anomaly, filtered shell curtain density contrast and Mohorovich discontinuity are changed averagely deeply with Moho relief underneath
Degree is basic data, obtains Moho depth using frequency domain density-interface inversion method Parker-Oldenburg inverting；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect CRUSTAL THICKNESS
Variation characteristic, smooth hypsography is subtracted from Moho depth later, obtain research area's CRUSTAL THICKNESS.
S42 in the step can also be substituted with following scheme, the shell curtain density contrasts of specifically vertical cross directional variations it
Afterwards, density-interface inversion inversion method Moho depth is utilized:
S421, its average shell curtain density contrast Δ ρ is calculated according to the shell curtain density contrast in entire research area0, calculate shell curtain density
The deviation δ ρ of differencei=Δ ρi-Δρ0；
S422, Moho depth data and Mohorovich discontinuity mean depth and shell curtain density according to CRUST1.0 database
The deviation of difference is caused using gravity anomaly caused by following formula forward modelling shell curtain Density inhomogeneity, and from Moho relief underneath variation
Gravity anomaly in eliminate:
Wherein G is universal gravitational constant；Indicate two-dimentional wave number field； It is heavy calculating
When gravity anomaly caused by lamination thickness change, the density contrast of sedimentary and the earth's crustIt is established according to CRUST1.0 data；It
Afterwards using the frequency spectrum of gravitational field caused by frequency domain density interface gravity anomaly forward modelling deposit thickness, then pass through Fourier
Inverse transformation obtains GRAVITY ANOMALIES；
S423, to eliminate gravity anomaly after shell curtain Density inhomogeneity and average shell curtain density contrast Δ ρ0Based on number
According to obtaining Moho depth using frequency domain density-interface inversion method Parker-Oldenburg inverting；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect CRUSTAL THICKNESS
Variation characteristic, smooth hypsography is subtracted from Moho depth later, can equally obtain research area's CRUSTAL THICKNESS.
In order to make those skilled in the art more fully understand the technical solution in this specification, below with the Indian Ocean southeast
It is right for mid-oceanic ridge and adjacent area Moho depth and CRUSTAL THICKNESS inverting (research range is 75 ° of E~150 ° E, 30 ° of S~65 ° S)
Technical solution of the present invention explains.
CRUST1.0 data file is downloaded by global Crustal Model CRUST1.0 database first
Crust1.0.tar.gz, and utilize getCN1xyz.f text in Microsoft Visual Studio platform operation downloading file
Part, with it is middle read each density layer covering the whole world (water layer, ice sheet, upper sedimentary, middle sedimentary, lower sedimentary, upper crust,
The middle earth's crust, lower crust) thickness and this 8 density layers and outer mantle density data.Indian Ocean southeast midocean is extracted respectively
It (75 ° of E~150 ° E, 30 ° of S~65 ° S) each density layer thickness and density data and is arranged in ridge and adjacent area research range, it
Projective transformation is carried out afterwards obtains each density layer thickness of Indian Ocean southeast mid-oceanic ridge and adjacent area and density under plane right-angle coordinate
Data.
Download Indian Ocean southeast mid-oceanic ridge and adjacent area landform and bouguer gravity anomaly data, this example mesorelief data and cloth
Lattice gravity anomaly data are both from global Gravity Models WGM2012 (network address: http://bgi.omp.obs-mip.fr/
data-products/Grids-and-models/wgm2012).Fig. 2 is embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area
Topographic map, wherein the overwhelming majority depth of water in sea area has been more than 3500m, the sea-floor relief at mid-oceanic ridge in apparent protuberance, be by
Caused by being gushed in oceanic crust expansion, earth mantle substance.Fig. 3 is the Bouguer map of Indian Ocean southeast mid-oceanic ridge and adjacent area, Bouguer
Gravity anomaly overall trend reflects the variation of Moho relief underneath, therefore in land area (such as Australia of map sheet northwest corner)
Gravity anomaly negative value is relatively low, and higher in ocean region (region of large area in the middle part of such as map sheet) gravity anomaly negative value, instead
It is larger to have reflected land area Moho depth, and ocean region Moho depth is relatively shallower.
Indian Ocean southeast mid-oceanic ridge and each horizontal position Dian Chu in adjacent area are then calculated according to weighted average density calculation formula
Shell curtain density contrast, formula is as follows:
Gridding is carried out to calculated result, is obtained close with the shell curtain of the research same size of area's bouguer gravity anomaly data grids
Poor grid data is spent, and to eliminate the influence of local density's variation, shell curtain density difference data is filtered, local density is eliminated
The shell curtain density contrast of inhomogeneities, finally obtained Indian Ocean southeast mid-oceanic ridge and adjacent area is as shown in Figure 4.Shell curtain density contrast is in land
Ground region is larger, is -0.5g/cm3Left and right, and the shell curtain density contrast in ocean region is then smaller, mostly less than -0.4g/cm3, this
It is since land area crustal structure is more complete, and caused by the general absence upper crust of ocean region.
Since embodiment Indian Ocean southeast mid-oceanic ridge and adjacent area lack the data such as shallow earthquake, drilling well, can not be obtained according to it
Must study area's deposit thickness data, thus by U.S.National Oceanic and Atmospheric Administration (NOAA) (network address:
Www.ngdc.noaa.gov/mgg/sedthick/ sea area deposit thickness data) are obtained, area's land-based area deposit thickness number is studied
It is supplemented according to by the deposit thickness data in CRUST1.0 data file.By data split and gridding, India is obtained
Foreign southeast mid-oceanic ridge and the deposit thickness of adjacent area are as shown in Figure 5.The maximum region of deposit thickness is located at Australian coastal waters
Region, maximum gauge are more than 7km, other area deposition thickness degree are smaller, and ocean region is generally less than 100m.
The density contrast of sedimentary and substrate is obtained as shown in fig. 6, using ground using the identical method of shell curtain density contrast is calculated
Shape and deposit thickness data obtain bottom boundary of sedimentary stratum depth data, respectively using landform and bottom boundary of sedimentary stratum as bound
Face calculates gravity anomaly caused by deposit thickness changes using frequency domain density interface gravity anomaly Forward Formula, and formula is such as
Under:
Wherein G is universal gravitational constant；Indicate two-dimentional wave number field； Utilize the public affairs
Formula calculates the frequency spectrum of gravitational field caused by deposit thickness, then is obtained caused by deposit thickness variation by Fourier inversion
Gravity anomaly is as shown in Figure 7.Since ocean area deposition layer is very thin, thus its caused gravity anomaly in most of region not
Foot -10mGal, the maximum place of deposit thickness, the caused reachable -160mGal of gravity anomaly maximum.
Gravity anomaly caused by deposit thickness changes is subtracted from bouguer gravity anomaly, the sedimentary that is eliminated influences it
Gravity anomaly afterwards is filtered the exception, eliminates the gravity anomaly of local density's heterogeneous body in the earth's crust, utilizes later
The corresponding GRAVITY ANOMALIES of part Mohorovich discontinuity known depth point carries out background to above-mentioned filtered gravity anomaly in CRUST1.0
It is as shown in Figure 8 to obtain gravity anomaly caused by Moho relief underneath changes for adjustment.This gravity anomaly and bouguer gravity anomaly entirety area
Domain is almost the same, only in regional area difference, is mainly manifested in the region near mid-oceanic ridge.
The mean depth of research area's Mohorovich discontinuity is determined according to CRUST1.0 data etc., weight caused by changing with Moho relief underneath
Power exception, shell curtain density contrast and Mohorovich discontinuity mean depth are basic data, (i.e. using frequency domain density-interface inversion method
Parker-Oldenburg method) to obtain Moho depth as shown in Figure 9 for inverting.According to inversion result, the Indian Ocean of embodiment
Southeast mid-oceanic ridge and adjacent area Moho depth are 10~38km, and the maximum region of depth is located at Australia, and Moho depth is general
All over being more than 30km, most areas Moho depth is close to 10km, individual pocket Moho depth phases in the Indian Ocean
It may be the reflection in micro- plot of big midocean close to 20km to larger.
Indian Ocean southeast mid-oceanic ridge and adjacent area terrain data are filtered, local relief variation, Zhi Houcong are eliminated
Smooth hypsography is subtracted in Moho depth, and it is as shown in Figure 10 to obtain research area's CRUSTAL THICKNESS.According to calculated result, implement
The Indian Ocean southeast mid-oceanic ridge and adjacent area CRUSTAL THICKNESS of example are 6~38km, and Australian CRUSTAL THICKNESS has generally been more than 32km,
In the Indian Ocean most areas CRUSTAL THICKNESS be 7km, individual fritter CRUSTAL THICKNESSs are more than 18km, may be big midocean micro-ly
The reflection of block.In addition, CRUSTAL THICKNESS at mid-ocean rise is 8~9km, slightly larger than the Crust thickness of its side, thus it is speculated that its be by
Mid-oceanic ridge slow expansion, earth mantle material buildup and caused by.
Above-described embodiment, which is further demonstrated, establishes the earth's crust and earth mantle according to the density data of earth's crust different layers position variation
Interface (i.e. Mohorovich discontinuity) is deep using force density Interface Inversion technology inverting Mohorovich discontinuity up and down with the density contrast of change in location
Degree, and then obtain the thickness of the earth's crust.As it can be seen that technology provided by the present invention has simple, easy to operate, and accuracy is higher
Feature.Compared with the existing inversion method using Chang Midu, inversion result accuracy of the present invention is higher, more meets the ground of the earth's crust
Matter feature can provide key technology support for infrastructure research and regional study.
The limitation that technical solution of the present invention is not limited to the above specific embodiments, it is all to do according to the technique and scheme of the present invention
Technology deformation out, falls within the scope of protection of the present invention.
Claims (7)
1. a kind of CRUSTAL THICKNESS gravitational inversion method based on variable density, which comprises the following steps:
S1, obtain Crustal Model CRUST1.0 CRUSTAL THICKNESS and the density data that spatial resolution covering the whole world is 1 ° × 1 ° with
And density of upper mantle data；
S2, it is established according to the upper, middle and lower CRUST DENSITY and thickness and density of upper mantle data of Crustal Model CRUST1.0 with water
The shell curtain density contrast of flat change in location；
S3, gravity anomaly caused by Moho relief underneath changes is obtained by bouguer gravity anomaly；
S4, caused gravity anomaly inversion Moho depth is changed by Moho relief underneath and obtains CRUSTAL THICKNESS.
2. the CRUSTAL THICKNESS gravitational inversion method according to claim 1 based on variable density, which is characterized in that the S1's
Specific steps are as follows:
S11, spatial resolution is opened for 1 ° × 1 ° global Crustal Model CRUST1.0 database, in Download Section
CRUST1.0 data file crust1.0.tar.gz is downloaded, includes using formula translation in file crust1.0.tar.gz
The code write, to extract CRUSTAL THICKNESS and density data and density of upper mantle data；
S12, getCN1xyz.f file in downloading file is run using Microsoft Visual Studio platform, from data text
The thickness of 179.5 ° of W~179.5 ° E covering the whole world, 8 density layers within the scope of 89.5 ° of S~89.5 ° N longitudes and latitudes is read in part
Spend file and the density data file of this 8 density layers and outer mantle；
S13, according to research area specific longitude and latitude range, from the data file of S12 respectively extract research area's range Nei Gemi
Thickness degree and density data are spent, and is organized into each density layer thickness and each density layer density both of these documents respectively.
3. the CRUSTAL THICKNESS gravitational inversion method according to claim 2 based on variable density, which is characterized in that the S2's
Specific steps are as follows:
S21, projective transformation is carried out to each density layer thickness and density data put in order in S13, geographical coordinate is projected as putting down
Face rectangular co-ordinate；
S22, the shell curtain density contrast that each point at research area's different level position is calculated using following formula:
Wherein i indicates different level location point, Δ ρiFor the shell curtain density contrast at the i-th point, Point
Not Wei at the i-th point upper crust, the middle earth's crust, lower crust and outer mantle density value,At respectively i-th point
Upper crust, the middle earth's crust, lower crust thickness；
S23, to Δ ρiData gridding is carried out, is obtained close with the shell curtain of the research same size of area's bouguer gravity anomaly data grids
Poor grid data is spent, then shell curtain density difference data is filtered, obtains filtered shell curtain density difference data.
4. the CRUSTAL THICKNESS gravitational inversion method according to claim 3 based on variable density, which is characterized in that the S3's
Specific steps are as follows:
S31, the deposit thickness data that research area is obtained according to shallow earthquake, drilling well or comprehensive geophysical inversion result, with
And by U.S.National Oceanic and Atmospheric Administration NOAA acquisition sea area deposit thickness data and by CRUST1.0 data file
In deposit thickness data obtain land-based area deposit thickness data；The sea area deposit thickness data resolution be 5 ' ×
5 ', contain the sea area data in the whole world；
S32, projective transformation is carried out to deposit thickness data, geographical coordinate is projected as plane rectangular coordinates, later to deposition
Layer thickness data carries out gridding, and sizing grid and bouguer gravity anomaly data grids are in the same size；
S33, gravity anomaly caused by deposit thickness changes is calculated, terrain data and deposit thickness data is utilized when calculating
Obtain bottom boundary of sedimentary stratum depth data；The terrain data obtains research by U.S.National Oceanic and Atmospheric Administration NOAA
The terrain elevation data ETOPO1 in area, which is 1 ' × 1 ', and is converted into plane rectangular coordinates, logarithm
According to gridding is carried out, sizing grid and bouguer gravity anomaly data grids are in the same size；Frequency domain density interface weight is utilized later
Power exception Forward Formula calculates gravity anomaly caused by deposit thickness changes:
Wherein G is universal gravitational constant；Indicate two-dimentional wave number field； Calculating sedimentary
Caused by thickness change when gravity anomaly, the density contrast of sedimentary and the earth's crustIt is established according to CRUST1.0 data；It is sharp later
The frequency spectrum of the gravitational field caused by frequency domain density interface gravity anomaly forward modelling deposit thickness, then pass through Fourier's contravariant
Get GRAVITY ANOMALIES in return；
S34, gravity anomaly caused by deposit thickness changes is subtracted from bouguer gravity anomaly, the sedimentary that is eliminated influences it
Gravity anomaly afterwards is filtered the gravity anomaly, eliminates the gravity anomaly of local density's heterogeneous body in the earth's crust, Zhi Houli
With the corresponding GRAVITY ANOMALIES of Mohorovich discontinuity known depth point obtained by deep reflection seismic profile data to above-mentioned filtered heavy
Power carries out background adjustment extremely, obtains gravity anomaly caused by Moho relief underneath changes.
5. the CRUSTAL THICKNESS gravitational inversion method according to claim 4 based on variable density, which is characterized in that the S4's
Specific steps are as follows:
S41, being averaged for research area's Mohorovich discontinuity is determined according to deep reflection seismic profile, submarine earthquake section OBS, CRUST1.0 data
Depth；
S42, it is with gravity anomaly, filtered shell curtain density contrast and Mohorovich discontinuity mean depth caused by Moho relief underneath variation
Basic data obtains Moho depth using frequency domain density-interface inversion method Parker-Oldenburg inverting；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect the change of CRUSTAL THICKNESS
Change feature, subtract smooth hypsography from Moho depth later, obtains research area's CRUSTAL THICKNESS.
6. the CRUSTAL THICKNESS gravitational inversion method according to claim 1 based on variable density according to claim 4,
It is characterized in that, the specific steps of the S4 are as follows:
S41, being averaged for research area's Mohorovich discontinuity is determined according to deep reflection seismic profile, submarine earthquake section OBS, CRUST1.0 data
Depth；
S42, it establishes after the shell curtain density contrast of cross directional variations, utilizes density-interface inversion inversion method Moho depth；
S421, its average shell curtain density contrast Δ ρ is calculated according to the shell curtain density contrast in entire research area0, calculate the inclined of shell curtain density contrast
Poor δ ρi=Δ ρi-Δρ0；
S422, according to the Moho depth data and Mohorovich discontinuity mean depth of CRUST1.0 database and shell curtain density contrast
Deviation, using gravity anomaly caused by formula (1) forward modelling shell curtain Density inhomogeneity, and from Moho relief underneath variation caused by
It is eliminated in gravity anomaly；
S423, to eliminate gravity anomaly after shell curtain Density inhomogeneity and average shell curtain density contrast Δ ρ0For basic data, utilize
Frequency domain density-interface inversion method Parker-Oldenburg inverting obtains Moho depth；
S43, research area's terrain data is filtered, eliminates local relief variation objectively to reflect the change of CRUSTAL THICKNESS
Change feature, subtract smooth hypsography from Moho depth later, obtains research area's CRUSTAL THICKNESS.
7. the CRUSTAL THICKNESS gravitational inversion method according to claim 1 based on variable density according to claim 2,
It is characterized by: 8 density layers are specially water layer, ice sheet, upper sedimentary, middle sedimentary, lower sedimentary, Shangdi in the S12
Shell, the middle earth's crust, lower crust.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910502882.4A CN110244352B (en) | 2019-06-11 | 2019-06-11 | Variable density-based crust thickness gravity inversion method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910502882.4A CN110244352B (en) | 2019-06-11 | 2019-06-11 | Variable density-based crust thickness gravity inversion method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110244352A true CN110244352A (en) | 2019-09-17 |
| CN110244352B CN110244352B (en) | 2020-12-25 |
Family
ID=67886681
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910502882.4A Expired - Fee Related CN110244352B (en) | 2019-06-11 | 2019-06-11 | Variable density-based crust thickness gravity inversion method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110244352B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111045098A (en) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | Method for picking up underground deep structure information |
| CN111337993A (en) * | 2020-03-30 | 2020-06-26 | 中国地质科学院 | Variable density and variable depth constraint-based gravity density interface inversion method |
| CN111999778A (en) * | 2020-07-13 | 2020-11-27 | 国家海洋信息中心 | Antarctic continental Mohuo surface depth inversion method based on satellite gravity gradient data |
| CN112666607A (en) * | 2019-10-16 | 2021-04-16 | 中国石油天然气集团有限公司 | Method and device for gravity inversion of loess layer thickness distribution |
| CN113985490A (en) * | 2021-09-22 | 2022-01-28 | 中国人民解放军战略支援部队信息工程大学 | Method and device for performing surface gravity simulation by using terrain and crust density data |
| CN115373024A (en) * | 2022-08-09 | 2022-11-22 | 中国科学院南海海洋研究所 | Method and device for inverting passive land edge crustal structure based on stratum recording settlement |
| WO2023280122A1 (en) * | 2021-07-05 | 2023-01-12 | 中国科学院地质与地球物理研究所 | Density determination method and apparatus, and electronic device |
| CN117555025A (en) * | 2024-01-11 | 2024-02-13 | 应急管理部国家自然灾害防治研究院 | Multi-layer crust structure inversion method based on gravity data |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013087881A1 (en) * | 2011-12-15 | 2013-06-20 | 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 |
| CN104835201A (en) * | 2015-05-08 | 2015-08-12 | 华东师范大学 | Method for simulating and displaying global crustal structure on digital earth software platform |
| CN106886047A (en) * | 2017-02-28 | 2017-06-23 | 中国地质大学（北京） | A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio |
| WO2017131848A1 (en) * | 2016-01-26 | 2017-08-03 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
| CN108919338A (en) * | 2018-05-28 | 2018-11-30 | 中国地震局地震预测研究所 | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake |
-
2019
- 2019-06-11 CN CN201910502882.4A patent/CN110244352B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013087881A1 (en) * | 2011-12-15 | 2013-06-20 | 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 |
| CN104835201A (en) * | 2015-05-08 | 2015-08-12 | 华东师范大学 | Method for simulating and displaying global crustal structure on digital earth software platform |
| WO2017131848A1 (en) * | 2016-01-26 | 2017-08-03 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
| CN106886047A (en) * | 2017-02-28 | 2017-06-23 | 中国地质大学（北京） | A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio |
| CN108919338A (en) * | 2018-05-28 | 2018-11-30 | 中国地震局地震预测研究所 | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake |
Non-Patent Citations (6)
| Title |
|---|
| LASKE G等: ""Update on CRUST1.0-A1-degree Global Model of Earth’s Crust", 《GEOPHYS. RES. ABSTRACTS》 * |
| MOHAMMAD BAGHERBANDI: ""MohoIso: A MATLAB program to determine crustal thickness by an isostatic and a global gravitational model"", 《COMPUTERS & GEOSCIENCES》 * |
| PAVLIS N K等: ""The development and evaluation of the Earth Gravitational Model 2008（EGM2008）"", 《J. GEOPHYS. RES》 * |
| 冯旭亮 等: ""渭河盆地基底三维密度重力反演"", 《石油地球物理勘探》 * |
| 冯旭亮: ""空间域密度界面反演方法及其进展"", 《地球科学进展》 * |
| 陈国安 等: ""长江中下游成矿带及邻区Moho深度与成矿背景探讨"", 《地质学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112666607A (en) * | 2019-10-16 | 2021-04-16 | 中国石油天然气集团有限公司 | Method and device for gravity inversion of loess layer thickness distribution |
| CN111045098A (en) * | 2019-12-27 | 2020-04-21 | 核工业北京地质研究院 | Method for picking up underground deep structure information |
| CN111337993A (en) * | 2020-03-30 | 2020-06-26 | 中国地质科学院 | Variable density and variable depth constraint-based gravity density interface inversion method |
| CN111999778A (en) * | 2020-07-13 | 2020-11-27 | 国家海洋信息中心 | Antarctic continental Mohuo surface depth inversion method based on satellite gravity gradient data |
| WO2023280122A1 (en) * | 2021-07-05 | 2023-01-12 | 中国科学院地质与地球物理研究所 | Density determination method and apparatus, and electronic device |
| CN113985490A (en) * | 2021-09-22 | 2022-01-28 | 中国人民解放军战略支援部队信息工程大学 | Method and device for performing surface gravity simulation by using terrain and crust density data |
| CN115373024A (en) * | 2022-08-09 | 2022-11-22 | 中国科学院南海海洋研究所 | Method and device for inverting passive land edge crustal structure based on stratum recording settlement |
| CN115373024B (en) * | 2022-08-09 | 2023-04-18 | 中国科学院南海海洋研究所 | Method and device for inverting passive land edge crustal structure based on stratum recording settlement |
| WO2023098441A1 (en) * | 2022-08-09 | 2023-06-08 | 中国科学院南海海洋研究所 | Method and apparatus for inverting crustal structure of passive continental margin on basis of subsidence in stratigraphic record |
| CN117555025A (en) * | 2024-01-11 | 2024-02-13 | 应急管理部国家自然灾害防治研究院 | Multi-layer crust structure inversion method based on gravity data |
| CN117555025B (en) * | 2024-01-11 | 2024-04-02 | 应急管理部国家自然灾害防治研究院 | Multi-layer crust structure inversion method based on gravity data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110244352B (en) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110244352A (en) | A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density | |
| Hinze | Bouguer reduction density, why 2.67? | |
| Samrock et al. | 3-D analysis and interpretation of magnetotelluric data from the Aluto-Langano geothermal field, Ethiopia | |
| Tian et al. | Mantle transition zone topography and structure beneath the central Tien Shan orogenic belt | |
| CN104765065B (en) | Method and device for determining sedimentary stratum ancient landform | |
| Kopp et al. | Deep structure of the central Lesser Antilles Island Arc: relevance for the formation of continental crust | |
| Sircar et al. | Geothermal exploration in Gujarat: case study from Dholera | |
| Zaher et al. | Exploration and assessment of the geothermal resources in the Hammam Faraun hot spring, Sinai Peninsula, Egypt | |
| Gautam et al. | 2D Geo-electrical imaging for shallow depth investigation in Doon Valley Sub-Himalaya, Uttarakhand, India | |
| Bailey et al. | VS and density structure beneath the Colorado Plateau constrained by gravity anomalies and joint inversions of receiver function and phase velocity data | |
| Pfeffer et al. | Evaluating surface and subsurface water storage variations at small time and space scales from relative gravity measurements in semiarid Niger | |
| CN105425315A (en) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range | |
| Tiede et al. | Modeling the density at Merapi volcano area, Indonesia, via the inverse gravimetric problem | |
| Abdelzaher et al. | A coastal aquifer study using magnetotelluric and gravity methods in Abo Zenema, Egypt | |
| Karpiah et al. | Crustal structure and basement-cover relationship in the Dangerous Grounds, offshore North-West Borneo, from 3D joint CSEM and MT imaging | |
| Lodolo et al. | Gravity map of the Isla Grande de Tierra del Fuego, and morphology of Lago Fagnano | |
| Li et al. | Sedimentary basin analysis using airborne gravity data: a case study from the Bohai Bay Basin, China | |
| CN114518605A (en) | Low-altitude, shallow water and deep water integrated geological measurement method based on electromagnetic method | |
| Cruz-Falcón et al. | Depth estimation to crystalline basement in the valley of La Paz, Baja California Sur, Mexico | |
| Joel et al. | Identification of geothermal reservoirs in south cameroon from a combined landsat 8 and amt data | |
| Zhao et al. | The implication of the geophysical data and drilling records for the formation and paleoenvironment of the Hobq Desert, China | |
| Corbo et al. | Structure of the Guaraní aquifer in the surroundings of the Uruguay river from magnetotelluric soundings | |
| Pospil | Crustal and Lithospheric Structure of the Carpathian–Pannonian Region: A Geophysical PerspectiveRegional Geophysical Data on the Carpathian–Pannonian Lithosphere | |
| Gisel Peri et al. | Shallow geophysical evaluation of the transition zone between the Guaraní and Yrendá-Toba-Tarijeño aquifer systems (Argentine Gran Chaco) | |
| Rao et al. | Electrical conductance map for the Kachchh Rift Basin: constraint on tectonic evolution and seismotectonic implications |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201225 Termination date: 20210611 |