CN-115775598-B - Seismic anomaly extraction method for multi-source satellite methane column concentration product

CN115775598BCN 115775598 BCN115775598 BCN 115775598BCN-115775598-B

Abstract

The invention discloses an earthquake anomaly extraction method for a multi-source satellite methane column concentration product, which comprises the following steps of S1, selecting a research area and determining applicability of the methane earthquake anomaly information extraction method, S2, preparing a multi-source satellite data methane gas product, S3, producing a multi-source satellite data fused methane column concentration product, S4, calculating an anomaly index of each pixel methane column concentration in the research area based on the fused methane column concentration product, and S5, analyzing an anomaly high value and predicting a dangerous area based on the anomaly index. The invention mainly focuses on the advantages of a multi-source satellite data channel and an orbit, provides a novel multi-source satellite-based methane column concentration seismic anomaly extraction method, solves the problems of low precision of different channel satellite methane column concentration data products and low single satellite coverage, can provide support for remote sensing gas geochemical seismic service application extraction, and is beneficial to promoting domestic satellite products to enter seismic monitoring service.

Inventors

CUI JING
ZOU MINGMIN
XU SONG
CHU WEI
WANG XU
WANG LIN
HUANG YALAN
SHEN XUHUI

Assignees

应急管理部国家自然灾害防治研究院

Dates

Publication Date: 20260508
Application Date: 20221129

Claims (5)

1. The seismic anomaly extraction method for the multi-source satellite methane column concentration product is characterized by comprising the following steps of: S1, selecting a research area and determining the applicability of a methane seismic anomaly information extraction method; s2, preparing a multi-source satellite data methane gas product; s3, producing a methane column concentration product fused with multi-source satellite data; S4, calculating abnormal indexes of the concentration of the methane columns of each pixel in the research area based on the fused methane column concentration product; s5, analyzing an abnormal high value and predicting a dangerous area based on the abnormality index; The specific method of the step S3 comprises the following substeps: s3-1, space-time disk matching multivariate data; s3-2, carrying out methane product convolution fusion based on the load parameters; The specific method of the step S3-2 is that based on space-time matching data, an average kernel function convolution of short wave infrared channel load is adopted to process methane profile products, and after convolution processing, the methane profile products are subjected to convolution processing according to the formula: ; and obtaining the concentration Xch of the methane mixing ratio of the vertical dry air after fusion, wherein, Representing the earth pressure, Representing the atmospheric head pressure.
2. The method for seismic anomaly extraction of multi-source satellite methane column concentration products of claim 1, wherein the specific method of step S1 comprises the sub-steps of: S1-1, selecting a research area; s1-2, carrying out brief analysis on structural topography of a research area; S1-3, preliminarily determining the regional applicability of the abnormality extraction method from the regional geological angle.
3. The method for seismic anomaly extraction of multi-source satellite methane column concentration products of claim 1, wherein the specific method of step S2 comprises the sub-steps of: s2-1, preparing a near-surface methane concentration product inverted by a short-wave infrared channel; S2-2, preparing a methane profile product for thermal infrared channel inversion.
4. The method for seismic anomaly extraction of multi-source satellite methane column concentration products of claim 1, wherein in step S4, according to the formula: ; Obtaining abnormality index Wherein, the method comprises the steps of, Representing the time of observation Corresponding longitude is Latitude of A value of CH 4 content; Representing the time of observation Corresponding longitude is Latitude of Background value at CH 4 , sigma (x, y, t) is the same location (longitude is Latitude of ) At the same time ) Corresponding standard deviation, where i is year and n=year.
5. The method for seismic anomaly extraction of multi-source satellite methane column concentration products of claim 1, wherein the specific method of step S5 comprises the sub-steps of: S5-1, analyzing the space-time evolution characteristics of the concentration change of the methane column in the research area; S5-2, screening a high-value center of the concentration of the methane column related to the construction; s5-3, based on the steps, focusing on the high-value area as a dangerous area.

