CN108919338A - Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake - Google Patents
Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake Download PDFInfo
- Publication number
- CN108919338A CN108919338A CN201810520390.3A CN201810520390A CN108919338A CN 108919338 A CN108919338 A CN 108919338A CN 201810520390 A CN201810520390 A CN 201810520390A CN 108919338 A CN108919338 A CN 108919338A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- earth
- lithosphere
- tectonic stress
- carrying
- 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
- 230000005484 gravity Effects 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000010276 construction Methods 0.000 claims abstract description 11
- 241001191378 Moho Species 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 23
- 230000001960 triggered effect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 23
- 238000012937 correction Methods 0.000 description 6
- 238000011144 upstream manufacturing Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 241000209094 Oryza Species 0.000 description 3
- 235000007164 Oryza sativa Nutrition 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005429 filling process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 235000009566 rice Nutrition 0.000 description 3
- 240000006108 Allium ampeloprasum Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000005945 translocation Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000001976 improved effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000002453 shampoo Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012876 topography 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/01—Measuring or predicting earthquakes
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method of reservoir-induced earthquake is predicted based on earth's surface gravity and GNSS observation data, that is, calculates the vertical tectonic stress of proposed reservoir location/section lithosphere carrying；If the vertical tectonic stress of lithosphere carrying, far from the earth's core (direction is upward), which is suitble to build large reservoir, reservoir macroseism will not be triggered；If the vertical tectonic stress of lithosphere carrying is directed toward the earth's core (direction is downward), which is not suitable for construction large reservoir, there is a possibility that triggering reservoir macroseism.The present invention builds the relationship between adjoint earth's surface load by the vertical tectonic stress field and large reservoir of analysis lithosphere carrying, and the new method of prediction reservoir-induced earthquake is provided from the angle of lithospheric isostasy.The present invention can be in the reservoir proposed stage, the seismic risk that reservoir filling may induce is predicted earlier according to earth's surface gravity and GNSS observation data, it not only contributes to mitigate reservoir-induced earthquake bring disaster, can also provide theoretical foundation for large-scale water conservancy construction addressing.
Description
Technical field
The present invention relates to a kind of methods for predicting reservoir-induced earthquake, specifically, the present invention relates to one kind to be based on earth's surface
The method of gravity and GNSS observation data prediction reservoir-induced earthquake.
Background technique
As global reservoir is built with annual 11 or so speed, reservoir-induced earthquake disaster and its risk are increasingly
Cause the extensive concern of people.Reservoir-induced earthquake is a kind of special seismic activity, mainly due to reservoir filling/sluiced
The reservoir area of Cheng Yinfa and neighbouring seismic activity.
According to incompletely statistics, the whole world has occurred reservoir-induced earthquake 134, wherein 6 grades or more earthquake 4, in the world most
Big reservoir-induced earthquake is the 6.5 grades of earthquakes of India Ge Yina reservoir, and the maximum reservoir-induced earthquake in China is March 19 in 1962
6.1 grades of earthquakes of Xinfengjiang Reservoir that day occurs.The generation of these reservoir-induced earthquakes is brought safely to the people's lives and property
Immeasurable loss.
As the country of water resource relative abundance, so far, research of the China in terms of reservoir-induced earthquake mainly collects
In around it is following eight aspect expansion：1. reservoir height of dam, library type, storage capacity, storage (let out) correlation of water process with seismic activity；②
The correlation of stratum media, geological structure, topography and geomorphology and seismic activity；3. the correlation of fault activities and seismic activity；④
The correlation of crustal deformation and seismic activity；5. the correlation of hydrogeological environment and seismic activity；6. crustal stress, pore pressure
With the correlation of seismic activity；7. reservoir-induced earthquake hair shake mark and Forming Mechanism；8. reservoir-induced earthquake hazard assessment is monitored and is answered
Anxious prediction scheme.
The studies above is mostly the correlation inquired into from different perspectives after reservoir is built up with seismic activity, obtains difference
The cognition about reservoir-induced earthquake of degree has important scientific value.But if with can proposing a kind of effective reservoir
Earthquake prediction method can predict increase and decrease and the Reservoir induced seismicity of seismicity caused by reservoir filling before reservoir construction
The presence or absence of earthquake is sent out, will be worth with more scientific meaning and mitigation.
Summary of the invention
To reduce the earthquake disaster induced due to reservoir construction, the object of the present invention is to provide one kind based on earth's surface gravity with
The method that GNSS observes data prediction reservoir-induced earthquake.This method can predict related reservoir before reservoir is formally built
The seismic risk that water storage process may induce, and providing is suitable for the location/section suggestion for building reservoir, to mitigate earthquake
Disaster and casualties resulting from and property loss.
To achieve the above object, the present invention uses following technical scheme：One kind observing data based on earth's surface gravity and GNSS
The method for predicting reservoir-induced earthquake, it is characterised in that：Calculate the vertical construction in proposed reservoir location or the carrying of section lithosphere
Stress；If the vertical tectonic stress of lithosphere carrying, far from the earth's core, which is suitble to build large reservoir, it is not easy to touch
Shampoo library macroseism；If the vertical tectonic stress of lithosphere carrying is directed toward the earth's core, which is not suitable for building large-scale water
There is a possibility that triggering reservoir macroseism in library.
The vertical tectonic stress F of the lithosphere carryingIt is verticalFor：
Wherein, h1For isostasy face depth, h2For Crust Moho plane depth, Δ ρ is the density of earth mantle and crust materials
Difference, different regions Δ ρ value is different, and g is acceleration of gravity；
Isostasy face depth h is calculated according to GNSS observed result1, calculation formula is as follows：
Wherein, h*It is the surface observation station height above sea level that GNSS is observed, ρtIt is the density of the above substance in sea level, Δ ρ is
The density contrast of earth mantle and crust materials, H are that this area's b.s.l. is averaged CRUSTAL THICKNESS；
Crust Moho plane depth h is calculated according to earth's surface gravity and GNSS observation data2, calculation formula is as follows：
Wherein, G is universal gravitational constant, and Δ ρ is the density contrast of earth mantle and crust materials；h0It is the Moho with reference to observation point
Face depth, Δ g0For the bouguer gravity anomaly with reference to observation point；Δg*For the bouguer gravity anomaly at surface observation station.
It is positive when observing the vertical tectonic stress that the lithosphere that data determine carries based on earth's surface gravity and GNSS, shows this
For the vertical tectonic stress that location or section lithosphere are born far from the earth's core, direction is upward, opposite with earth's surface gravity direction；The rock
The vertical tectonic stress of stone circle carrying is negative, and the vertical tectonic stress for showing that the location or section lithosphere are born is directed toward the earth's core,
Direction is downward, identical as earth's surface gravity direction.
The present invention builds adjoint earth's surface by the vertical tectonic stress field and large reservoir of analysis lithosphere carrying and loads
Between relationship, give from the angle of lithospheric isostasy based on earth's surface gravity and GNSS observation data, prediction Reservoir-induced
The method of shake.Water storage process can be unfolded by earth's surface gravity and GNSS joint observation before reservoir is formally built in the present invention
A possibility that inducing macroseism is assessed, and predicts the seismic risk that related reservoir filling may induce earlier, while being capable of root
It is suitable for building location/section of reservoir, therefore the present invention not only contributes to mitigate reservoir-induced earthquake bring calamity according to needing to provide
Casualties and property loss harmful and resulting from, the decision-making work that can also build addressing for large-scale water conservancy (such as reservoir) mention
For the theoretical foundation and higher reference value of science.
Detailed description of the invention
Fig. 1 is the vertical tectonic stress variation schematic diagram of lithosphere carrying before and after reservoir filling；
Fig. 2 is the vertical tectonic stress computation model based on law of buoyancy；
Fig. 3 is Yarlung Zangbo Grand Canyon, Tibet, China surrounding area geologic structure diagram and earth's surface gravity and GNSS observational network；
Fig. 4 is Yarlung Zangbo Grand Canyon, Tibet, China surrounding area Lithospheric Structure stress field；
Fig. 5 is Yarlungzangbo Grand Canyon upstream and downstream river valley elevation and the vertical tectonic stress field pattern of lithosphere.
Specific embodiment
Structure and feature of the invention are described in detail with reference to the accompanying drawings and examples.It should be noted that can
To make various modifications to disclosed embodiments, therefore, embodiment disclosed in specification be should not be considered as to the present invention
Limitation, and only as the example of embodiment, the purpose is to keep feature of the invention apparent.
Reservoir-induced earthquake mainly due to reservoir filling or sluicing during, water body to surrounding rock circle generation act on,
To the reservoir area caused and neighbouring seismic activity.Reservoir-induced earthquake is generally divided into four kinds of situations:1. there is no historical earthquake before water storage
It records, occurs apparent seismic activity after water storage；2. the magnitude and frequency that occur after water storage are recorded higher than historical earthquake；③
The earthquake magnitude of earthquake is lower than the earthquake magnitude before water storage after water storage；4. Seismicity frequency and earthquake magnitude do not change significantly after water storage.Its
Middle first three phenomenon may be corresponded with three kinds of Main Morphologies of the vertical tectonic stress field of lithosphere, specifically, situation is 1.
It is likely to occur in the lesser area of vertical tectonic stress that lithosphere is substantially balanced, and lithosphere carries.Because the load of library water is made
With destroying lithospheric equilibrium state, therefore there is apparent seismic activity；2. situation is likely to occur in hanging down for lithosphere carrying
It is directed toward the downward area in the earth's core i.e. direction to tectonic stress, 3. situation is likely to occur in the vertical tectonic stress of lithosphere far from the earth's core
That is the upward area in direction.Situation is 4. planted as the, it may be possible to which, since reservoir is small, water body is self-possessed on lithospheric influence
It is too low caused.In short, the vertical tectonic stress field distributional pattern of lithosphere can explain above-mentioned four kinds of reservoir-induced earthquake well
Phenomenon.
Fig. 1 is the vertical tectonic stress variation schematic diagram of lithosphere carrying before and after reservoir filling.As shown, the present invention will
Reservoir construction is considered as the load caused by mankind's activity to earth's surface, and the water storage process of reservoir is considered as water body to lithosphere surface
A kind of stress loading effect for being directed toward the earth's core (i.e. direction is downward)；The earth's surface stress loading formed by comparing reservoir filling process
With the relationship of the vertical tectonic stress field of proposed Reservoir Area lithosphere carrying, determine that reservoir filling process answers Lithospheric Structure
The field of force adds unloading effect, to predict whether reservoir filling process may induce earthquake, and whether provides this location/section
It is suitable for the judgement for building reservoir.
As shown in Figure 1, if the vertical tectonic stress that carries of the lithosphere of proposed reservoir section far from the earth's core, that is, direction to
On, after reservoir filling, the self weight of extensive water storage will offset the vertical tectonic stress of part lithosphere carrying, generally form
" unloading " effect, so that the lithosphere of the section is more balanced, so that macroseism will not be triggered., whereas if proposed reservoir section
The vertical tectonic stress of lithosphere carrying be directed toward that the earth's core i.e. direction is downward or the lithosphere of proposed reservoir section is in basic
Equilibrium state, after reservoir filling, the self weight (i.e. the water storage process of reservoir) of extensive water storage hangs down formation to lithosphere carrying
It is acted on to " load " of tectonic stress, keeps lithosphere more unbalanced, to increase regional seismic activity, in some instances it may even be possible to trigger
Macroseism.
Therefore, from the point of view of lithospheric isostasy, the method that the present invention predicts reservoir-induced earthquake is：Based on earth's surface gravity
Data are observed with GNSS, calculate the vertical tectonic stress in proposed reservoir location or the carrying of section lithosphere, if lithosphere carries
Vertical tectonic stress it is upward far from the earth's core, that is, direction, the location or section are suitble to build large reservoir, are not easy to trigger reservoir strong
Shake；If to be directed toward the earth's core i.e. direction downward for the vertical tectonic stress of lithosphere carrying, it is large-scale that the location or section are not suitable for construction
Reservoir has a possibility that triggering reservoir macroseism.Specifically, the reservoir filling in the location will make lithosphere more unbalanced, from
And enhance regional seismic activity, in some instances it may even be possible to trigger macroseism.
The earth is the sphere with Circle Structure, can substantially be divided into the earth's crust, earth mantle and three layers of the earth's core.From very long ground
From the point of view of the matter age, the relatively light and hard earth's crust is swum on heavy and soft earth mantle, as solid ice swims in liquid water
On.As shown in Fig. 2, the earth's crust has elasticity for the practical earth, the tectonic stress of all directions can be carried, including
Vertical tectonic stress.Buoyancy caused by density layer between some areas isostatic surface and Mohorovich discontinuity is generally held by the earth's crust
It carries and balances.Theoretically, buoyancy caused by the density layer between isostatic surface and Mohorovich discontinuity is carried vertical outer with the earth's crust
Power is equal in magnitude contrary.Therefore, it can be quantified by the buoyancy that the density layer between isostasy face and Mohorovich discontinuity generates
The vertical tectonic stress of the earth's crust (i.e. lithosphere) carrying is calculated, calculation formula is：
Wherein, h1For isostasy face depth, h2For Crust Moho plane depth, g is acceleration of gravity；Δ ρ is earth mantle and ground
The density contrast of shell matter, different regions Δ ρ value is different, and the density of earth mantle generally can be taken as 3.3g/cm3, the density of the earth's crust is not
It is different with region, it is general to be obtained according to Bouguer gravity data inversion.
According to GNSS observed result calculably shell isostatic surface depth h1, calculation formula is：
Wherein, h1CRUSTAL THICKNESS when being at equilibrium, that is, the isostasy face to be calculated depth；h*It is that GNSS is seen
The surface observation point height above sea level measured, H are that this area's b.s.l. is averaged CRUSTAL THICKNESS, ρtIt is the above substance in sea level
Density, Δ ρ are the density contrasts of earth mantle and crust materials.
Data calculably shell Moho depth h is observed according to earth's surface gravity and GNSS2, calculation formula is：
Wherein, G is universal gravitational constant, and Δ ρ is the density contrast of earth mantle and crust materials；h0For the Moho with reference to observation point
Face depth, Δ g0For the bouguer gravity anomaly with reference to observation point；Δg*For the bouguer gravity anomaly at surface observation station.Bouguer gravity
It is abnormal to be obtained after adjustment, topographic correction, middle layer correction and altitude correction by earth's surface gravity and GNSS observation data, it calculates
Method is the common sense of those skilled in the art, and details are not described herein.
Isostasy face depth h is calculated according to formula (2) and (3)1With Moho depth h2, calculated further according to formula (1)
The vertical tectonic stress of the location/section lithosphere carrying shows this if the vertical tectonic stress of lithosphere carrying is positive out
The vertical tectonic stress of location or the carrying of section lithosphere is far from the earth's core, and direction is upward, with earth's surface gravity direction on the contrary, the location
Or section is suitble to build large reservoir, is not easy to trigger reservoir macroseism；If the vertical tectonic stress of lithosphere carrying is negative, show
The vertical tectonic stress of the location or the carrying of section lithosphere is directed toward the earth's core, and direction is downward, the ground identical as earth's surface gravity direction
Section or section are not suitable for construction large reservoir, and in the section, construction reservoir will enhance regional seismic activity, in some instances it may even be possible to trigger
Reservoir macroseism.
Below by specific embodiment, invention is further described in detail.
Since 20th century, 6 grades or more earthquakes 28 times occur altogether for Yarlungzangbo Grand Canyon surrounding area, including nineteen forty-seven
The earthquake of 7.7 grades of Lang County and the earthquake of 8.6 grades of nineteen fifty Chayu show that this area has the tectonic setting that violent earthquake occurs.2017
Year 6.9 grades of Milin earthquake occurs, and display this area's seismic activity has the tendency that active.
Under the premise of fully considering survey region transportation condition, as shown in figure 3, we are in Yarlung Zangbo Grand Canyon, Tibet, China week
Region In Yanbian establishes the high density gravimetric network being made of 107 observation stations, and complete a phase earth's surface gravity and
GNSS joint observation.The area Fromlingzhi, tibet Ba Yi, Milin County, Bowo County and Motuo County are covered in specific observation, to improve Grand Canyon week
Region In Yanbian gravity anomaly field accuracy and spatial resolution, and then obtain the vertical tectonic stress of the lithosphere carrying of this area
?.In Fig. 3, NB and GP respectively represent southern character used in proper names and in rendering some foreign names bar watt and the top Galle Bai Lei position, black dot and black thin item
Respectively newly-built gravimetric observation station and translocation route, the small box of black are " land state network " gravimetric observation station.We create gravity
Above three " land state network " gravimetric observation station is crossed by Netcom, is connected with national gravimetric network, to obtain absolute gravity observation base
The control at quasi- (such as Lhasa absolute gravity station).
Earth's surface gravity and the GNSS observation of each observation station are substantially all to be completed on the same day.Survey station spacing is generally 5 kilometers
Left and right, alpine terrain is complicated, and survey station spacing is relatively small, to guarantee the reliability of Regional Gravity Field.Every station GNSS observation time
No less than 40 minutes, the sampling interval was 30 seconds, and instrument is Leica GX1230 double frequency GNSS receiver.Earth's surface gravimetric observation is taken
The round-trip series connection observed pattern of A → B → C →...→ C → B → A form, two instrument simultaneous observations, to reduce the sight of earth's surface gravity
Transmission error is surveyed, accuracy of observation is improved.Gravity apparatus uses CG5/Burris relative gravity instrument.
The isostasy face depth h of each observation station is calculated according to above-mentioned observation data1With Crust Moho plane depth h2。
Isostasy face depth h1Calculation method：
Wherein：h*It is the height above sea level for each observation station that Leica GX1230 double frequency GNSS receiver observes.ρtFor Hai Ping
The density of the above substance in face, Δ ρ be earth mantle and crust materials density contrast, the two values be not it is fixed, different regions are different
Sample needs to be chosen according to the actual density of local lithospheric structure and earth mantle and the earth's crust, and the density of earth mantle is generally taken as 3.3g/
cm3If ρtValue is 2.67g/cm3, then Δ ρ value is 0.63g/cm3.H is that the b.s.l. of reference station surrounding area is flat
Equal CRUSTAL THICKNESS, different zones values is different, and e.g., the Sichuan Basin is about 40 kilometers, this research institute Qinghai-xizang Plateau Region take
Value is 53.6 kilometers.
Crust Moho plane depth h2Calculating it is more relative complex, calculation formula is as follows：
Wherein, it is with reference to observation point, h that we, which choose Lhasa absolute gravity station,0For the Moho face thickness of Lhasa observation station, take
Value is 69km；Δg0For the bouguer gravity anomaly at Lhasa absolute gravity station, value is -471.367mGal；Δg*For surface observation
The bouguer gravity anomaly value stood is calculated by earth's surface gravity/GNSS observation data.G is universal gravitational constant, and Δ ρ is earth mantle
It is 0.63g/cm in survey region value with the density contrast of crust materials3。
What is observed by relative gravity instrument CG5/Burris is certain earth's surface gravimetric observation station relative to reference observation point
Gravitational difference, to obtain the bouguer gravity anomaly i.e. Δ g at the surface observation station*, need to carry out data processing to observational data.Place
Steps are as follows for reason：1. the absolute gravity that each observation station is calculated by observing data information is seen on the basis of with reference to observation point
Measured value；2. carrying out altitude correction, middle layer correction and topographic correction, the bouguer gravity anomaly calculated result of each observation station is obtained,
What the bouguer gravity anomaly reflected is the variation of b.s.l. material density.It should be pointed out that being seen from earth's surface gravity/GNSS
The data processing of the bouguer gravity anomaly of measured data calculating observation point, is data processing basic known to gravity measurement personnel
Means.
The b.s.l. that table 1 gives survey region is averaged CRUSTAL THICKNESS, and in four observation stations such as Milin, rice is auspicious
The elevation observed, the isostatic surface depth being calculated, bouguer gravity anomaly, Moho depth and vertical tectonic stress knot
Fruit.
Table 1：Isostatic surface, Mohorovich discontinuity and the vertical tectonic stress calculated result of four observation stations such as Milin
According to the above-mentioned each survey station isostasy face depth h being calculated1With Crust Moho plane depth h2, counted by formula (1)
Calculation obtains each vertical tectonic stress field of observation point lithosphere in Yarlungzangbo Grand Canyon surrounding area, as shown in Figure 4.Wherein, it observes
The lithosphere vertical configuration stress value F of pointIt is verticalIt is immediately arrived at by above-mentioned calculating, the result in the place other than observation point is produced by interpolation
It is raw.Therefore, the result precision near observation station is higher, is indicated with gray scale；Apart from observation station, farther away local computational accuracy is opposite
It is lower, it is indicated with isopleth.
Fig. 4 shows that the vertical structural stress tectonic of biggish negative value occurs in Linzhi western part regional area, the result major embodiment
The influence of massif, our survey grid pass through more than 4600 meters of height above sea level of Mount Sejila mouth on the ground, and the above-mentioned vertical construction of negative value is answered
Power is exactly the response of massif self weight.It is maximum to be distributed in Motuo surrounding area, the lithosphere of this area on the occasion of vertical tectonic stress
Carry the vertical tectonic stress of up to 50MPa or so now, so vertical tectonic stress positive value of huge lithosphere should be with
The higher earth's surface denudation degree in Grand Canyon Lower Reaches is related.In addition, surveying near the Jialing Rive basin where the area northwestward, easy tribute
Lithosphere also carry the vertical tectonic stress of 30-40MPa, formation mechenism needs further to be confirmed.Milin and rice it is auspicious it
Between Grand Canyon upstream area lithosphere carry faint negative value tectonic stress, result estimation and a large amount of silt in the section
Deposit related, the ground river valley sedimentation in the presence of up to 500 meters or so now.
The vertical tectonic stress distribution of lithosphere in order to more clearly show Yarlungzangbo Grand Canyon upstream and downstream area, we
The vertical tectonic stress field data near the river of the Yarlung Zangbo River is intercepted, and they are plotted in section shown in fig. 5.Fig. 5
Water flow distance with the gravimetric observation station (endpoint of survey grid southwest shown in Fig. 3) near Milin for starting point, with river away from the point
River is launched into one and is cutd open with observation station elevation (or vertical tectonic stress value of lithosphere carrying) for the longitudinal axis for horizontal axis
Face shows the size and Orientation of the vertical tectonic stress of elevation variation and the lithosphere carrying in the section river valley in detail.Fig. 5
In, downward arrow indicates that the vertical tectonic stress of lithosphere carrying is negative value, is directed toward the earth's core；Upward arrow indicates lithosphere
The vertical tectonic stress of carrying is positive value, far from the earth's core；Grey square frame indicates Yarlungzangbo Grand Canyon position.Fig. 5 is clear
It has been shown that, the lithosphere of the Grand Canyon upstream area such as auspicious surrounding area in Milin, rice carry the downward tectonic stress of relative weak,
And carry the upward tectonic stress of up to 50MPa or so if Grand Canyon Lower Reaches such as Motuo surrounding area.
Fig. 5 gives the vertical tectonic stress of Yarlungzangbo Grand Canyon upstream and downstream section now and is distributed in detail.It obtains accordingly：
Yarlung Zangbo Grand Canyon, Tibet, China lower reache such as Motuo surrounding area is suitable for building large reservoir, because the lithosphere in the region is now
The upward vertical tectonic stress of up to 50MPa or so is carry, the self weight of reservoir filling will make lithospheric loading effect
Lithosphere is more balanced, is not easy to initiate earthquake；Grand Canyon upper reach is not suitable for building large reservoir, because of the region lithosphere
The downward vertical tectonic stress of relative weak is carry now, if as reservoir is built in Milin surrounding area at this, it is contemplated that can
The earthquake of certain scale can be triggered, regional seismic activity is increased.
The present invention is built by the vertical tectonic stress field of analysis lithosphere and large reservoir between adjoint earth's surface load
Relationship gives a kind of based on earth's surface gravity and GNSS translocation data prediction reservoir-induced earthquake from the angle of lithospheric isostasy
New method, and Applied D emonstration is unfolded in Yarlungzangbo Grand Canyon upstream and downstream area.This method can be before reservoir be formally built
The seismic risk that related reservoir filling may induce is predicted, while can provide as needed is suitable for the river for building reservoir
Section, therefore, the higher application value of this method and effectiveness of disaster reduction.
Claims (3)
1. a kind of method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake, it is characterised in that：It calculates quasi-
Build the vertical tectonic stress in reservoir location or the carrying of section lithosphere；If the vertical tectonic stress of lithosphere carrying is remotely
The heart, the location/section are suitble to build large reservoir, are not easy to trigger reservoir macroseism；If the vertical tectonic stress of lithosphere carrying
It is directed toward the earth's core, which is not suitable for construction large reservoir, there is a possibility that triggering reservoir macroseism.
2. the method according to claim 1 based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake,
It is characterized in that：The vertical tectonic stress F of the lithosphere carryingIt is verticalFor：
Wherein, h1For isostasy face depth, h2For Crust Moho plane depth, Δ ρ is the density contrast of earth mantle and crust materials, no
Different with area Δ ρ value, g is acceleration of gravity；
Isostasy face depth h is calculated according to GNSS observed result1, calculation formula is as follows：
Wherein, h*It is the surface observation station height above sea level that GNSS is observed, ρtIt is the density of the above substance in sea level, Δ ρ is earth mantle
With the density contrast of crust materials, H is that this area's b.s.l. is averaged CRUSTAL THICKNESS；
Crust Moho plane depth h is calculated according to earth's surface gravity and GNSS observation data2, calculation formula is as follows：
Wherein, G is universal gravitational constant, and Δ ρ is the density contrast of earth mantle and crust materials；h0It is the Mohorovich discontinuity depth with reference to observation point
Degree, Δ g0For the bouguer gravity anomaly with reference to observation point；Δg*For the bouguer gravity anomaly at surface observation station.
3. the method according to claim 2 based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake,
It is characterized in that：It is positive based on the vertical tectonic stress for the lithosphere carrying that earth's surface gravity and GNSS observation data determine, shows this
For the vertical tectonic stress that location or section lithosphere are born far from the earth's core, direction is upward, opposite with earth's surface gravity direction；
The vertical tectonic stress of the lithosphere carrying is negative, and shows the location or the vertical tectonic stress that section lithosphere is born
It is directed toward the earth's core, direction is downward, identical as earth's surface gravity direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520390.3A CN108919338A (en) | 2018-05-28 | 2018-05-28 | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520390.3A CN108919338A (en) | 2018-05-28 | 2018-05-28 | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108919338A true CN108919338A (en) | 2018-11-30 |
Family
ID=64409915
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810520390.3A Pending CN108919338A (en) | 2018-05-28 | 2018-05-28 | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108919338A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244352A (en) * | 2019-06-11 | 2019-09-17 | 西安石油大学 | A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070118292A1 (en) * | 2005-11-22 | 2007-05-24 | Geomechanics International, Inc. | Stress and pore pressure limits using a strength of materials approach |
| CN204649979U (en) * | 2015-02-13 | 2015-09-16 | 厦门地震勘测研究中心 | Reassembling type convertible rifle battle array focus platform |
| CN105631155A (en) * | 2016-01-12 | 2016-06-01 | 昆明理工大学 | Reservoir-induced earthquake probability calculation method |
-
2018
- 2018-05-28 CN CN201810520390.3A patent/CN108919338A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070118292A1 (en) * | 2005-11-22 | 2007-05-24 | Geomechanics International, Inc. | Stress and pore pressure limits using a strength of materials approach |
| CN204649979U (en) * | 2015-02-13 | 2015-09-16 | 厦门地震勘测研究中心 | Reassembling type convertible rifle battle array focus platform |
| CN105631155A (en) * | 2016-01-12 | 2016-06-01 | 昆明理工大学 | Reservoir-induced earthquake probability calculation method |
Non-Patent Citations (3)
| Title |
|---|
| 付广裕，等: "汶川MW7.9地震周边地区布格重力异常与岩石圈垂向构造应力场", 《地震》 * |
| 彭汉兴: "《环境工程水文地质学》", 31 May 1998 * |
| 高尚华，等: "利用重力/GPS联合观测数据计算地壳垂向构造应力的新方法", 《地球物理学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244352A (en) * | 2019-06-11 | 2019-09-17 | 西安石油大学 | A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tang et al. | Geohazards in the three Gorges Reservoir Area, China–Lessons learned from decades of research | |
| Van Camp et al. | Geophysics from terrestrial time‐variable gravity measurements | |
| CN108254782B (en) | A kind of acquisition methods and system of side slope eaerthquake damage failure probability | |
| De Natale et al. | The Campi Flegrei caldera: unrest mechanisms and hazards | |
| Pan et al. | GPS imaging of vertical bedrock displacements: quantification of two‐dimensional vertical crustal deformation in China | |
| CN108678034B (en) | Monitoring method based on cubic bridge foundation scouring safety monitoring device | |
| de Zeeuw-van Dalfsen et al. | Net gravity decrease at Askja volcano, Iceland: constraints on processes responsible for continuous caldera deflation, 1988–2003 | |
| Jousset et al. | Elastic models for the magma intrusion associated with the 2000 eruption of Usu Volcano, Hokkaido, Japan | |
| Ferguson et al. | The 4D microgravity method for waterflood surveillance: Part 3—4D absolute microgravity surveys at Prudhoe Bay, Alaska | |
| Hamiel et al. | Seismic and aseismic slip evolution and deformation associated with the 2009–2010 northern Malawi earthquake swarm, East African Rift | |
| Winterwerp et al. | Mud-induced wave damping and wave-induced liquefaction | |
| Zandt et al. | Crustal flexure associated with normal faulting and implications for seismicity along the Wasatch Front, Utah | |
| CN108918682B (en) | Entrenched valley Slope Rock Mass natural crustal stress indoor test analysis method now | |
| Itaba et al. | Continuous observation of groundwater and crustal deformation for forecasting Tonankai and Nankai earthquakes in Japan | |
| Sturkell et al. | Seismic and geodetic insights into magma accumulation at Katla subglacial volcano, Iceland: 1999 to 2005 | |
| CN108919338A (en) | Method based on earth's surface gravity and GNSS observation data prediction reservoir-induced earthquake | |
| Imaeva et al. | Seismotectonic activation of modern structures of the Siberian Craton | |
| Plateaux et al. | Implications of fracturing mechanisms and fluid pressure on earthquakes and fault slip data in the east Iceland rift zone | |
| Bagnardi | Dynamics of magma supply, storage and migration at basaltic volcanoes: Geophysical studies of the Galápagos and Hawaiian volcanoes | |
| Romano et al. | Non-isothermal momentum transfer and ground displacements rate at Campi Flegrei caldera (Southern Italy) | |
| Jousset et al. | Post-eruptive volcanic dome evolution as revealed by deformation and microgravity observations at Usu volcano (Hokkaido, Japan) | |
| Prasad et al. | Co-seismic gravity changes in the Koyna-Warna region: Implications of mass redistribution | |
| Berrino | Gravity changes and present-day dynamics of the Island of Pantelleria (Sicily Channel—Italy) | |
| Namdar | Tsunami and liquefaction resistance of subsoil | |
| CN105911591B (en) | A kind of seabed is with respect to land surveying method |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181130 |
|
| WD01 | Invention patent application deemed withdrawn after publication |