CN107272061B - A kind of automatic correlation method of infrasound signals and seismic events - Google Patents

A kind of automatic correlation method of infrasound signals and seismic events Download PDF

Info

Publication number
CN107272061B
CN107272061B CN201710513724.XA CN201710513724A CN107272061B CN 107272061 B CN107272061 B CN 107272061B CN 201710513724 A CN201710513724 A CN 201710513724A CN 107272061 B CN107272061 B CN 107272061B
Authority
CN
China
Prior art keywords
seismic events
infrasound signals
signals
event
infrasonic sound
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
CN201710513724.XA
Other languages
Chinese (zh)
Other versions
CN107272061A (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.)
Ctbt Beijing National Data Center
Original Assignee
Ctbt Beijing National Data Center
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 Ctbt Beijing National Data Center filed Critical Ctbt Beijing National Data Center
Priority to CN201710513724.XA priority Critical patent/CN107272061B/en
Publication of CN107272061A publication Critical patent/CN107272061A/en
Application granted granted Critical
Publication of CN107272061B publication Critical patent/CN107272061B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6161Seismic or acoustic, e.g. land or sea measurements

Abstract

The present invention relates to the automatic correlation methods of a kind of infrasound signals and seismic events, it is first determined infrasonic sound station inventory to be associated;For each infrasonic sound station to be associated, calculate separately its epicentral distance between seismic events to be associated, the station is to the azimuth of seismic events, the retrieval beginning and ending time of infrasound signals is determined according to earthquake Time To Event, the P phase and S phase signals for determining infrasonic sound station record obtain its infrasound signals search beginning and ending time;Infrasound signals in the azimuth coverage of each infrasonic sound station are scanned for, qualified infrasound signals is found out and is associated with as the infrasound signals induced by seismic events, and with seismic events, acoustic shock is formed and combines correlating event.The comprehensive multi-signal interconnection constraint condition of the present invention and event region historical information, can identify the infrasound signals induced by seismic events, improve event analysis accuracy;It can also will determine that part excludes in infrasonic monitoring data sequence source signal, reduce signal and judge risk by accident.

Description

