阅读说明：本技术 一种基于平均值约束的层析反演方法 (Average value constraint-based chromatographic inversion method ) 是由 郭昌放 杨真 于 2020-07-07 设计创作，主要内容包括：本发明公开了一种基于平均值约束的层析反演方法,先采用无线电波坑透获取工作面的原始数据；将工作面进行网格离散化,每个网格作为一个像素,利用直射线理论获得每条射线在工作面内所有网格的截距；根据原始无线电波坑透数据、网格截距以及直射线理论,构建无线电波层析反演矩阵方程组；通过已知的反演方法结合反演矩阵方程组,求得每个像素的反演吸收系数值；最后将已揭露的多源先验数据作为平均值约束的数据引入到层析反演计算过程中,不断修正得出的每个像素的反演吸收系数值,最终将修正后的数据通过图像显示,从而得出工作面内的异常分布情况。因此本发明能提高反演得出的工作面内异常分布情况的准确性。(The invention discloses a chromatographic inversion method based on average value constraint, which comprises the steps of firstly, acquiring original data of a working face by adopting a radio wave pit; carrying out grid discretization on the working surface, taking each grid as a pixel, and obtaining the intercept of each ray in all the grids in the working surface by using a straight ray theory; constructing a radio wave chromatography inversion matrix equation set according to original radio wave pit penetration data, grid intercept and a straight ray theory; solving the inversion absorption coefficient value of each pixel by combining an inversion matrix equation set through a known inversion method; and finally, introducing the revealed multi-source prior data serving as data constrained by an average value into a chromatographic inversion calculation process, continuously correcting the obtained inversion absorption coefficient value of each pixel, and finally displaying the corrected data through an image so as to obtain the abnormal distribution condition in the working surface. Therefore, the method can improve the accuracy of the abnormal distribution condition in the working surface obtained by inversion.)

1. A chromatographic inversion method based on mean value constraint is characterized by comprising the following specific steps:

acquiring original data of a working surface by adopting a radio wave pit;

II, carrying out grid discretization on the working surface, and obtaining the intercept of all grids of each ray in the working surface by using a straight ray theory; wherein the number of grids is B, and the number of rays is A; each grid is called a pixel and is represented by x_jThe true absorption coefficient of the jth grid, j ═ 1,2,3 … B;

suppose any ray y_iThe ray has an intercept d in the j-th grid_ij(ii) a If the number of the set rays is A, the ith ray path comprises the following components: 1,2,3 … a;

due to the fact that:

H_i＝H₀-20y_ilg e-20lg r_i(2)

equations (1) and (2) are combined, namely:

in the above formula:

H_iacquired for the ith rayRadio wave field strength data;

H₀is the initial field intensity value of the emission of radio waves;

r_iis the length of the ith ray;

III, constructing a radio wave chromatography inversion matrix equation set according to original radio wave pit penetration data, grid intercept and a straight ray theory;

since the whole process of the coal mine radio wave pit penetration detection is to generate radio waves by using a plurality of transmitting points, wherein the radio waves generated by each transmitting point are excited by a plurality of receiving points to acquire data, the formula (1) can be further generalized as follows:

after simplification, we obtain:

[D][X]＝[Y]

in the formula:

[D] the coefficient matrix is an A x B order coefficient matrix, namely the intercept of each ray passing through each grid, wherein the number of the rays is A, and the number of the grids is B;

[ X ] is a B1-order unknown coefficient matrix, which is a radio wave absorption coefficient of each pixel;

[ Y ] is an A x 1 order constant matrix, and constants related to the actually measured field intensity in each observation mode are known data;

wherein, [ X ]]Obtaining the value of the inverted absorption coefficient x 'of each pixel by known inversion method according to the formula'_j；

IV, in the process of inverting the absorption coefficient of the propagation medium in the working face through the radio wave penetration technology, inverting the absorption coefficient value x 'of each obtained grid'_jAnd (4) carrying out constraint processing:

x 'in the area where no geological abnormality is revealed in the coal mine working face'_jAverage absorption coefficient value x of all grids obtained by current inversion_avgX'_jAnd does not change; x'_jGreater than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

② if x 'in the area where the geological abnormality exists in the coal mine working face'_jThe average absorption coefficient value x of all grids obtained by the current inversion is larger than or equal to_avgX'_jDo not change; x'_jLess than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

The concrete formula is as follows:

wherein:

x'_jthe value of the inversion absorption coefficient corresponding to the jth grid of the working surface;

x_avgtaking the average value of all grid absorption coefficients of the current inversion;

Q₁is larger than x in the area where no geological abnormal area is revealed on the coal mine working face_avgA set of grid sequence numbers of;

