CN113325473B - Vibration region determination method and system - Google Patents

Vibration region determination method and system Download PDF

Info

Publication number
CN113325473B
CN113325473B CN202110592068.3A CN202110592068A CN113325473B CN 113325473 B CN113325473 B CN 113325473B CN 202110592068 A CN202110592068 A CN 202110592068A CN 113325473 B CN113325473 B CN 113325473B
Authority
CN
China
Prior art keywords
vibration
area
region
vertexes
seismic source
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
Application number
CN202110592068.3A
Other languages
Chinese (zh)
Other versions
CN113325473A (en
Inventor
郑凡
凤一鸣
杨大鹏
杨泓渊
刘铭扬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN202110592068.3A priority Critical patent/CN113325473B/en
Publication of CN113325473A publication Critical patent/CN113325473A/en
Application granted granted Critical
Publication of CN113325473B publication Critical patent/CN113325473B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention belongs to the seismic positioning technology, in particular to a vibration area judgment method and a system, which are used for acquiring the seismic wave velocity in a monitoring range; acquiring arrival time of one or more rectangular vibration region vertexes in a monitoring range, and obtaining the distance difference between every two vertexes and an unknown seismic source in the same vibration region according to the seismic wave velocity and arrival time differences between every two vertexes; calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region; comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region; and judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitudes. The invention can judge whether the vibration generated by an unknown seismic source occurs in the vibration region by calculating the corresponding area of the rectangular vibration region. And according to the condition of each vibration area, whether the vibration occurs in the whole monitoring range can be known. The method is simple and efficient.

Description

Vibration region determination method and system
Technical Field
The invention belongs to the field of seismic positioning technology, and particularly relates to a vibration area determination method and a vibration area determination system.
Background
With the development of science and technology, the vibration positioning technology is no longer limited to the monitoring of natural earthquakes. The method has wide application in other fields such as military, civil security and the like. There are many methods of seismic localization, such as a target signal time of arrival (TOA) based localization method, an object signal angle of arrival (AOA) based localization method, and a target signal time difference of arrival (TDOA) based localization method, which are commonly used for passive localization, and some methods based on machine learning. In some seismic localization applications, it is desirable to determine whether the seismic source is within the observed seismic region. For example, in some systems, it may not make sense that the source is outside the monitoring range. It is therefore essential to determine whether the vibrations generated by the unknown seismic source occur within or outside the monitoring range. When the vibration target is judged to be out of the monitoring range, the system can only prompt or not perform any treatment. When a vibration target comes within the monitoring range, other operations of the system are required.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a method and a system for determining a vibration region, which are used to determine whether the vibration generated by an unknown seismic source occurs in a corresponding monitoring range.
The present invention is achieved in such a way that,
a vibration region determination method, the method comprising:
acquiring the seismic wave velocity in the monitoring range;
the arrival time at the vertexes of one or more rectangular vibration regions in the monitoring range is acquired,
obtaining the distance difference between every two vertexes and an unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region;
comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region;
and judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitudes.
Further, when the area of one triangle is larger than half of the area of the vibration region, judging that the vibration generated by the unknown seismic source occurs outside the vibration region;
when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with the area equal to half of the area of the vibration region exists, the vibration generated by the unknown seismic source is judged to be generated on the boundary of the vibration region.
Further, when the vibration generated by the unknown seismic source is judged to be outside the vibration region, the arrival time of four vertexes of the vibration region is judged, two vertexes with the shortest arrival time are obtained, the sides of a rectangle where the two vertexes are located are determined, and the direction of the unknown seismic source outside the region is determined.
Further, calculating the area of the triangle formed by the shock point and the two adjacent vertexes by using the distance difference and the side length of the corresponding shock region comprises calculating the area of the triangle by using a Helen formula, and converting the area of the triangle into a distance difference form.
And further, splicing and covering a plurality of vibration areas with the same size in a monitoring range, and after splicing, independently judging whether the vibration of the unknown seismic source occurs in the local vibration area or not by each vibration area according to a calculation result, wherein if the vibration of the unknown seismic source occurs in the local vibration area, the unknown seismic source is positioned in the local vibration area.
A vibration region determination system, the system comprising:
the wave velocity acquisition module is used for acquiring the seismic wave velocity in the monitoring range;
the arrival time acquisition module is used for acquiring arrival time at the vertexes of one or more rectangular vibration areas in the monitoring range,
the processing module is used for obtaining the distance difference between each two vertexes and the unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region; comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region; and judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitudes.
Further, the processing module is further configured to determine that the vibration generated by the unknown seismic source occurs outside the vibration region when there is a triangle whose area is greater than half of the area of the vibration region to which the triangle belongs; when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with the area equal to half of the area of the vibration region exists, the vibration generated by the unknown seismic source is judged to be generated on the boundary of the vibration region.
Further, the processing module is further configured to, when it is determined that the vibration generated by the unknown seismic source occurs outside the vibration region, determine the arrival time of four vertices of the vibration region, obtain two vertices with the shortest arrival time, determine the sides of the rectangle where the two vertices are located, and determine the direction of the unknown seismic source outside the region.
Further, calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region comprises: the area of the triangle is calculated by adopting a Helen formula, and then the triangle is converted into a distance difference form.
Compared with the prior art, the invention has the beneficial effects that:
the invention can judge whether the vibration generated by an unknown seismic source occurs in the vibration region by calculating the corresponding area of the rectangular vibration region. And according to the condition of each vibration area, whether the vibration occurs in the whole monitoring range can be known. The method can be used independently as a ground vibration positioning method, and can also be used in combination with other ground vibration positioning algorithms to judge whether the vibration generated by the unknown seismic source is in a corresponding monitoring range. It may be determined in which seismic region the unknown source seismic vibrations specifically occur, if they occur within the respective monitoring range. A suitable method for determining the vibration region often determines the reliability of many vibration location systems, making them more excellent in practical applications.
Drawings
FIG. 1 is a schematic diagram of a target vibration in a vibration region according to an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating a target vibration occurring outside a vibration region according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a vibration region splicing method according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The embodiment of the invention provides a vibration region judgment method, which comprises the following steps:
acquiring the seismic wave velocity in the monitoring range;
acquiring arrival time at the vertexes of one or more rectangular vibration areas in a monitoring range, wherein the arrival time refers to the time for transmitting generated seismic waves to the vertexes of the vibration areas when an earthquake occurs, and the vibration areas can be one or a plurality of rectangles with the same size or different sizes; one or more are arranged in the monitoring range. When the size is the same, the plurality of vibration regions form an array structure. The number of rectangles in the longitudinal and transverse directions of the formed array may be the same or different.
Obtaining the distance difference between every two vertexes and an unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes; the same calculation is carried out on other vibration areas, and the distance difference can be obtained by the product of the seismic wave velocity and the arrival time difference between every two seismic waves in one vibration area.
Calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region;
comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region;
and judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitudes.
The method needs to utilize a ground acquisition station to acquire data of the vibration signal, and calculates and processes the vibration signal generated by the unknown seismic source according to the data. And (3) placing a detector at each vertex as a collection channel, and dividing the whole area into two parts by using a collection channel array: outside the monitoring range and within the monitoring range. It is necessary to determine whether the vibrations generated by the unknown seismic source occur outside the monitoring range or in a particular vibration region within the monitoring range. It is defined here that every four acquisition traces (geophones) constitute a rectangular seismic region. The four acquisition traces are located at the vertices of a matrix.
Calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region, wherein the area comprises the following steps: the area of the triangle is calculated by adopting a Helen formula, and then the triangle is converted into a distance difference form. And calculating the areas of the four triangles by utilizing a Helen formula and converting the areas into a form of distance difference between two sides of the triangle. And multiplying the arrival time difference between the acquisition channels by the seismic wave velocity to obtain the distance difference between each two acquisition channels and the seismic source. The distance difference is substituted into the area of four triangles, and the logic is used for judging whether the unknown seismic source vibration occurs in the vibration region or not and which vibration region the unknown seismic source vibration occurs in. The Helen formula converts the area of four triangles into the difference between the distances of two sides of a triangle as follows: as shown with reference to figure 1 of the drawings,
four line segments AO, BO, EO and DO with lengths a, B, E and D are respectively made for connecting the collection tracks A, B, E and D of each vertex of the rectangle with the vibration occurrence point O (hypothesis), and two side lengths AB (DE) and c (c) of the rectangle1,AD=BE=c2Is an arbitrary value;
to OH1Perpendicular rectangle AD is bordered by H1、OH2Vertical rectangle BE is bordered by H2As shown in FIG. 1, let DH1=EH2=x,AH1=BH2=c2-x. The following can be obtained:
d2-x2=a2-(c2-x)2 (1)
e2-x2=b2-(c2-x)2 (2)
then (1) - (2) get:
a2-b2=d2-e2 (3)
using a sum-difference product formula and arranging to obtain:
Figure BDA0003089969250000051
from the formula (4):
Figure BDA0003089969250000052
and a is expressed by the difference of the distance between each vertex acquisition track and the seismic source, namely the difference of a, b, d and e:
Figure BDA0003089969250000061
b is expressed as the difference between a and a, b:
b=a+(b-a) (7)
from the Helen formula, the area of the triangular AOB is:
Figure BDA0003089969250000062
substituting (7) into (8) to obtain:
Figure BDA0003089969250000063
finishing to obtain:
Figure BDA0003089969250000064
substituting (6) into (10) for 2a yields:
Figure BDA0003089969250000065
finishing to obtain:
Figure BDA0003089969250000066
this time equation (12) has the triangular area as a distance difference, and is obtained from each arrival time difference obtained by multiplying the measured wave velocity v by the arrival time difference: d-a, e-a, b-a, a-b, d-e five distance differences (constant) are taken into the formula (12) to obtain the area of the triangular AOB;
similarly, the areas of the other three triangular DOEs, AODs and BOEs can be calculated;
using the area of the four triangles and half of the area of the rectangle
Figure BDA0003089969250000067
The comparison is carried out in such a way that,
when the area of one triangle is larger than half of the area of the vibration region, judging that the vibration generated by the unknown seismic source occurs outside the vibration region;
when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with the area equal to half of the area of the vibration region exists, the vibration generated by the unknown seismic source is judged to be generated on the boundary of the vibration region.
When the vibration generated by the unknown seismic source is judged to be outside the vibration area, the arrival time of four vertexes of the vibration area is judged, two vertexes with the shortest arrival time are obtained, the sides of a rectangle where the two vertexes are located are determined, and the direction of the unknown seismic source outside the vibration area is determined.
Specifically, the method comprises the following steps: when four triangles are all in the same direction
Figure BDA0003089969250000071
Judging that the vibration generated by the unknown seismic source occurs in the vibration area;
when a triangle exists, the
Figure BDA0003089969250000072
Judging that the vibration generated by the unknown seismic source occurs on the boundary of the vibration area;
when a triangle exists, the
Figure BDA0003089969250000073
And judging that the vibration generated by the unknown seismic source occurs outside the vibration area.
As shown in fig. 2, the unknown sources are calculated in the same way inside and outside the seismic region, and only the area values obtained are different.
When the vibration generated by the unknown seismic source occurs outside the vibration region, two acquisition tracks with the smallest arrival time can be found by comparing the arrival times of the four acquisition track data. And determining the sides of the two rectangles, and further determining the directions of the unknown seismic sources outside the area, namely the directions outside the vibration area, namely the directions outside the edges (upper, lower, left and right) of the vibration area.
A certain monitoring range can be covered by splicing the vibration areas (two adjacent vibration areas can share two acquisition channels), as shown in fig. 3. After splicing, each vibration area independently judges whether the vibration of the unknown seismic source occurs in the vibration area according to the calculation result. If yes, the unknown seismic source is located in the local seismic region, and if not, no processing is performed.
The judgment method for the invasion vibration target in the specific vibration area is as follows:
1. when all vibration regions judge that the unknown seismic source is outside the self region, the unknown seismic source appears outside the monitoring range;
2. when one vibration area determines that the unknown seismic source is in the self area, and all other vibration areas determine that the unknown seismic source is out of the self area, the unknown seismic source is in the vibration area;
3. when two vibration regions determine that the unknown seismic source is on one side of the self region, and all other vibration regions determine that the unknown seismic source is outside the self unit, the unknown seismic source appears on a boundary line of the two vibration regions.
The method can be used as a ground vibration positioning method for independent use, and can also be used in combination with other ground vibration positioning algorithms for judging whether the vibration generated by the unknown seismic source is in a corresponding monitoring range. It may be determined in which seismic region the unknown source seismic vibrations specifically occur, if they occur within the respective monitoring range. A suitable method for determining the vibration region often determines the reliability of many vibration location systems, making them more excellent in practical applications.
The present invention provides a vibration region determination system, including:
the wave velocity acquisition module is used for acquiring the seismic wave velocity in the monitoring range;
and the arrival time acquisition module acquires arrival time at the vertexes of one or more rectangular vibration areas in the monitoring range and acquires signals of acquisition channels.
The processing module is used for obtaining the distance difference between each two vertexes and the unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region; comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region; and judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitudes.
The processing module is further used for judging that the vibration generated by the unknown seismic source is out of the vibration region when the area of one triangle is larger than half of the area of the vibration region; when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with the area equal to half of the area of the vibration region exists, the vibration generated by the unknown seismic source is judged to be generated on the boundary of the vibration region.
The processing module is further used for judging the arrival time of four vertexes of the vibration area when the vibration generated by the unknown seismic source is judged to be outside the vibration area, obtaining two vertexes with the shortest arrival time, determining the sides of the rectangle where the two vertexes are located, and determining the direction of the unknown seismic source outside the vibration area.
Calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region, wherein the area comprises the following steps: the area of the triangle is calculated by adopting a Helen formula, and then the triangle is converted into a distance difference form.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (3)

1. A vibration region determination method, comprising:
acquiring the seismic wave velocity in the monitoring range;
the arrival time at the vertexes of one or more rectangular vibration regions in the monitoring range is acquired,
obtaining the distance difference between every two vertexes and an unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region;
comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region;
judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitude;
when the area of one triangle is larger than half of the area of the vibration region, judging that the vibration generated by the unknown seismic source occurs outside the vibration region;
when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with an area equal to half of the area of the vibration region exists, judging that the vibration generated by the unknown seismic source occurs on the boundary of the vibration region;
when the vibration generated by the unknown seismic source is judged to be outside the vibration region, judging the arrival time of four vertexes of the vibration region, obtaining two vertexes with the shortest arrival time, determining the sides of a rectangle where the two vertexes are located, and determining the direction of the unknown seismic source outside the region;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region, wherein the area comprises the following steps: the area of the triangle is calculated by adopting a Helen formula, and then the triangle is converted into a distance difference form.
2. The method of claim 1, wherein a monitoring area is covered by a plurality of vibration regions of the same size in a splicing manner, and after the splicing, each vibration region independently determines whether the vibration of the unknown seismic source occurs in the local vibration region according to the calculation result, and if so, the unknown seismic source is located in the local vibration region.
3. A vibration region determination system, characterized by comprising:
the wave velocity acquisition module is used for acquiring the seismic wave velocity in the monitoring range;
the arrival time acquisition module is used for acquiring arrival time at the vertexes of one or more rectangular vibration areas in the monitoring range,
the processing module is used for obtaining the distance difference between each two vertexes and the unknown seismic source in the same vibration region according to the seismic wave velocity and the arrival time difference between every two vertexes;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region; comparing the areas of the four triangles in the same vibration region with half of the area of the vibration region; judging whether the vibration generated by the unknown seismic source occurs outside or inside the vibration region according to the comparison result of the magnitude;
the processing module is further used for judging that the vibration generated by the unknown seismic source is out of the vibration region when the area of one triangle is larger than half of the area of the vibration region; when the areas of the four triangles are smaller than half of the area of the vibration area, judging that the vibration generated by the unknown seismic source occurs in the vibration area; when a triangle with an area equal to half of the area of the vibration region exists, judging that the vibration generated by the unknown seismic source occurs on the boundary of the vibration region;
the processing module is further used for judging the arrival time of four vertexes of the vibration area when the vibration generated by the unknown seismic source is judged to be outside the vibration area, obtaining two vertexes with the shortest arrival time, determining the sides of a rectangle where the two vertexes are located, and determining the direction of the unknown seismic source outside the vibration area;
calculating the area of a triangle formed by the shock point and two adjacent vertexes by using the distance difference and the side length of the corresponding shock region, wherein the area comprises the following steps: the area of the triangle is calculated by adopting a Helen formula, and then the triangle is converted into a distance difference form.
CN202110592068.3A 2021-05-28 2021-05-28 Vibration region determination method and system Active CN113325473B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110592068.3A CN113325473B (en) 2021-05-28 2021-05-28 Vibration region determination method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110592068.3A CN113325473B (en) 2021-05-28 2021-05-28 Vibration region determination method and system

Publications (2)

Publication Number Publication Date
CN113325473A CN113325473A (en) 2021-08-31
CN113325473B true CN113325473B (en) 2022-03-22

Family

ID=77422159

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110592068.3A Active CN113325473B (en) 2021-05-28 2021-05-28 Vibration region determination method and system

Country Status (1)

Country Link
CN (1) CN113325473B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634476A (en) * 2013-12-10 2014-03-12 南京信大气象装备有限公司 Positioning method based on recording device of mobile phone
CN103700154A (en) * 2013-12-27 2014-04-02 付诚 Beidou satellite-based special article positioning and supervising method and system
CN107036770A (en) * 2017-04-18 2017-08-11 浙江理工大学 The leakage detection and localization method of air cooler fin tube bank
CN107533409A (en) * 2015-06-19 2018-01-02 株式会社艾拉博 Coordinate detecting device and coordinate detection method
CN107765294A (en) * 2017-10-20 2018-03-06 重庆泛嘉地质勘查有限公司 Method of seismic source monitoring
CN110084986A (en) * 2019-04-29 2019-08-02 西人马(厦门)科技有限公司 A kind of circumference safety protection method and device
CN111044975A (en) * 2019-12-10 2020-04-21 北京无线电计量测试研究所 Method and system for positioning earth vibration signal
KR102218181B1 (en) * 2020-11-06 2021-02-23 케이아이티밸리(주) Method for determining real-time alarm for local target earthquake early warning using mechanical learning

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634476A (en) * 2013-12-10 2014-03-12 南京信大气象装备有限公司 Positioning method based on recording device of mobile phone
CN103700154A (en) * 2013-12-27 2014-04-02 付诚 Beidou satellite-based special article positioning and supervising method and system
CN107533409A (en) * 2015-06-19 2018-01-02 株式会社艾拉博 Coordinate detecting device and coordinate detection method
CN107036770A (en) * 2017-04-18 2017-08-11 浙江理工大学 The leakage detection and localization method of air cooler fin tube bank
CN107765294A (en) * 2017-10-20 2018-03-06 重庆泛嘉地质勘查有限公司 Method of seismic source monitoring
CN110084986A (en) * 2019-04-29 2019-08-02 西人马(厦门)科技有限公司 A kind of circumference safety protection method and device
CN111044975A (en) * 2019-12-10 2020-04-21 北京无线电计量测试研究所 Method and system for positioning earth vibration signal
KR102218181B1 (en) * 2020-11-06 2021-02-23 케이아이티밸리(주) Method for determining real-time alarm for local target earthquake early warning using mechanical learning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于MET振动传感器的二维定位技术研究;王春晖;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20190115(第1期);全文 *

Also Published As

Publication number Publication date
CN113325473A (en) 2021-08-31

Similar Documents

Publication Publication Date Title
US8575929B1 (en) Magnetic anomaly surveillance system using spherical trilateration
US8094833B2 (en) Sound source localization system and sound source localization method
CN105093299A (en) Observation system optimization method based on offset vector tile technology and apparatus thereof
US20090058993A1 (en) Cmos stereo camera for obtaining three-dimensional image
US7690250B2 (en) Method for measuring surface currents using a long-range single station high frequency ground wave radar system
CN107918156B (en) Detect the method and device of subsea node acquisition earthquake data polarity
US10094911B2 (en) Method for tracking a target acoustic source
CN105044669A (en) Three-station time-difference-measuring stereoscopic positioning method
CN113009559B (en) Earthquake evaluation method based on multi-type geophysical data
CN105866636B (en) Transformer substation positioning method based on time difference positioning
CN111399066B (en) Method for processing scalar magnetic anomaly gradient signal based on orthogonal basis function
Raab Quasi-static magnetic-field technique for determining position and orientation
WO2018131946A1 (en) Method and apparatus for performing drive test in mobile communication system
CN113325473B (en) Vibration region determination method and system
CN113960532B (en) Microseism positioning method based on secondary positioning calculation of virtual source
JP3127945B2 (en) Moving sound source identification method and apparatus
JP7030496B2 (en) Underground cavity detection device and underground cavity detection method
CN112525407B (en) Shock wave single multi-sensor array testing method and device
CN113325465B (en) Method and system for quickly positioning ground intrusion target
CN102435866A (en) Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar
US9880008B2 (en) Indoor positioning method and indoor positioning system
CN112857699B (en) Long oil and gas pipeline leakage detection and positioning method of acoustic vector sensor
JP4499617B2 (en) Lightning position limiting system and method
KR100371793B1 (en) A Method for Determining the Location of Sound Source
Cai et al. KCF‐Based Identification Approach for Vibration Displacement of Double‐Column Bents under Various Earthquakes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant