CN114004104A - CORS site selection method based on checkerboard test - Google Patents
CORS site selection method based on checkerboard test Download PDFInfo
- Publication number
- CN114004104A CN114004104A CN202111318048.3A CN202111318048A CN114004104A CN 114004104 A CN114004104 A CN 114004104A CN 202111318048 A CN202111318048 A CN 202111318048A CN 114004104 A CN114004104 A CN 114004104A
- Authority
- CN
- China
- Prior art keywords
- cors
- fracture zone
- inversion
- nodes
- test
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 34
- 238000010187 selection method Methods 0.000 title claims abstract description 9
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 title claims abstract 17
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000011160 research Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 6
- 201000004569 Blindness Diseases 0.000 abstract description 2
- 230000000903 blocking effect Effects 0.000 description 8
- 239000013505 freshwater Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/12—Symbolic schematics
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a CORS site selection method based on checkerboard test. The CORS station used for inversion in the traditional interseismic deformation research is selected from stations on blocks around a fracture zone, and blindness is achieved. The method comprises the steps of symmetrically selecting a plurality of virtual line segments on two sides of a fracture zone by utilizing GIS software; selecting a plurality of characteristic points on the virtual line segment as simulated CORS sites; two rows of nodes are arranged along the trend and the inclination direction of the fracture zone; and (3) performing positive and negative evolution checkerboard test by adopting a DEFNODE negative dislocation inversion program, wherein the smaller the color difference of the black-white gray-scale image of the checkerboard is, the higher the resolution and inversion reliability of the test node model on CORS site data is, and the higher the resolution and inversion reliability is taken as a point selection basis. The research result of the invention provides important reference and reference for earthquake fault activity researchers in fault distribution, and saves expensive cost for building CORS sites.
Description
Technical Field
The invention belongs to the technical field of inversion of inter-seismic deformation by a block negative dislocation model in geodetic survey, geodynamics and seismology research, and particularly relates to a CORS site selection method based on checkerboard test.
Background
In addition to the wide application in the field of geological science, the GNSS technology has the advantages of high observation precision and reliability, wide observation range, all-weather performance, real-time performance and the like in the aspect of data acquisition, and is particularly applied to other subjects such as seismology, geodynamics, geophysics, geology and the like in the current research field. With the continuous improvement of a multimode foundation enhancement system and the continuous improvement of the accuracy of a multisource modification model in data processing, the GNSS is more deeply developed in the fields of high-accuracy monitoring of global plate motion, regional crust deformation, inversion of earthquake dangerousness and the like.
For the research of earthquake dangerousness, particularly in the induction stage of an earthquake, a DEFINODE negative dislocation program specially used for inverting earthquake deformation is developed by McCaffrey and the like, a fracture zone and an adjacent block velocity field thereof are obtained by calculating long-time observation data of a continuous operation Reference station CORS (continuous Operating Reference states), and the block negative dislocation model is combined to invert the locking degree of each section of the earthquake discontinuous fracture zone and the spatial distribution characteristics of the sliding loss rate, so that corresponding earthquake danger analysis is carried out.
Most of traditional inter-seismic deformation inversion researches are carried out on a fracture zone by utilizing negative dislocation, wherein the CORS station for inversion directly selects stations on blocks around the fracture zone, certain blindness is achieved, and research preparation work is not fine enough. When aiming at complex engineering projects, such as large project research range or long strip-shaped research area, more and complex fracture zones are involved, and at the moment, the selection work of the CORS sites needs to be more precise, reasonable, scientific and reliable aiming at different fracture zones.
Disclosure of Invention
In order to make up for the defects of the prior art, the invention provides a CORS site selection method based on checkerboard test, which is used for selecting CORS sites in a reasonable range.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the CORS site selection method based on the checkerboard test comprises the following steps:
the method comprises the following steps: acquiring geological data of a fracture zone, and fitting and simplifying the geometric form of the fracture zone by using GIS software to form a simplified fracture zone line segment;
step two: in GIS software, a plurality of virtual line segments are symmetrically selected on two sides of a fracture zone by taking the fracture zone as a center, and the geometric form of the virtual line segments is the same as that of the fracture zone;
step three: selecting a plurality of feature points on the virtual line segment determined in the second step according to the set distance and the equal interval, taking the feature points as simulated CORS sites, taking the axial symmetric regions wrapping the selected virtual sites as the range of virtual plates on the left side and the right side of the fracture zone, and acquiring longitude and latitude information of the simulated CORS sites and the plates;
step four: arranging 4 nodes at equal intervals of 30km along the trend of the fracture zone to complete the arrangement of the first row of nodes; selecting a first node in a second row of nodes along the direction of the inclination of the fracture zone, and arranging the other 3 nodes at equal intervals of 30km along the trend of the fracture zone to complete the arrangement of the nodes in the second row;
step five: carrying out positive and negative error checkerboard test by using a DEFNODE negative dislocation inversion program;
step six: comparing and analyzing the latching degree of 8 nodes of the fracture zone after constraint in the forward test with the latching degree of 8 nodes of the fracture zone obtained in the inversion test, wherein the smaller the color difference of the black-white checkerboard gray-scale image is, the higher the resolution and inversion reliability of the CORS site data of the test node model are, and the higher the resolution and inversion reliability is taken as a point selection basis;
step seven: and testing all simulated CORS sites according to the fifth step and the sixth step, and determining the selection range of the final CORS sites.
Specifically, the step five comprises the following steps:
the method comprises the following steps: forward test
Setting a speed field of a simulated CORS site, forcibly restricting the locking degree of 8 nodes of a fracture zone, and obtaining a predicted value of the speed field of the simulated CORS site through a DEFNODE negative dislocation inversion program by combining longitude and latitude information of the simulated CORS site;
step two: inversion test
And taking the predicted value of the velocity field of the simulated CORS station obtained in the forward modeling as constraint information, not performing forced constraint on 8 nodes of the fractured zone, and obtaining the locking degree of the 8 nodes of the fractured zone through a DEFINODE negative dislocation inversion program.
The invention has the beneficial effects that:
1) according to the method, a DEFNODE negative dislocation inversion program is utilized, sensitivity analysis is carried out by simulating CORS sites and plates and utilizing a checkerboard test, the smaller the color difference of a checkerboard black-white gray-scale image is, the higher the resolution and inversion reliability of a test node model to CORS site data is, and the reliability of the method is high;
2) the research result of the invention provides important reference and reference for earthquake fault activity researchers in fault distribution, and saves expensive cost for building CORS sites;
3) the method aims at the banded regions, and has scientific, reasonable and reliable reference value and significance for selecting the existing CORS sites when the fault activity characteristics and the dynamic benchmark of the regions are researched to be maintained.
Drawings
FIG. 1 is a diagram of virtual site setup results;
FIG. 2 is a graph of node blocking coefficient results during forward evolution;
FIG. 3 is a graph of node locking coefficient results during inversion;
FIG. 4 is a schematic diagram of a CORS site and velocity field near the fresh water river fault zone;
FIG. 5 is a schematic diagram of dividing adjacent blocks of a fresh water river fracture zone;
FIG. 6 is a result graph of node locking degree in inversion of sites within 25km of fresh water river fracture zone;
FIG. 7 is a graph showing the results of node locking degree in the inversion of a fresh water river fault zone between 25km and 100 km.
Detailed Description
The present invention will be described in detail with reference to specific embodiments.
A CORS site selection method based on a checkerboard test is characterized in that a plurality of simulation sites with the same distance from a fracture zone are arranged on the basis of a DEFNODE negative dislocation model inversion program, a checkerboard test method is adopted to carry out a forward and backward experiment, and sensitivity tests with different distances are carried out, and the method specifically comprises the following steps:
1. and acquiring the trend data of the characteristic nodes of the fracture zone to be researched to form a shape format file.
2. Setting a series of CORS stations with different distances from a fracture zone to obtain a series of virtual CORS station longitudes and latitudes, and specifically comprising the following steps of:
1) fitting and simplifying a shape file of the fracture zone in GIS software to form a simplified fracture zone segment;
2) in GIS software, a plurality of virtual line segments are symmetrically selected on two sides of a fracture zone by taking the fracture zone as a center; the geometric form of the virtual line segment is the same as that of the fracture zone, and the fracture zone line segment is copied in a translation mode;
3) manually adding virtual sites on the newly-built virtual line segment according to a certain distance, so that the virtual sites are uniformly distributed on the virtual line segment;
4) the final effect is shown in fig. 1, with the axially symmetric region containing the selected virtual site as the extent of the virtual plate on the left and right sides of the fracture zone;
5) and acquiring the longitude and latitude of all the virtual sites and the virtual plate nodes.
3. And performing positive and negative error check board test by using the block model and the negative dislocation model through a DEFNODE inversion program, and judging the setting effect of the virtual station.
The DEFNODE program inputs data including one or more of GPS velocity field data, tide data, level data, earth surface inclination velocity, seismic slip vector, fault long-term slip velocity, InSAR data and fault conversion azimuth angle for constraint inversion during inversion. In this patent, the GPS horizontal velocity field is selected as the inversion constraint because the GPS vertical velocity field is not accurate enough in the inversion process.
If the uniform strain inside the block is not considered, the formula of the inverse model is as follows:
in the above formula, X represents the station position; vi(X) is the velocity of station X; b is the number of blocks; Δ b is the region boundary of the b block model, if X is inside the block b, H is 1, otherwise H is 0;RΩbeuler pole for block b relative to a reference block; i is a unit vector of the speed in the i direction; f is the number of broken layers; nk is the number of sections of the fault k; phi is ankIs the latching value of the nth node of the fault k; gij(X,Xnk) Is a response function;hΩfthe Euler pole of the lower disc relative to the upper disc; xnkIs the position of the nth node on the fault k; j is the unit vector of the trend or dip along the fault plane. The left side of the middle sign in the upper formula represents the observed speed of the earth surface measuring station, the first part on the right side of the formula represents the motion of the block and is represented by an Euler pole, and the Euler pole of the block is used for obtaining the rigid motion rate of the station on the block. The second part on the right side of the formula represents fault blocking effect and is firstly measured by parametershΩfAnd XnkAnd obtaining the sliding rate, and multiplying the sliding rate by the blocking coefficient of the fault at the position to obtain the sliding loss rate corresponding to the node position.
According to the DEFNODE program inversion fault blocking degree principle, virtual plate positions are set, and a chessboard test is used for testing the setting rationality of virtual stations, and the method specifically comprises the following steps:
(1) firstly, selecting one node every 30km along the trend of a fracture zone, setting 4 nodes in total, and selecting one node every 30km along the inclination direction, and two rows of nodes in total. In the forward modeling process, the assumed speed field and longitude and latitude information of a simulation CORS site are input, the blocking coefficient of a constraint node is forced, white represents complete blocking (the blocking coefficient is 1), black represents complete creeping (the blocking coefficient is 0), and the forward modeling result is shown in FIG. 2;
(2) in the inversion process, the predicted value of the velocity field obtained by the model in forward modeling is used as constraint information, the 8 nodes are not constrained forcibly, and the locking degree of the 8 nodes can be obtained in the inversion process. FIG. 3 is a locking coefficient of a node model obtained by inversion at a CORS station;
and comparing and analyzing the locking degrees of all the nodes after the forward modeling constraint with the locking degrees of the nodes obtained by inversion, so that the resolution and the inversion reliability of the node model on CORS site data can be tested. Comparing the chessboard test result, the smaller the color difference of the black and white gray level image of the chessboard is, the higher the resolution and inversion reliability of the test node model to CORS station data is, and the better distance of the selected result is used as the point selection basis for researching the best CORS station distance of the deformation between the earthquakes.
Taking the fresh water river fracture zone as an example, the specific steps are as follows:
1) firstly, the velocity field of all available China continental construction environment monitoring networks (CMONOC) observation stations near the fresh water river fault zone is obtained, and the station position and the velocity field movement trend are shown in figure 4.
2) Respectively screening available stations within 25km and between 25km and 100km around the fresh water river fracture zone to obtain longitude and latitude coordinates and a horizontal direction velocity field of the stations.
3) According to geological data, dividing adjacent blocks of the fresh water river fracture zone, dividing a research area into Bayan karya blocks and northwest Sichuan blocks, and acquiring longitude and latitude coordinates of the blocks, as shown in fig. 5.
4) The inversion sensitivities of CORS sites within 25km and CORS sites between 25km and 100km were examined using a checkerboard test, the procedure being identical to that of the previous embodiment, and the final results are shown in FIGS. 6 and 7. The result shows that CORS sites within 25km have better inversion effect.
The invention is not limited to the examples, and any equivalent changes to the technical solution of the invention by a person skilled in the art after reading the description of the invention are covered by the claims of the invention.
Claims (2)
1. The CORS site selection method based on the checkerboard test is characterized by comprising the following steps of: the method comprises the following steps:
the method comprises the following steps: acquiring geological data of a fracture zone, and fitting and simplifying the geometric form of the fracture zone by using GIS software to form a simplified fracture zone line segment;
step two: in GIS software, a plurality of virtual line segments are symmetrically selected on two sides of a fracture zone by taking the fracture zone as a center, and the geometric form of the virtual line segments is the same as that of the fracture zone;
step three: selecting a plurality of feature points on the virtual line segment determined in the second step according to the set distance and the equal interval, taking the feature points as simulated CORS sites, taking the axial symmetric regions wrapping the selected virtual sites as the range of virtual plates on the left side and the right side of the fracture zone, and acquiring longitude and latitude information of the simulated CORS sites and the plates;
step four: arranging 4 nodes at equal intervals of 30km along the trend of the fracture zone to complete the arrangement of the first row of nodes; selecting a first node in a second row of nodes along the direction of the inclination of the fracture zone, and arranging the other 3 nodes at equal intervals of 30km along the trend of the fracture zone to complete the arrangement of the nodes in the second row;
step five: carrying out positive and negative error checkerboard test by using a DEFNODE negative dislocation inversion program;
step six: comparing and analyzing the latching degree of 8 nodes of the fracture zone after constraint in the forward test with the latching degree of 8 nodes of the fracture zone obtained in the inversion test, wherein the smaller the color difference of the black-white checkerboard gray-scale image is, the higher the resolution and inversion reliability of the CORS site data of the test node model are, and the higher the resolution and inversion reliability is taken as a point selection basis;
step seven: and testing all simulated CORS sites according to the fifth step and the sixth step, and determining the selection range of the final CORS sites.
2. The CORS site selection method based on checkerboard testing as claimed in claim 1, wherein: the fifth step comprises the following steps:
the method comprises the following steps: forward test
Setting a speed field of a simulated CORS site, forcibly restricting the locking degree of 8 nodes of a fracture zone, and obtaining a predicted value of the speed field of the simulated CORS site through a DEFNODE negative dislocation inversion program by combining longitude and latitude information of the simulated CORS site;
step two: inversion test
And taking the predicted value of the velocity field of the simulated CORS station obtained in the forward modeling as constraint information, not performing forced constraint on 8 nodes of the fractured zone, and obtaining the locking degree of the 8 nodes of the fractured zone through a DEFINODE negative dislocation inversion program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111318048.3A CN114004104A (en) | 2021-11-09 | 2021-11-09 | CORS site selection method based on checkerboard test |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111318048.3A CN114004104A (en) | 2021-11-09 | 2021-11-09 | CORS site selection method based on checkerboard test |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114004104A true CN114004104A (en) | 2022-02-01 |
Family
ID=79928295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111318048.3A Pending CN114004104A (en) | 2021-11-09 | 2021-11-09 | CORS site selection method based on checkerboard test |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114004104A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115201825A (en) * | 2022-09-16 | 2022-10-18 | 眉山环天智慧科技有限公司 | Atmospheric delay correction method in InSAR (interferometric synthetic aperture radar) inter-seismic deformation monitoring |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7254091B1 (en) * | 2006-06-08 | 2007-08-07 | Bhp Billiton Innovation Pty Ltd. | Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs |
WO2021008630A1 (en) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | Method, apparatus, and system for quantitative analysis of sealing strength of trap fault |
-
2021
- 2021-11-09 CN CN202111318048.3A patent/CN114004104A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7254091B1 (en) * | 2006-06-08 | 2007-08-07 | Bhp Billiton Innovation Pty Ltd. | Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs |
WO2021008630A1 (en) * | 2019-07-18 | 2021-01-21 | 中国石油天然气股份有限公司 | Method, apparatus, and system for quantitative analysis of sealing strength of trap fault |
Non-Patent Citations (1)
Title |
---|
赵静;江在森;武艳强;刘晓霞;魏文薪;王阅兵;李强;徐晶;: "Defnode负位错模型反演结果的可靠性和稳定性分析", 大地测量与地球动力学, no. 01, 15 February 2013 (2013-02-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115201825A (en) * | 2022-09-16 | 2022-10-18 | 眉山环天智慧科技有限公司 | Atmospheric delay correction method in InSAR (interferometric synthetic aperture radar) inter-seismic deformation monitoring |
CN115201825B (en) * | 2022-09-16 | 2023-01-17 | 眉山环天智慧科技有限公司 | Atmospheric delay correction method in InSAR (interferometric synthetic aperture radar) inter-seismic deformation monitoring |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Draxler | Sensitivity of a trajectory model to the spatial and temporal resolution of the meteorological data during CAPTEX | |
Marques et al. | GPS and tectonic evidence for a diffuse plate boundary at the Azores Triple Junction | |
Wang et al. | A stable reference frame for landslide monitoring using GPS in the Puerto Rico and Virgin Islands region | |
Zanutta et al. | Monitoring geodynamic activity in the Victoria Land, East Antarctica: Evidence from GNSS measurements | |
Tobita et al. | 3‐D surface deformation of the 2000 Usu eruption measured by matching of SAR images | |
CN106443765A (en) | Urban engineering seismic detection comprehensive imaging method based on multi-component observation system | |
CN109521444B (en) | Self-adaptive least square fitting estimation algorithm for GPS horizontal velocity field of crustal movement | |
CN105242328B (en) | The determination method and device of ancient hot Lithospheric Thickness | |
CN106199662A (en) | A kind of Big Dipper high-precision locating method based on analytic hierarchy process (AHP) | |
CN114004104A (en) | CORS site selection method based on checkerboard test | |
CN115964915A (en) | Kelatong rock-ring three-dimensional material framework tracing method | |
Sanchez et al. | A 3-D lithospheric model of the Caribbean-South American plate boundary | |
Ezati et al. | Tectonic activity level evaluation using geomorphic indices in the Shekarab Mountains, Eastern Iran | |
Fojtíková et al. | Quantifying capability of a local seismic network in terms of locations and focal mechanism solutions of weak earthquakes | |
CN113189660B (en) | Method and system for observing array type land time-varying gravity and gradient field | |
CN109033459B (en) | Spatial gridding construction method for soil volume weight data | |
Agata et al. | Detectability analysis of interplate fault slips in the Nankai subduction thrust using seafloor observation instruments | |
McElfresh et al. | Stress modeling of tectonic blocks at Cape Kamchatka, Russia using principal stress proxies from high-resolution SAR: new evidence for the Komandorskiy Block | |
CN113834875A (en) | Elastic wave tomography detection method and system based on three-dimensional hexagonal body survey line arrangement | |
CN216013697U (en) | Array type land time-varying gravity and gradient field observation device | |
Premužić et al. | Accuracy analysis of point velocities determined by different software packages and GNSS measurement processing methods | |
Shengtao et al. | Preliminary horizontal co-seismic displacements caused by the 2023 Mw 7.8 and Mw 7.5 Türkiye earthquakes estimated using high-rate GPS observations | |
CN114779352B (en) | Geologic body survey system based on air-ground-well multi-dimensional electromagnetic method | |
Zahran et al. | Geodynamics implication of GPS and satellite altimeter and gravity observations to the Eastern Mediterranean | |
CN108132486A (en) | A kind of heavy magnetic gradient combines the Optimized Simulated annealing method of Interface Inversion with seismic data |
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 |