CN114004104B - CORS site selection method based on checkerboard test - Google Patents
CORS site selection method based on checkerboard test Download PDFInfo
- Publication number
- CN114004104B CN114004104B CN202111318048.3A CN202111318048A CN114004104B CN 114004104 B CN114004104 B CN 114004104B CN 202111318048 A CN202111318048 A CN 202111318048A CN 114004104 B CN114004104 B CN 114004104B
- Authority
- CN
- China
- Prior art keywords
- cors
- fracture zone
- nodes
- inversion
- 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.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 33
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 title claims abstract 16
- 238000010187 selection method Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 238000011160 research Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 6
- 201000004569 Blindness Diseases 0.000 abstract description 2
- 239000013505 freshwater Substances 0.000 description 8
- 230000000903 blocking effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 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
- 238000011081 inoculation Methods 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
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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 method for selecting a CORS site based on checkerboard testing. The CORS stations used for inversion in the traditional inter-seismic deformation research are stations on blocks around fracture zones, and blindness is achieved. The invention utilizes GIS software to symmetrically select a plurality of virtual line segments on two sides of the fracture zone; selecting a plurality of characteristic points on the virtual line segment as simulation CORS sites; two rows of nodes are arranged along the trend and the trend direction of the fracture zone; and performing forward and backward checkerboard test by adopting DEFNODE negative dislocation inversion program, wherein the smaller the color difference of the checkerboard black-white gray image is, the higher the resolution of the test node model on CORS site data and the reliability of inversion are, and taking the result as a point selection basis. The research result of the invention provides important reference and reference for the seismic fault activity researchers in the fault distribution process, and saves the expensive cost of constructing the CORS site.
Description
Technical Field
The invention belongs to the technical field of inversion inter-seismic deformation of a block negative dislocation model in geodetic and earth dynamics and seismology research, and particularly relates to a method for selecting a CORS (coordinated rotation system) station based on a checkerboard test.
Background
In addition to the wide application of the GNSS technology in the field of geology due to the advantages of high observation precision and reliability, wide observation range, all weather, real-time performance and the like in the aspect of data acquisition, the current research field is particularly deep into other subjects such as seismology, geodynamics, geophysics, geology and the like. With the continuous improvement of multimode foundation enhancement systems and the continuous improvement of the precision of multisource modification models in data processing, GNSS (Global navigation satellite System) is deeply developed in the fields of high-precision monitoring of global plate movement, regional crust deformation, inversion of earthquake dangers and the like.
For the research of earthquake dangers, especially for the inoculation stage of the earthquake, MCCAFFREY and the like, a DEFNODE negative dislocation program special for inversion of the deformation between the earthquakes is developed, a fracture zone and an adjacent block velocity field thereof obtained by calculation by using long-time observation data of a continuously operated reference station CORS (Continuously Operating Reference Stations) are utilized, and the spatial distribution characteristics of the locking degree and the sliding loss rate of each section of the earthquake fracture zone are inverted by combining a block negative dislocation model, so that corresponding earthquake dangers are analyzed.
The traditional inter-seismic deformation research is mostly aimed at a fracture zone by utilizing negative dislocation inversion, wherein CORS sites for inversion directly select sites on blocks around the fracture zone, so that the method has certain blindness, and the research preparation work is not fine enough. When the complex engineering project is aimed at, for example, the project research range is larger or the band-shaped research area is longer, the related fracture bands are more and complex, and at this time, the selection work of the CORS site is required to be finer, more reasonable, more scientific and more reliable aiming at different fracture bands.
Disclosure of Invention
In order to make up for the defects of the prior art, the invention provides a method for selecting a CORS site based on a checkerboard test, which is used for selecting a CORS site in a reasonable range.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
The method for selecting the CORS site based on the checkerboard test comprises the following steps:
step one: obtaining geological data of the fracture zone, and fitting and simplifying the geometric form of the fracture zone by using GIS software to form a simplified fracture zone segment;
Step two: in GIS software, taking a broken belt as a center, symmetrically selecting a plurality of virtual line segments on two sides of the broken belt, wherein the geometric form of the virtual line segments is the same as that of the broken belt;
step three: selecting a plurality of characteristic points on the virtual line segment determined in the second step according to a set distance and equal intervals, taking the characteristic points as simulation CORS sites, taking axisymmetric areas wrapping the selected virtual sites at the left side and the right side of a fracture zone as the range of virtual plates, and acquiring longitude and latitude information of the simulation CORS sites and the plates;
Step four: 4 nodes are arranged at equal intervals of 30km along the trend of the fracture zone, so that the arrangement of the first row of nodes is completed; selecting a first node in the second row of nodes along the direction of the tendency of the fracture zone, and setting the rest 3 nodes at equal intervals of 30km along the direction of the fracture zone to finish the setting of the second row of nodes;
Step five: performing a forward and backward tessellation test by adopting DEFNODE negative dislocation inversion program;
Step six: comparing and analyzing the locking degree of 8 nodes of the fracture zone after constraint in the forward test with the locking degree of 8 nodes of the fracture zone obtained in the inversion test, wherein the smaller the color difference of the checkerboard black-white gray-scale image is, the higher the resolution of the test node model on CORS site data and the reliability of inversion are, and taking the comparison result as a point selection basis;
Step seven: and step five, step six, testing all the simulated CORS sites, and determining the final CORS site selection range.
Specifically, the fifth step comprises the following steps:
step one: forward testing
Setting a speed field of a simulated CORS station, and forcedly restricting the locking degree of 8 nodes of a fracture zone, and acquiring a predicted value of the speed field of the simulated CORS station through DEFNODE negative dislocation inversion program by combining longitude and latitude information of the simulated CORS station;
Step two: inversion test
And taking the predicted value of the speed field of the simulated CORS station obtained in forward modeling as constraint information, and obtaining the locking degree of 8 nodes of the fracture zone through DEFNODE negative dislocation inversion program without forcibly constraining the 8 nodes of the fracture zone.
The invention has the beneficial effects that:
1) According to the invention, by using DEFNODE negative dislocation inversion program, sensitivity analysis is performed by simulating CORS sites and plates and using checkerboard test, the smaller the color difference of checkerboard black-white gray images is, the higher the resolution of the test node model on CORS site data and the reliability of inversion are, and the reliability of the method is high;
2) The research result of the invention provides important references and references for seismic fault activity researchers in the fault point distribution process, and saves the expensive cost of constructing CORS sites;
3) The method aims at the banded region, and has scientific, reasonable and reliable reference value and meaning for selecting the existing CORS site when the fault activity characteristic of the region and the dynamic reference are researched and maintained.
Drawings
FIG. 1 is a diagram of virtual site setup results;
FIG. 2 is a graph of the results of the node locking coefficients during forward modeling;
FIG. 3 is a graph of the results of the node locking coefficients during inversion;
FIG. 4 is a schematic diagram of a CORS station and a velocity field near a fresh water river fracture zone;
FIG. 5 is a schematic view of the division of adjacent blocks of fresh water river fracture zones;
FIG. 6 is a graph of the results of node locking degree in the inversion of a station in a fresh water river fracture zone of 25 km;
FIG. 7 is a graph showing the result of node locking degree in the inversion of the stations between 25km and 100km of the fracture zone of the fresh water river.
Detailed Description
The present invention will be described in detail with reference to the following embodiments.
A method for selecting CORS sites based on checkerboard test is to set a plurality of simulation sites with the same distance from a fracture zone on the basis of DEFNODE negative dislocation model inversion program, and to carry out forward and backward experiments by adopting the checkerboard test method and to carry out sensitivity tests with different distances, and comprises the following specific steps:
1. and acquiring trend data of 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 stations longitude and latitude, wherein the specific steps are as follows:
1) Fitting and simplifying the broken band shape file in GIS software to form simplified broken band segments;
2) In GIS software, taking a broken belt as a center, and symmetrically selecting a plurality of virtual line segments on two sides of the broken belt; the geometric form of the virtual line segment is the same as that of the broken belt, and the broken belt segment is copied in a translation mode;
3) Manually adding virtual sites on the newly-built virtual line segments according to a certain distance to ensure that the virtual sites are uniformly distributed on the virtual line segments;
4) The axisymmetric areas containing the selected virtual sites are used as the ranges of the virtual plates on the left side and the right side of the fracture zone, and the final effect is shown in figure 1;
5) And acquiring the longitude and latitude of all the virtual sites and the virtual plate nodes.
3. And performing forward and backward tessellation test by using the block model and the negative dislocation model through DEFNODE inversion program, and judging the setting effect of the virtual site.
The DEFNODE program performs constraint inversion on one or more data selected from GPS velocity field data, tidal data, level data, earth surface inclination rate, seismic slip vector, long-term fault slip rate, inSAR data and fault conversion azimuth angle 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 inversion model is:
In the above formula, X represents a station position; v i (X) is the speed of station X; b is the number of blocks; Δb is the region boundary of the b block model, h=1 if X is inside block b, otherwise h=0; RΩb Euler pole for block b relative to the reference block; i is a unit vector of i direction speed; f is the number of broken layers; nk is the number of knots of fault k; phi nk is the blocking value of the nth node of the fault k; g ij(X,Xnk) is a response function; hΩf The Euler pole is the lower disc relative to the upper disc; x nk is the position of the nth node on fault k; j is a trend or trend unit vector along the fault plane. The left side of the middle sign in the above formula represents the speed obtained by observation of a surface measuring station, the first part on the right side of the formula represents the movement of a block, the movement is represented by Euler poles, and the rigid body movement rate of the station on the block is obtained by using the Euler poles of the block. The second part on the right side of the formula represents the fault locking effect, the sliding rate is obtained by using the parameters hΩf and X nk, and then the sliding loss rate corresponding to the node position can be obtained by multiplying the sliding rate by the locking coefficient of the fault at the position.
According to the DEFNODE program inversion fault locking degree principle, setting a virtual plate position and carrying out experiments on the setting rationality of a virtual site by using a chessboard test, the specific steps are as follows:
(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 in the trend direction, wherein the total is two rows of nodes. In the forward modeling process, the assumed speed field and longitude and latitude information of the simulation CORS station are input, the blocking coefficient of the node is constrained forcedly, white represents complete blocking (blocking coefficient is 1), black represents complete creep (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 forcedly, and the locking degree of the 8 nodes can be obtained in the inversion process. FIG. 3 is a graph of the locking coefficients of a node model obtained by inversion of a CORS station;
And comparing and analyzing all node locking degrees after forward constraint with node locking degrees obtained by inversion, so that the resolution ratio of the node model to CORS site data and the reliability degree of inversion can be tested. And comparing the chessboard test results, wherein the smaller the color difference of the black-and-white gray-scale graph of the chessboard is, the higher the resolution of the test node model on the CORS site data and the reliability of inversion are, taking the test node model as a point selection basis, and selecting the distance with the better result as the most CORS site distance for researching the inter-seismic deformation.
Taking a fresh water river fracture zone as an example, the specific steps are as follows:
1) Firstly, acquiring the speed field of all available China continental construction environment monitoring networks (CMONOC) around a fresh water river fracture zone to observe a site, wherein the site position and the motion trend of the speed field are shown in figure 4.
2) And respectively screening available sites within 25km and between 25km and 100km around the broken zone of the fresh water river, and acquiring longitude and latitude coordinates and a horizontal speed field.
3) According to geological data, dividing adjacent blocks of a fresh water river fracture zone, dividing a research area into blocks of Ba Yan Kala and blocks of Sichuan northwest, and acquiring longitude and latitude coordinates of the blocks, as shown in fig. 5.
4) The inversion sensitivity of the CORS sites within 25km and between 25km and 100km was checked by using a checkerboard test, the procedure was identical to that of the previous embodiment, and the final results are shown in FIGS. 6 and 7. The result shows that the CORS station within 25km has better inversion effect.
The content of the invention is not limited to the examples listed, and any equivalent transformation to the technical solution of the invention that a person skilled in the art can take on by reading the description of the invention is covered by the claims of the invention.
Claims (1)
1. The method for selecting the CORS site based on the checkerboard test is characterized by comprising the following steps of: the method comprises the following steps:
step one: obtaining geological data of the fracture zone, and fitting and simplifying the geometric form of the fracture zone by using GIS software to form a simplified fracture zone segment;
Step two: in GIS software, taking a broken belt as a center, symmetrically selecting a plurality of virtual line segments on two sides of the broken belt, wherein the geometric form of the virtual line segments is the same as that of the broken belt;
step three: selecting a plurality of characteristic points on the virtual line segment determined in the second step according to a set distance and equal intervals, taking the characteristic points as simulation CORS sites, taking axisymmetric areas wrapping the selected virtual sites at the left side and the right side of a fracture zone as the range of virtual plates, and acquiring longitude and latitude information of the simulation CORS sites and the plates;
Step four: 4 nodes are arranged at equal intervals of 30km along the trend of the fracture zone, so that the arrangement of the first row of nodes is completed; selecting a first node in the second row of nodes along the direction of the tendency of the fracture zone, and setting the rest 3 nodes at equal intervals of 30km along the direction of the fracture zone to finish the setting of the second row of nodes;
Step five: performing a forward and backward tessellation test by adopting DEFNODE negative dislocation inversion program;
Step six: comparing and analyzing the locking degree of 8 nodes of the fracture zone after constraint in the forward test with the locking degree of 8 nodes of the fracture zone obtained in the inversion test, wherein the smaller the color difference of the checkerboard black-white gray-scale image is, the higher the resolution of the test node model on CORS site data and the reliability of inversion are, and taking the comparison result as a point selection basis;
Step seven: according to the fifth step and the sixth step, testing is carried out on all the simulated CORS stations, and a final CORS station selection range is determined;
The fifth step comprises the following steps:
step one: forward testing
Setting a speed field of a simulated CORS station, and forcedly restricting the locking degree of 8 nodes of a fracture zone, and acquiring a predicted value of the speed field of the simulated CORS station through DEFNODE negative dislocation inversion program by combining longitude and latitude information of the simulated CORS station;
Step two: inversion test
And taking the predicted value of the speed field of the simulated CORS station obtained in forward modeling as constraint information, and obtaining the locking degree of 8 nodes of the fracture zone through DEFNODE negative dislocation inversion program without forcibly constraining the 8 nodes of the fracture zone.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111318048.3A CN114004104B (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 CN114004104B (en) | 2021-11-09 | 2021-11-09 | CORS site selection method based on checkerboard test |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114004104A CN114004104A (en) | 2022-02-01 |
| CN114004104B true CN114004104B (en) | 2024-06-07 |
Family
ID=79928295
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111318048.3A Active CN114004104B (en) | 2021-11-09 | 2021-11-09 | CORS site selection method based on checkerboard test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114004104B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115201825B (en) * | 2022-09-16 | 2023-01-17 | 眉山环天智慧科技有限公司 | Atmospheric delay correction method in InSAR (interferometric synthetic aperture radar) inter-seismic deformation monitoring |
| CN117454578A (en) * | 2023-09-11 | 2024-01-26 | 兰州交通大学 | Fault sliding speed inversion method based on block-fault negative dislocation model |
| CN117761729A (en) * | 2024-01-26 | 2024-03-26 | 中国地震局地震预测研究所 | GNSS site layout method for monitoring fault occlusion depth |
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/CN114004104B/en active Active
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负位错模型反演结果的可靠性和稳定性分析;赵静;江在森;武艳强;刘晓霞;魏文薪;王阅兵;李强;徐晶;;大地测量与地球动力学;20130215(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114004104A (en) | 2022-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114004104B (en) | CORS site selection method based on checkerboard test | |
| CN102927934B (en) | Method for obtaining mining area earth surface three-dimensional deformation fields through single interferometric synthetic aperture radar (InSAR) interference pair | |
| Marques et al. | GPS and tectonic evidence for a diffuse plate boundary at the Azores Triple Junction | |
| Tobita et al. | 3‐D surface deformation of the 2000 Usu eruption measured by matching of SAR images | |
| Cloetingh et al. | Neotectonics and intraplate continental topography of the northern Alpine Foreland | |
| CN104729486A (en) | Bathymetric surveying method without tide observation based on quasigeoid refinement | |
| CN113341476B (en) | Method for improving ocean gravity spatial resolution based on submarine topography-gravity combination | |
| CN104200036A (en) | Method for acquiring noise models of coordinate time series of regional GPS (global positioning system) reference stations | |
| CN109521444B (en) | Self-adaptive least square fitting estimation algorithm for GPS horizontal velocity field of crustal movement | |
| CN107247927A (en) | A kind of remote sensing images coastline information extracting method and system based on K-T Transformation | |
| Cui et al. | Crustal thickness (H) and Vp/Vs ratio (κ) images beneath the central Tien Shan revealed by the H-κ-c method | |
| Sanchez et al. | A 3-D lithospheric model of the Caribbean-South American plate boundary | |
| Liu et al. | The spatial response pattern of coseismic landslides induced by the 2008 Wenchuan earthquake to the surface deformation and Coulomb stress change revealed from InSAR observations | |
| CN113960532A (en) | Microseism positioning method based on secondary positioning calculation of imaginary source | |
| Aboud et al. | Using a 3D gravity inversion technique to image the subsurface density structure in the Lunayyir volcanic field, Saudi Arabia | |
| CN113189660B (en) | Method and system for observing array type land time-varying gravity and gradient field | |
| Leng et al. | Satellite derived active-passive fusion bathymetry based on GRU model | |
| 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 | |
| CN115018973A (en) | Low-altitude unmanned-machine point cloud modeling precision target-free evaluation method | |
| van der Hoeven et al. | GPS probes the kinematics of the Vrancea seismogenic zone | |
| Yang et al. | A novel phase unwrapping method for low coherence interferograms in coal mining areas based on a fully convolutional neural network | |
| Rybin et al. | New geophysical data on the depth structure of the Tien Shan intermontane depressions | |
| CN117236518B (en) | Prediction method for geological deformation of subway along line | |
| CN113591012B (en) | Submersible vehicle positioning method and system based on multidimensional space gravity gradient information | |
| Tendero-Salmerón | Recent and active deformation structures in the central-eastern sector of the Betic Cordillera and the Alboran Sea: indentation processes and roll-back |
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 |