Q₂less than x in the area where geological abnormality is revealed on the coal mine working face_avgA set of grid sequence numbers of;

and processing the absorption coefficient value of each pixel obtained by inversion through the average value constraint, and finally displaying the processed data through an image so as to obtain the abnormal distribution condition in the working surface.

Technical Field

The invention relates to a chromatographic inversion method based on average value constraint, and belongs to the technical field of structural inversion in a coal mine working face.

Background

Coal is the main energy and basic industry of China, and geological abnormalities such as faults and collapse columns in the coal mine working face directly affect working face extraction scheme formulation, extraction efficiency and life safety of underground workers. The conditions encountered will become more complex as the depth of coal mining increases in the future. Therefore, the method has very important practical significance for the safe and efficient production of the coal mine by accurately inverting and predicting the geological abnormal structure in the working face.

As a non-contact detection method, the underground radio wave pit penetration technology has the advantages of convenience and quickness in construction, high efficiency, strong resolving power, low cost, no destructiveness and the like, so that the underground radio wave pit penetration technology becomes a main method for detecting and inverting the geological anomaly in the coal mine working face at home and abroad at present. In the process of analyzing and explaining the geological anomaly inversion result in the working face by the radio wave penetration technology, the comprehensive judgment is often carried out by combining multi-source data such as roadway exploration, drilling and even recovery disclosure. However, these published multisource prior data are only used for qualitative analysis of auxiliary results, and do not directly and quantitatively participate in the calculation process of the chromatographic inversion. Not only has certain influence on the accuracy of the prediction result, but also generates certain data resource waste.

Disclosure of Invention

Aiming at the problems in the prior art, the invention provides a chromatographic inversion method based on average value constraint, which introduces the revealed multi-source prior data as average value constraint into the chromatographic inversion calculation process, and continuously corrects the inversion calculation process, thereby improving the accuracy of the abnormal distribution condition in the working surface obtained by inversion.

In order to achieve the purpose, the invention adopts the technical scheme that: a chromatographic inversion method based on average value constraint comprises the following specific steps:

acquiring original data of a working surface by adopting a radio wave pit;

II, carrying out grid discretization on the working surface, and obtaining the intercept of all grids of each ray in the working surface by using a straight ray theory; wherein the number of grids is B, and the number of rays is A; each grid is called a pixel and is represented by x_jThe true absorption coefficient of the jth grid, j ═ 1,2,3 … B;

suppose any ray y_iThe ray has an intercept d in the j-th grid_ij. Wherein if the number of the set rays is A, the ith ray path isThe method comprises the following steps: 1,2,3 … a;

due to the fact that:

H_i＝H₀-20y_ilge-20lgr_i(2)

equations (1) and (2) are combined, namely:

in the above formula:

H_ithe field intensity data of the radio wave collected aiming at the ith ray;

H₀is the initial field intensity value of the emission of radio waves;

r_ithe length of the ith ray is a known value obtained by initial calculation;

III, constructing a radio wave chromatography inversion matrix equation set according to original radio wave pit penetration data, grid intercept and a straight ray theory;

since the whole process of the coal mine radio wave pit penetration detection is to generate radio waves by using a plurality of transmitting points, wherein the radio waves generated by each transmitting point are excited by a plurality of receiving points to acquire data, the formula (1) can be further generalized as follows:

after simplification, we obtain:

[D][X]＝[Y]

in the formula:

[D] the coefficient matrix is an A x B order coefficient matrix, namely the intercept of each ray passing through each grid, wherein the number of the rays is A, and the number of the grids is B;

[ X ] is a B1-order unknown coefficient matrix, which is a radio wave absorption coefficient of each pixel;

[ Y ] is an A x 1 order constant matrix, and constants related to the actually measured field intensity in each observation mode are known data;

wherein, [ X ]]Obtaining the value of the inverted absorption coefficient x 'of each pixel by known inversion method according to the formula'_j；

IV, in the process of inverting the absorption coefficient of the propagation medium in the working face through the radio wave penetration technology, inverting the absorption coefficient value x 'of each obtained grid'_jAnd (4) carrying out constraint processing:

x 'in the area where no geological abnormality is revealed in the coal mine working face'_jAverage absorption coefficient value x of all grids obtained by current inversion_avgX'_jAnd does not change; x'_jGreater than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

② if x 'in the area where the geological abnormality exists in the coal mine working face'_jThe average absorption coefficient value x of all grids obtained by the current inversion is larger than or equal to_avgX'_jDo not change; x'_jLess than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

The concrete formula is as follows:

wherein:

x'_jthe value of the inversion absorption coefficient corresponding to the jth grid of the working surface;

x_avgtaking the average value of all grid absorption coefficients of the current inversion;