Description

Seismic anomaly extraction method for multi-source satellite methane column concentration product Technical Field The invention belongs to the technical field of seismic monitoring, and particularly relates to a seismic anomaly extraction method for a multi-source satellite methane column concentration product. Background Internal and external researches show that the major earthquake can cause the change of trace gas components in the atmosphere, such as CH 4、CO2 and the like, and the parameter changes can indicate the intensity of construction activity and have certain relevance with earthquake pregnancy period. At present, the remote sensing monitoring of the earthquake gas at home and abroad has made important progress, but the application analysis is mainly carried out based on foreign data (AIRS, GOSAT, TOMS, MOPIT, etc), and the application of domestic satellite data is less. And satellite methane column concentration data products aiming at a single data source are mainly used, but all have some problems. For example, short wave infrared channels are more sensitive to near-formation atmosphere CH 4, but cannot invert the profile. Thermal infrared channels are highly sensitive to the middle and high layer atmospheres, but not to the bottom layer atmospheres. This results in an undesirable accuracy of the single channel satellite methane column concentration product. The satellite of 28 minutes GF-5 (01) at 5.9.2.2018 successfully transmits in a Taiyuan satellite transmitting center, a near infrared band of 1.65 mu m of a main greenhouse Gas Monitor (GMI) is carried for observing the vertical column concentration of inversion CH 4, and an in-orbit test result shows that the acquired Level1 spectrum data is better in quality than that of an international similar satellite at shoulder. The short wave infrared channel corresponding to the sensor has higher sensitivity to the low-level atmosphere CH 4, can make up for the shortages of some foreign satellite thermal infrared satellites, and ensures that the data quality of methane gas products is more perfect. But its data acquisition mode determines its low coverage. Disclosure of Invention Aiming at the defects in the prior art, the invention provides a seismic anomaly extraction method for a multi-source satellite methane column concentration product, which solves the problems of low precision of different channel satellite methane column concentration data products and low single satellite coverage. In order to achieve the aim of the invention, the invention adopts the following technical scheme: The seismic anomaly extraction method for the multi-source satellite methane column concentration product comprises the following steps: S1, selecting a research area and determining the applicability of a methane seismic anomaly information extraction method; s2, preparing a multi-source satellite data methane gas product; s3, producing a methane column concentration product fused with multi-source satellite data; S4, calculating abnormal indexes of the concentration of the methane columns of each pixel in the research area based on the fused methane column concentration product; S5, analyzing the abnormal high value and predicting the dangerous area based on the abnormality index. Further, the specific method of step S1 comprises the following sub-steps: S1-1, selecting a research area; s1-2, carrying out brief analysis on structural topography of a research area; S1-3, preliminarily determining the regional applicability of the abnormality extraction method from the regional geological angle. Further, the specific method of step S2 comprises the following sub-steps: s2-1, preparing a near-surface methane concentration product inverted by a short-wave infrared channel; S2-2, preparing a methane profile product for thermal infrared channel inversion. Further, the specific method of step S3 comprises the following sub-steps: s3-1, space-time disk matching multivariate data; s3-2, carrying out methane product convolution fusion based on the load parameters. Further, the specific method of the step S3-2 is that based on space-time matching data, an average kernel function convolution processing methane profile product of short wave infrared channel load is adopted, and after convolution processing, the methane profile product is subjected to convolution processing according to the formula: And obtaining the concentration Xch4 of the mixed ratio of the vertical dry air methane, wherein P 0 represents the surface pressure and P T represents the top pressure of the atmosphere. Further, in step S4, according to the formula: ALICE(x,y,t)=[G(x,y,t)-Gref(x,y,t)]/σ(x,y,t) An abnormal index ALICE is obtained by selecting a certain period before and after an earthquake, wherein G i (x, y, t) represents observation time t, the corresponding longitude is x, the latitude is a CH4 content value at y, G ref (x, y, t) represents the observation time t, the corresponding longitu