A kind of automatic correlation method of infrasound signals and seismic events
Technical field
The invention belongs to earthquake monitoring fields, and in particular to a kind of correlating method of earthquake-infrasound signals, it is especially a kind of The correlating method of infrasound signals and seismic events based on earthquake kind subevent.
Background technique
Infrasonic sound is in addition to for atmospheric event monitoring, event quick-fried for region mine, or even at a distance to larger earthquake magnitude Shake event also has certain monitoring capability, not only can recorde the signal that seismic wave directly acts on the station, moreover it is possible to record shake The source infrasound signal that nearby vibrations such as high mountain, island generate.By comprehensive analysis seismic signal and infrasound signals, it can be achieved that compared with Single means it is more abundant analysis as a result, be more advantageous to event property analysis and differentiate.Due to atmosphere and meteorological condition Complexity, the transmission speed of infrasound signals can not accurately measure as seismic wave, be used for the associated lattice point of seismic network The traditional technologies such as search method can not be suitable for infrasound signals and be associated with.The potential infrasonic sound station relevant to specific seismic events is really Determine, the search matching of phase infrasound signals in unison is not the technological difficulties for needing to solve in infrasonic sound-earthquake association.
Summary of the invention
The object of the present invention is to provide the automatic correlation method of a kind of infrasound signals and seismic events, solve extensive time Sound platform net monitor infrasound signals and global earthquakes event can not the associated problem of Auto-matching, signal erroneous judgement wind can also be reduced Danger;It can be automatically processed in seismic data or on the basis of man-machine interactively is analyzed, realize not phase infrasound signals and seismic events in unison Efficient association, the Wave data Conjoint Analysis for two kinds of technological means of seismic monitoring and infrasonic monitoring provides technical support.
The technical solution of the present invention is as follows: the automatic correlation method of a kind of infrasound signals and seismic events, it is characterised in that: packet Include following steps:
S1, it determines that the corresponding infrasound signals of seismic events to be associated can monitor epicentral distance mag_th: extracting first specifically Body wave magnitude mb, the local magnitude M_L of shake event choose maximum value m therein as the station and screen earthquake magnitude, determine specific earthquake The relationship that the corresponding infrasound signals of event can monitor epicentral distance mag_th and m is as follows:
According to above-mentioned relation, determines and meet earthquake magnitude-epicentral distance being associated with monitoring data in infrasonic monitoring platform net Infrasonic sound station list 1;
S2, it is searched for from history infrasonic sound event database with seismic events distance to be associated less than specific epicentral distance threshold value With reference to event, determine that once the associated infrasonic sound station, formation can be associated with infrasonic sound station list 2 to this area's history time sound events, it is comprehensive Can be associated with infrasonic sound station list 1 with can be associated with infrasonic sound station list 2, form infrasonic sound station inventory to be associated;
S3, for each infrasonic sound station to be associated, calculate separately its epicentral distance between seismic events to be associated disti,j, the azimuth az of the station to seismic eventsi,j;Wherein i is seismic events serial number to be associated, and j is the infrasonic sound station to be associated Serial number;
Utilize epicentral distance disti,jCalculate signal frequency upper limit fmaxi,j, the signal theory under different equivalent sonic conditions walks When, choose walking time error delt_timei,j
S4, association infrasonic sound station aspect angle error;
S5, the retrieval beginning and ending time t_start that infrasound signals are determined according to earthquake Time To Eventi,jAnd t_endi,j
S6, the P phase and S phase signals for determining infrasonic sound station record obtain its infrasound signals search beginning and ending time;
S7, the infrasound signals in the azimuth coverage of each infrasonic sound station are scanned for, finds out qualified infrasonic sound letter It number is associated with as the infrasound signals induced by seismic events, and with seismic events, forms acoustic shock and combine correlating event.
The utility model has the advantages that on the one hand the present invention can be known by establishing the correlating methods of a kind of infrasound signals and seismic events The infrasound signals that Chu not be induced by seismic events, improve event analysis accuracy;It on the other hand can be by infrasonic monitoring data sequence It determines that part excludes in source signal, reduces signal and judge risk by accident.The invention proposes infrasound signals can monitoring distance with earthquake The variation relation of earthquake magnitude gives infrasound signals frequency content with the changing rule of epicentral distance, it is determined that phase signals are not examined in unison Infrasonic sound inspection may be implemented by the comprehensive utilization of a variety of constraint conditions and event region historical information in rope time interval Survey signal is associated with earthquake kind subevent.The present invention can be not only used for infrasonic sound-seismic events association of man-machine interactively auditing events, Also it can be used for infrasonic sound-seismic events association that computer automatically processes event.By using event correlation station analysis on its rationality, The false detection seismic event identification technology such as event correlation station type analysis rejects the false association thing for automatically processing generation Part, to improve the associated accuracy of infrasonic sound-seismic events.It is the infrasonic sound-for driving seed that the present invention, which is realized with seismic events, Seismic events correlating method can be not only used for infrasonic sound-earthquake association of earthquake event, it can also be used to time of ground burst event Sound-earthquake association.It solves extensive infrasonic sound platform net monitoring infrasound signals and is associated with problem with the Auto-matching of global earthquakes event, Auxiliary waveform analysis personnel's progress infrasound signals are quick, accurate correlation is analyzed;The present invention is also multiple waveforms monitoring means simultaneously The beneficial trial of Conjoint Analysis seismic events.
Detailed description of the invention
Fig. 1 is that infrasound signals association flow chart is merged in acoustic shock of the invention.
Specific embodiment
To keep the purpose of the present invention, content and advantage clearer, below with reference to embodiment, to specific implementation of the invention Mode is described in further detail.
The present invention provides the automatic correlation method of a kind of infrasound signals and seismic events, it is characterised in that: including following step It is rapid:
If the seismic events that data processing system automatically generates S1, are carried out with infrasonic sound-earthquake association, carry out first automatic Handle false event identification, the method for use include: event correlation station analysis on its rationality, event correlation station type analysis, Event correlation signal analysis etc.;
If audit seismic events are carried out with infrasonic sound-earthquake association directly since step 2;
S2, it determines that the corresponding infrasound signals of specific seismic events can monitor epicentral distance mag_th: extracting specific earthquake first Body wave magnitude mb, the local magnitude M_L of event choose maximum value MAX (mb, ML) therein as the station and screen earthquake magnitude m, determine It is as follows with the relationship of m that the corresponding infrasound signals of specific seismic events can monitor epicentral distance mag_th (unit is radian):
According to the relationship of above-mentioned epicentral distance and m, determine that meeting earthquake magnitude-epicentral distance in infrasonic monitoring platform net has monitoring number According to be associated with infrasonic sound station list 1;
S3, search and specific earthquake incident distance to be associated are less than specific epicentral distance threshold from history infrasonic sound event database The reference event of value determines that once the associated infrasonic sound station, formation can be associated with infrasonic sound station list 2 to this area's history time sound events, Synthesis can be associated with infrasonic sound station list 1 with can be associated with infrasonic sound station list 2, form infrasonic sound station inventory to be associated;Specific earthquake centre It is preset according to demand away from threshold value, such as 100km;
S4, due to seismic events position to be associated it is known that for each infrasonic sound station to be associated, calculate separately its with to Epicentral distance dist between correlating eventi,j, the azimuth az of the station to seismic eventsi,j;Wherein i is event serial number to be associated, and j is Infrasonic sound station serial number to be associated;
Utilize disti,jCalculate the signal frequency upper limit, when the signal theory under different equivalent sonic conditions is walked and selection is walked When error:
f maxi,j=4-0.055*disti,j (2)
t1i,j=111.1949*disti,j/vh (3)
t2i,j=111.1949*disti,j/vm (4)
t3i,j=111.1949*disti,j/vl (5)
Wherein vh、vm、vlRespectively indicate the upper limit, average value and the lower limit of velocity equivalent;delt_time,dealt_time1, Dealt_time2 is respectively that infrasound signals walking time error threshold value, the low error coefficient of walking time error identification and walking time error identification are high Error coefficient;Dealt_time1, dealt_time2 can be preset according to demand;
S5, association infrasonic sound station aspect angle error: as infrasonic sound station epicentral distance disti,jGreater than preset base When quasi- epicentral distance, amplify azimuth angle reference error:
dealt_azi,j'=dealt_azi,j*dealt_az_times (7)
Wherein dealt_azi,j' it is amplified azimuth angle reference error, dealt_azi,jFor azimuth angle reference error, Dealt_az is azimuth angle reference error, and dealt_az_times is the amplification factor of benchmark error;Consider infrasound signals long distance It is influenced from transmission by atmosphere, presets certain benchmark epicentral distance, the benchmark as setting azimuth angle error;
Otherwise, without being enlarged;
S6, the retrieval beginning and ending time t_start for determining infrasound signalsi,jAnd t_endi,j
Initial time is retrieved for infrasound signals, if t1i,j-0.5*dealt_timei,j> 0:
t_starti,j=t1i,j+toridi-0.5*dealt_timei,j (8)
Otherwise:
t_starti,j=toridi (9)
Infrasound signals retrieve the end time are as follows:
t_endi,j=toridi+t3i,j*dealt_timei,j (10)
Wherein toridiFor Time To Event.
S7, the P phase and S phase signals for determining infrasonic sound station record, signal search beginning and ending time are respectively as follows:
tmin_pi,j=disti,j·pslow_min-ptime_delta+toridi (11)
tmax_pi,j=disti,j·pslow_max+ptime_delta+toridi (12)
tmin_Si,j=disti,j·Sslow_min+Stime_delta+toridi (13)
tmax_Si,j=disti,j·Sslow_max+Stime_delta+toridi (14)
Wherein, pslow_min is P phase signals slowness minimum value, pslow_max is P phase signals slowness maximum value, ptime_ Delta is P phase signals walking time error limit value, Sslow_min is S phase signals slowness minimum value, Sslow_max is that S phase signals are slow Spend maximum value, Stime_delta is S phase signals walking time error limit value;
S8, the infrasound signals search time section according to above-mentioned formula determination, signal frequency range, to each infrasonic sound station Azimuth [azi,j-dealt_azi,j′,azi,j+dealt_azi,j'] the detection infrasound signals in range scan for, find out symbol The infrasound signals of conjunction condition are associated with as the infrasound signals induced by seismic events, and with seismic events, are formed acoustic shock joint and are closed Connection event, acoustic shock fusion infrasound signals association process can refer to Fig. 1.
S9, using associated earthquake and infrasound signals to seismic events into analyzing again in detail, improve event analysis accuracy;
S10, it can will be excluded in infrasonic monitoring data sequence source signal by the infrasound signals that seismic events induce, reduce letter Number erroneous judgement risk.
The present invention retrieves matching associated infrasound signals using seismic events as seed driving source.Due to computer System detects seismic events automatically, and there may be more false events, if therefore when automatically processing event as kind of subevent, False event screening identification is carried out to it first, the associated efficiency of acoustic shock on the one hand can be improved, reduce infrasound signals and falseness The unnecessary association of seismic events;On the other hand it can prevent effective infrasound signals from " being accounted for due to false earthquake event correlation Have ", cause it that cannot be associated with real event source.
For any one earthquake kind to be associated subevent, first choice, which needs to differentiate, determines infrasonic sound station inventory to be associated;It is first First according to the intensity of seismic events, judge infrasound signals monitors epicentral distance, and event within the scope of epicentral distance can be monitored by finding out The infrasonic sound station that front and back there are monitoring data occurs, forms inventory to be associated.The relationship that seismic events earthquake magnitude-can monitor epicentral distance can It is statisticallyd analyze and is obtained according to historical events.Furthermore it is possible to the infrasound signals with reference to event region historical events are associated with situation, Such as historical events within the scope of the radius 100km of event epicenter periphery, the associated infrasonic sound station of historical events is added to be associated clear Dan Zhong.
According to the real information of event in secondary acoustic reference event (IRED), the equivalent velocity of sound of infrasound signals overland propagation is determined Valid interval, minimum value, average value, maximum value including the equivalent velocity of sound can avoid infrasound signals by using the equivalent velocity of sound Complicated refraction and reflection process in communication process simplify estimation process when infrasound signals are walked.In addition, according to history reference event Detection and analysis result, the relative positional relationship of the infrasonic sound station and event of infrasound signals, determine infrasound signals frequency content with biography The approximate attenuation relation of defeated distance, to avoid the unreasonable association of infrasound signals, such as remote event correlation high frequency infrasonic sound is believed Number.
Then for the infrasonic sound station in inventory to be associated, the orientation, distance between event are all it has been determined that can root According to seismic wave P wave, S wave transmission speed section in the earth's crust, the velocity equivalent section of I phase infrasound signals can determine three that earthquake generates The section when signal of kind sound phase (IPx, ISx, I) is walked.In conjunction with given signal walking time error and azimuth angle error, from infrasound signals Search matches qualified infrasound signals in testing result, and it is associated with earthquake kind subevent, to be formed Correlating event is merged in infrasonic sound-earthquake.
The infrasonic sound station includes but is not limited to the built international monitoring system infrasonic sound platform of Comprehensive Nuclear Test-Ban Treaty Organization It stands, refers to all kinds of installation micromanometers, micro-microphone etc. for recording the monitoring station of atmospheric pressure disturbance;The seismograph station Including but not limited to international monitoring system seismograph station refers to all kinds of installation displacement meters, speedometer and accelerometer and for surveying The monitoring station of amount ground motion conditions;Earthquake kind subevent refers to the seismic events for being associated with and being formed by seismic network;It is described The event of automatically processing refer to by computer system through signal detection, signature analysis, platform gateway connection and etc. the earthquake automatically generated Event;The magnitude is to characterize the measurement, including Near Shock Magnitude, surface wave magnitude, body wave magnitude etc. of earthquake power;The shake In away from referring to the distance of seismograph station or the infrasonic sound station to earthquake epicenter;The infrasonic sound P wave sound mutually refers to infrasonic sound station record Infrasound signals with seismic wave P wave characteristic;The S wave sound mutually refer to infrasonic sound station record with seismic wave S wave characteristic Infrasound signals;The I wave sound mutually refers to the infrasound signals with acoustic characteristic of infrasonic sound station record.It is characterized in that
The epicentral distance that monitors refers within the scope of a certain earthquake magnitude, can monitor the maximum distance of infrasound signals;
The infrasound signals azimuth angle error refers to the azimuth characteristic parameter and the infrasonic sound station to thing of detection infrasound signals The deviation of part earthquake centre angle;
The secondary acoustic reference event refers to confirmed and true history time sound events, including event information and associated station Information etc.;
The false event recognition methods refers to the method for automatically processing generated accidentally correlating event for identification, including But it is not limited to document " automatic detection seismic events STUDY ON SCREENING " (seismomagnetic observation and research, volume 34 of in August, 2013 3/4 phase) listed by method;
The epicentral distance frequency relation refers to the specific symbol between the centre frequency and station epicentral distance of detection infrasound signals Conjunction relationship;
The equivalent velocity of sound, which refers to, assumes the infrasound signals speed that cambered surface is transmitted at the earth's surface;
It is described be associated with infrasonic sound station inventory and refer to may all monitor earthquake infrasound signals and in event time period There is the infrasonic sound station list of data;
Embodiment:
Using infrasonic sound-seismic events correlating method of above-mentioned foundation, test is carried out to 7.6 grades of earthquakes of the Philippines and is said Bright, table 1 is that event International Data Center (IDC) waveform expert audits infrasound signals associated by bulletin, this result can be used as judgement A kind of one of the judgment criteria of earthquake infrasonic sound correlating method efficiency based on earthquake kind subevent of the present invention.Prohibited using comprehensive Only nuclear test International Data Center (IDC) of treaty organization carries out acoustic shock association analysis, association results such as table based on automatically processing event Shown in 2, the method for the invention is correctly associated with tri- kinds of signals of P, S, I, and the associated station includes I39PW, I30JP, I45RU, Due to epicentral distance range setting in algorithm, fail to be associated with the I59US station compared with REB.
1 International Data Center (IDC) of table audits the associated infrasound signals of bulletin
Table 2 is based on Philippine's seismic events infrasound signals auto-associating result of the invention
It can be seen that in conjunction with table 1 and table 2, can relatively accurately carry out infrasonic sound-seismic events using the present invention and be associated with, matching All suspicious infrasound signals are associated with out, after automatically processing modified result, can be obtained accurate event result.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations Also it should be regarded as protection scope of the present invention.

Claims (7)

1. the automatic correlation method of a kind of infrasound signals and seismic events, it is characterised in that: the following steps are included:
S1, it determines that the corresponding infrasound signals of seismic events to be associated can monitor epicentral distance mag_th: extracting earthquake to be associated first Body wave magnitude mb, the local magnitude M_L of event choose maximum value m therein as the station and screen earthquake magnitude, determine earthquake to be associated The relationship that the corresponding infrasound signals of event can monitor epicentral distance mag_th and m is as follows:
According to above-mentioned relation, determine that meeting earthquake magnitude-epicentral distance in infrasonic monitoring platform net is associated with infrasonic sound with monitoring data Station list 1;
S2, the reference for being less than specific epicentral distance threshold value with seismic events distance to be associated is searched for from history infrasonic sound event database Event determines that once the associated infrasonic sound station, formation can be associated with infrasonic sound station list 2 to this area's history time sound events, and synthesis can close Connection infrasonic sound station list 1 with can be associated with infrasonic sound station list 2, form infrasonic sound station inventory to be associated;
S3, for each infrasonic sound station to be associated, calculate separately its epicentral distance dist between seismic events to be associatedi,j, platform It stands to the azimuth az of seismic eventsi,j;Wherein i is seismic events serial number to be associated, and j is infrasonic sound station serial number to be associated;
Utilize epicentral distance disti,jCalculate signal frequency upper limit f maxi,j, the signal theory under different equivalent sonic conditions is when walking, Choose walking time error delt_timei,j
fmaxi,j=4-0.055*disti,j
t1i,j=111.1949*disti,j/vh
t2i,j=111.1949*disti,j/vm
t3i,j=111.1949*disti,j/vl
Wherein vh、vm、vlRespectively indicate the upper limit, average value and the lower limit of the equivalent velocity of sound;Dealt_time1, dealt_time2 points It Wei not the low error coefficient of infrasound signals walking time error identification and the high error coefficient of walking time error identification;dealt_time1,dealt_ Time2 is preset according to demand;
S4, association infrasonic sound station aspect angle error;
S5, the retrieval beginning and ending time t_start that infrasound signals are determined according to earthquake Time To Eventi,jAnd t_endi,j
Initial time is retrieved for infrasound signals, if t1i,j-0.5*dealt_timei,j> 0:
t_starti,j=t1i,j+toridi-0.5*dealt_timei,j
Otherwise:
t_starti,j=toridi
Infrasound signals retrieve the end time are as follows:
t_endi,j=toridi+t3i,j*dealt_timei,j
Wherein toridiFor seismic events time of origin;
S6, the P phase and S phase signals for determining infrasonic sound station record obtain its infrasound signals search beginning and ending time:
tmin_pi,j=disti,j·pslow_min-ptime_delta+toridi
tmax_pi,j=disti,j·pslow_max+ptime_delta+toridi
tmin_Si,j=disti,j·Sslow_min+Stime_delta+toridi
tmax_Si,j=disti,j·Sslow_max+Stime_delta+toridi
Wherein, pslow_min is P phase signals slowness minimum value, pslow_max is P phase signals slowness maximum value, ptime_ Delta is P phase signals walking time error limit value, Sslow_min is S phase signals slowness minimum value, Sslow_max is that S phase signals are slow Spend maximum value, Stime_delta is S phase signals walking time error limit value;
S7, the infrasound signals in the azimuth coverage of each infrasonic sound station are scanned for, finds out qualified infrasound signals and makees It for the infrasound signals induced by seismic events, and is associated with seismic events, forms acoustic shock and combine correlating event.
2. the automatic correlation method of a kind of infrasound signals and seismic events according to claim 1, it is characterised in that: if right The seismic events that data processing system automatically generates carry out infrasonic sound-earthquake association, then carry out automatically processing false event knowledge first S1 is not entered step afterwards;If carrying out infrasonic sound-earthquake to audit seismic events to be associated with, it is directly entered step S1.
3. the automatic correlation method of a kind of infrasound signals and seismic events according to claim 1, it is characterised in that: S2 In, specific epicentral distance threshold value is preset according to demand.
4. the automatic correlation method of a kind of infrasound signals and seismic events according to claim 1, it is characterised in that: S4 In, as infrasonic sound station epicentral distance disti,jWhen greater than preset benchmark epicentral distance, amplify azimuth angle reference error, otherwise, Without being enlarged;
dealt_azi,j'=dealt_azi,j*dealt_az_times
Wherein dealt_azi,j' it is amplified azimuth angle reference error, dealt_azi,jFor azimuth angle reference error, dealt_ Az_times is the amplification factor of benchmark error;Consider that infrasound signals are transmitted at a distance is influenced by atmosphere, presets one Determine benchmark epicentral distance, the benchmark as setting azimuth angle error.
5. the automatic correlation method of a kind of infrasound signals and seismic events according to claim 4, it is characterised in that: each time The azimuth coverage of the sound station is [azi,j-dealt_azi,j′,azi,j+dealt_azi,j′]。
6. the automatic correlation method of described in any item a kind of infrasound signals and seismic events according to claim 1~5, feature It is: using associated earthquake and infrasound signals seismic events progress is analyzed again in detail, improve event analysis accuracy.
7. the automatic correlation method of described in any item a kind of infrasound signals and seismic events according to claim 1~5, feature It is: can will be excluded in infrasonic monitoring data sequence source signal by the infrasound signals that seismic events induce, reduce signal erroneous judgement Risk.
CN201710513724.XA 2017-06-29 2017-06-29 A kind of automatic correlation method of infrasound signals and seismic events Active CN107272061B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710513724.XA CN107272061B (en) 2017-06-29 2017-06-29 A kind of automatic correlation method of infrasound signals and seismic events

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710513724.XA CN107272061B (en) 2017-06-29 2017-06-29 A kind of automatic correlation method of infrasound signals and seismic events

Publications (2)

Publication Number Publication Date
CN107272061A CN107272061A (en) 2017-10-20
CN107272061B true CN107272061B (en) 2019-02-05

Family

ID=60071336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710513724.XA Active CN107272061B (en) 2017-06-29 2017-06-29 A kind of automatic correlation method of infrasound signals and seismic events

Country Status (1)

Country Link
CN (1) CN107272061B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108844616B (en) * 2018-04-23 2020-06-30 中国人民解放军火箭军工程大学 Infrasound event detection method and system
CN109633548B (en) * 2018-12-10 2023-03-24 禁核试北京国家数据中心 Gateway connection method for underwater sound platform
CN109669185B (en) * 2018-12-10 2023-03-24 禁核试北京国家数据中心 Infrasound table network directional search correlation method
CN111208556B (en) * 2020-01-14 2022-04-26 禁核试北京国家数据中心 Method for estimating time-in interval of monitoring signals of underwater acoustic station
CN113607268B (en) * 2021-01-26 2024-01-09 禁核试北京国家数据中心 Regional infrasound event automatic association method
CN114721049B (en) * 2022-03-11 2024-04-09 新疆维吾尔自治区地震局 Virtual seismic station azimuth participating off-grid seismic positioning method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996018119A1 (en) * 1994-12-06 1996-06-13 Farnsworth David F Method for forecasting an earthquake from precusor signals
CN105572738A (en) * 2016-01-29 2016-05-11 禁核试北京国家数据中心 Method for detecting nuclear explosion earthquake event at specific area by means of single station
CN105716707A (en) * 2015-12-10 2016-06-29 中国航空工业集团公司北京长城计量测试技术研究所 Ultra-low-frequency abnormal infrasound signal judging method
CN106382981A (en) * 2016-08-26 2017-02-08 中国地震局地壳应力研究所 Single station infrasonic wave signal recognition and extraction method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996018119A1 (en) * 1994-12-06 1996-06-13 Farnsworth David F Method for forecasting an earthquake from precusor signals
CN105716707A (en) * 2015-12-10 2016-06-29 中国航空工业集团公司北京长城计量测试技术研究所 Ultra-low-frequency abnormal infrasound signal judging method
CN105572738A (en) * 2016-01-29 2016-05-11 禁核试北京国家数据中心 Method for detecting nuclear explosion earthquake event at specific area by means of single station
CN106382981A (en) * 2016-08-26 2017-02-08 中国地震局地壳应力研究所 Single station infrasonic wave signal recognition and extraction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
互相关算法在次声监测数据处理中的应用;唐伟 等;《环境工程》;20101231;第28卷(第6期);第83-86、88页

Also Published As

Publication number Publication date
CN107272061A (en) 2017-10-20

Similar Documents

Publication Publication Date Title
CN107272061B (en) A kind of automatic correlation method of infrasound signals and seismic events
CN102495425B (en) Energy-based method for automatically locating earthquake focus of microearthquake
CN105116440B (en) A kind of slope rock mass monitoring system and monitoring method
CN106154332B (en) A kind of borehole microseismic longitudinal and shear wave event first arrival recognition methods
CN107290787B (en) A kind of monitoring signals correlating method of earthquake infrasonic sound with the location station
Arosio et al. Analysis of microseismic signals collected on an unstable rock face in the Italian Prealps
CN104483700B (en) Formation fracture is monitored and early warning system and method
US11789173B1 (en) Real-time microseismic magnitude calculation method and device based on deep learning
CN110297271A (en) A kind of simple component probe P wave first arrival-time modification method for mine shake alarm
Puglia et al. Extensive characterization of Italian accelerometric stations from single-station ambient-vibration measurements
CN105093314B (en) A kind of method for determining microseism focus
Zeng et al. Turning a telecom fiber‐optic cable into an ultradense seismic array for rapid postearthquake response in an urban area
Scafidi et al. Robust picking and accurate location with RSNI‐Picker2: Real‐time automatic monitoring of earthquakes and nontectonic events
Hansen et al. Three recent larger earthquakes offshore Norway
US11906678B2 (en) Seismic observation device, seismic observation method, and recording medium on which seismic observation program is recorded
Takanami et al. The SIL seismological data acquisition system—As operated in Iceland and in Sweden—
CN115016006B (en) Microseism positioning method and positioning system based on mine safety monitoring
CN107884821B (en) A method of estimation microseism direction of arrival and distance
Wen et al. Earthquake early warning technology progress in Taiwan
Münchmeyer et al. Learning the deep and the shallow: Deep‐learning‐based depth phase picking and earthquake depth estimation
Wang et al. Cumulative absolute velocity prediction for earthquake early warning with deep learning
Kumar et al. Site classification of seismic recording stations of Garhwal region of earthquake early warning system for Uttarakhand, India
Stepnov et al. Earthworm-based automatic system for real-time calculation of local earthquake source parameters
CN104656135B (en) The accurate method of discrimination and system of a kind of earthquake shearing wave
Ladina et al. Strong-motion observations recorded in strategic public buildings during the 24 August 2016 Mw 6.0 Amatrice (central Italy) earthquake

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