Q₁is larger than x in the area where no geological abnormal area is revealed on the coal mine working face_avgA set of grid sequence numbers of;

Q₂less than x in the area where geological abnormality is revealed on the coal mine working face_avgA set of grid sequence numbers of;

and processing the absorption coefficient value of each pixel obtained by inversion through the average value constraint, and finally displaying the processed data through an image so as to obtain the abnormal distribution condition in the working surface.

Compared with the prior art, the method adopts the radio wave pit to obtain the original data of the working surface; carrying out grid discretization on the working surface, taking each grid as a pixel, and obtaining the intercept of each ray in all the grids in the working surface by using a straight ray theory; constructing a radio wave chromatography inversion matrix equation set according to original radio wave pit penetration data, grid intercept and a straight ray theory; solving the inversion absorption coefficient value of each pixel by combining an inversion matrix equation set through a known inversion method; and finally, introducing the revealed multi-source prior data serving as data constrained by an average value into a chromatographic inversion calculation process, continuously correcting the obtained inversion absorption coefficient value of each pixel, and finally displaying the corrected data through an image so as to obtain the abnormal distribution condition in the working surface. Therefore, the method can improve the accuracy of the abnormal distribution condition in the working surface obtained by inversion.

Drawings

Fig. 1 is a schematic diagram of grid discretization of radio wave tomography in the present invention.

Detailed Description

The present invention will be further explained below.

As shown in fig. 1, the method comprises the following specific steps:

acquiring original data of a working surface by adopting a radio wave pit;

II, carrying out grid discretization on the working surface, and obtaining the intercept of all grids of each ray in the working surface by using a straight ray theory; wherein the number of grids is B, and the number of rays is A; each grid is called a pixel and is represented by x_jThe true absorption coefficient of the jth grid, j ═ 1,2,3 … B;

suppose any ray y_iHair fromThe ray end passes through the working face to reach the receiving end, and the intercept of the ray at the j grid is d_ij. If the number of the set rays is A, the ith ray path comprises the following components: 1,2,3 … a;

due to the fact that:

H_i＝H₀-20y_ilge-20lgr_i(2)

equations (1) and (2) are combined, namely:

in the above formula:

H_ithe field intensity data of the radio wave collected aiming at the ith ray;

H₀is the initial field intensity value of the emission of radio waves;

r_ithe length of the ith ray is a known value obtained by initial calculation;

III, constructing a radio wave chromatography inversion matrix equation set according to original radio wave pit penetration data, grid intercept and a straight ray theory;

since the whole process of the coal mine radio wave pit penetration detection is to generate radio waves by using a plurality of transmitting points, wherein the radio waves generated by each transmitting point are excited by a plurality of receiving points to acquire data, the formula (1) can be further generalized as follows:

after simplification, we obtain:

[D][X]＝[Y]

in the formula:

[D] the coefficient matrix is an A x B order coefficient matrix, namely the intercept of each ray passing through each grid, wherein the number of the rays is A, and the number of the grids is B;

[ X ] is a B1-order unknown coefficient matrix, which is a radio wave absorption coefficient of each pixel;

[ Y ] is an A x 1 order constant matrix, and constants related to the actually measured field intensity in each observation mode are known data;

wherein, [ X ]]Obtaining the value of the inverted absorption coefficient x 'of each pixel by known inversion method according to the formula'_j；

IV, in the process of inverting the absorption coefficient of the propagation medium in the working face through the radio wave penetration technology, inverting the absorption coefficient value x 'of each obtained grid'_jAnd (4) carrying out constraint processing:

x 'in the area where no geological abnormality is revealed in the coal mine working face'_jAverage absorption coefficient value x of all grids obtained by current inversion_avgX'_jAnd does not change; x'_jGreater than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

② if x 'in the area where the geological abnormality exists in the coal mine working face'_jThe average absorption coefficient value x of all grids obtained by the current inversion is larger than or equal to_avgX'_jDo not change; x'_jLess than x_avgUsing the average absorption coefficient value x_avgReplacing the absorption coefficient value x 'of the grid'_j；

The concrete formula is as follows:

wherein:

x'_jthe value of the inversion absorption coefficient corresponding to the jth grid of the working surface;

x_avgtaking the average value of all grid absorption coefficients of the current inversion;

Q₁is larger than x in the area where no geological abnormal area is revealed on the coal mine working face_avgA set of grid sequence numbers of;

Q₂less than x in the area where geological abnormality is revealed on the coal mine working face_avgA set of grid sequence numbers of;

and processing the absorption coefficient value of each pixel obtained by inversion through the average value constraint, and finally displaying the processed data through an image so as to obtain the abnormal distribution condition in the working surface.