[go: up one dir, main page]

CN115935120B - Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data - Google Patents

Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data Download PDF

Info

Publication number
CN115935120B
CN115935120B CN202310022850.0A CN202310022850A CN115935120B CN 115935120 B CN115935120 B CN 115935120B CN 202310022850 A CN202310022850 A CN 202310022850A CN 115935120 B CN115935120 B CN 115935120B
Authority
CN
China
Prior art keywords
goaf
model
grid
density
population
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
CN202310022850.0A
Other languages
Chinese (zh)
Other versions
CN115935120A (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.)
XI'AN RESEARCH INSTITUTE OF CHINA COAL RESEARCH INSTITUTE
Beijing Institute of Spacecraft Environment Engineering
Original Assignee
XI'AN RESEARCH INSTITUTE OF CHINA COAL RESEARCH INSTITUTE
Beijing Institute of Spacecraft Environment Engineering
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 XI'AN RESEARCH INSTITUTE OF CHINA COAL RESEARCH INSTITUTE, Beijing Institute of Spacecraft Environment Engineering filed Critical XI'AN RESEARCH INSTITUTE OF CHINA COAL RESEARCH INSTITUTE
Priority to CN202310022850.0A priority Critical patent/CN115935120B/en
Publication of CN115935120A publication Critical patent/CN115935120A/en
Application granted granted Critical
Publication of CN115935120B publication Critical patent/CN115935120B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

The invention discloses a goaf scale measurement calculation method based on muir remote sensing data, which adopts a genetic algorithm solving mode, divides a model of an initial form of a coal mine goaf into grids, respectively calculates adaptation values of individuals through a Monte Carlo program Geant4 in a mode of constructing a model population according to a density normal (1) and a density loss (0), generates a new population through multiple hybridization variation of the individuals, obtains a model solution set meeting conditions through a genetic algorithm, gives out a final model solution and model uncertainty on the basis of the model solution set through a bootstrap algorithm, and finally gives out a most reasonable three-dimensional density structure model and uncertainty of the coal mine goaf by performing inverse coding according to the density normal (1) and the density loss (0). The invention overcomes the limitation of the traditional technology in complex goaf detection, and can give out a three-dimensional density structural model of the goaf of the coal mine and corresponding uncertainty.

Description

Goaf scale measurement calculation method based on muon remote sensing data
Technical Field
The invention belongs to the field of geological exploration treatment, and relates to a detection scale reconstruction method of a coal mine goaf based on muon remote sensing data.
Background
Coal is the most important one-time energy source in China for a long time, and the construction of environment-friendly, benefit-type and safety-type coal mines is the necessary premise of coal mine production, and the problems of safety, environment protection, coal mine benefit and the like are very important to the control of mine geological disasters.
The goaf of the coal mine is a common mine disaster, the old roof at the rear of the coal face collapses in the range of the falling step pitch, the old roof collapses under the action of mine pressure, the goaf presents dynamic continuity change, a large number of gaps, large holes and other void spaces with different scales are generated, and therefore, the grasping of a complex goaf geometric structure is extremely necessary.
In the past, when a goaf flow field is processed, the goaf is assumed to be a continuous medium, and is regarded as a porous medium category based on an REV concept. The method is influenced by various factors such as physical and mechanical properties of the coal and rock, joint, coal mining method, mining height and the like, the shapes of the coal and rock fragments are various, the formed cavity space is changed in a huge variety, and the cavity morphology of a newly-produced goaf of the same stope face is dynamically and randomly changed. It can be said that two identical new goaf void morphologies cannot be found either in time or in space. In addition, due to limitations of test conditions and equipment, no reasonable test method is found at present to know the real internal cavity structure of the goaf.
In the prior art, the goaf range and the residual settlement are evaluated and judged usually through drilling and geophysical prospecting, the method is only suitable for determining goafs in small ranges, and drilling and geophysical prospecting objects are blind and have no pertinence; specifically, the drilling means is to arrange 1 supplementary exploration line along goaf trend, and supplementary exploration drilling holes are constructed on the supplementary exploration line at intervals, and the drilling holes cannot cover the survey line and the region outside the survey line in a high density manner, and may cause goaf collapse.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention aims to provide a goaf scale measurement calculation method based on muon remote sensing data, which is characterized in that a cosmic ray receiving device is utilized to detect the change of a cosmic ray when the cosmic ray passes through holes with different density lengths, and a three-dimensional density inversion algorithm is adopted to detect the density abnormality of the goaf, so that the form of the internal density structure of the goaf is given.
In order to solve the technical problems, the invention adopts the following technical scheme:
a goaf scale measurement calculation method based on muon remote sensing data comprises the following steps:
step 1, laying a muon detector below a coal mine goaf model, and recording the incident position coordinates (x in ,y in ,z in ) And the exit position coordinates (x out ,y out ,z out ) Thereby obtaining the zenith angle theta and the azimuth angle phi of the track; dividing angle intervals on zenith angle theta and azimuth angle phi and counting the number of cosmic ray muon in each angle interval to obtain a flux-direction distribution histogram parameter muoncounts (theta, phi) of the cosmic ray muon in each position under the coal mine goaf model to be detected and a flux-direction distribution histogram parameter muonopenky (theta, phi) of the cosmic ray muon of the open sky above the coal mine goaf model to be detected; thereby obtaining the residual flux ratio of the cosmic ray muon Experiment (θ,φ)=muoncounts(θ,φ)/muonopensky(θ,φ),ratio Experiment (θ, φ) can characterize the opacity parameter of the muon detector in various directions across the coal mine goaf model;
step 2, based on the ratio observed in step 1 Experiment Reconstructing a goaf model, arranging the positions of muon detectors in a step 1, designating the accuracy dimensions grid_x, grid_y and grid_z of the reconstructed goaf model according to the size of an observation area, dividing the observation area into three-dimensional grids which are totally Nx multiplied by Ny multiplied by Nz grid solving areas to form a coal mine goaf three-dimensional grid lattice; then, randomly setting the grid of the three-dimensional grid lattice of the goaf of the coal mine, which is filled by the coal rock mass, as 1, and not being filledThe grid of the goaf is set to 0, and the initial model form of the coal mine goaf is randomly generated;
step 3, model population gene coding and goaf model population t generation: the goaf model constructed according to the step 2 has Nx Ny x Nz grid density parameters which are respectively replaced by 1 and 0, and the gene length is Nx Ny x Nz; generating a population t with N individuals by randomly changing the arrangement sequence of 0 and 1, wherein the population is the goaf initial morphology model population;
step 4, creating Nx, ny and Nz geometric models with the size of grid_x, grid_y, grid_z in the Geant4 simulation program, performing density assignment according to the population genes in the step 3 to form a goaf physical model, and setting a muon detector at the same observation position as the step 1; using CRY cosmic ray generation sub-model cosmic ray source to calculate muoncounts (theta, phi) received by muon detectors arranged at observation positions under goaf model to be measured, and then calculating ratio corresponding to each observation position Model (θ,φ);
Step 5: population adaptation value calculation: repeating the step 4 until the calculation of N model goaf individuals in the population is completed, and obtaining the histogram parameter ratio of N groups of model residual flux distribution Model (θ, φ); the adaptation value alpha is defined as the histogram parameter ratio of the model residual flux distribution Model Histogram parameter ratio of (θ, φ) to experimental residual flux distribution Experiment Variance of (θ, φ):
Figure 880387DEST_PATH_IMAGE001
wherein i and j are the number of angle interval divisions of zenith angle theta and azimuth angle phi respectively;
step 6, excellent population individual hybridization variation: the adaptation values of N individuals in the step 5 are sorted from small to large, 2N/3 individuals are removed, the density chromosomes of the rest N/3 individuals are subjected to hybridization variation to generate new 2N/3 individuals and the original winning N/3 individuals to form a new generation population t+1 of N individuals;
step 7, repeating the steps 4-6, outputting the density genes of N/3 partial population individuals when the adaptation values of N/3 individuals in the population are smaller than the designed relative error threshold, and decoding the density genes as a coal mine goaf three-dimensional density structure model meeting the conditions;
and 8, the N/3 three-dimensional density structure models obtained in the step 7 are not unique, so that the most reasonable coal mine goaf three-dimensional density structure model and corresponding uncertainty are obtained by continuously adopting a bootstrap resampling method.
The invention also comprises the following technical characteristics:
specifically, in the step 1, the coal mine goaf model is a natural coal mine goaf or a manually placed square block with different densities.
Specifically, the muon detector is a plastic scintillator array detector, a gas detector or a nuclear latex detector.
Specifically, in the step 1, the zenith angle θ and the azimuth angle Φ:
Figure 138193DEST_PATH_IMAGE002
Figure 235548DEST_PATH_IMAGE003
the method comprises the steps of carrying out a first treatment on the surface of the Wherein x=x in -x out ;y=y in -y out ;z=z in -z out
Specifically, in the step 2, the typical density of the grid corresponding to the coal rock mass 1 is 1.4g/cm 3 A typical density of the grid of 0 corresponding to the goaf is 0.5g/cm 3
Specifically, in the step 8, the bootstrap resampling method includes:
the first step is that K resampling is carried out on N/3 population data to form a plurality of bootstrap data sets, each resampling is random sampling, and the number of samples is the same as the number of the population data set elements;
the second step is that on the basis of an inversion frame, K bootstrap data sets are utilized to carry out K times of Geant4 modeling calculation to obtain an inversion result combination consisting of K density models, and the initial model, relevant inversion parameters and iteration times required by each calculation are the same;
and thirdly, carrying out statistical analysis and evaluation on the calculation result combination, indicating the uncertainty of the model according to the difference condition among the bootstrap calculation models, and finally giving out the most reasonable coal mine goaf three-dimensional density structure model and the corresponding uncertainty.
Compared with the prior art, the invention has the following technical effects:
according to the invention, a model most conforming to experimental observation data is searched by utilizing a principle of sensitivity of cosmic ray muon to opacity (density length) of a substance and utilizing a Meng Daka-roar simulation mode through a genetic algorithm; the method overcomes the limitation of the traditional technology in complex goaf detection, and can give out a three-dimensional density structural model of the goaf of the coal mine and corresponding uncertainty.
Drawings
FIG. 1 is a schematic model of an artificially constructed simulated goaf;
FIG. 2 is a side view of a coal mine goaf model;
FIG. 3 is a randomly generated goaf model;
FIG. 4 is a density inversion model obtained by decoding calculation.
Detailed Description
The invention provides a goaf scale measurement calculation method based on muir remote sensing data, which is a reconstruction method of the shape and position of a coal mine goaf, adopts a genetic algorithm solving mode, divides a model of the initial shape of the coal mine goaf into grids, carries out coding according to a density normal (1) and a density loss (0) to construct a model population, respectively calculates the adaptation value of an individual through a Monte Carlo program Geant4, generates a new population through multiple hybridization variation of the individual, obtains a model solution set meeting the condition through the genetic algorithm, gives out a final model solution and model uncertainty on the basis of the model solution set through the bootstrap algorithm, and finally carries out inverse coding according to the density normal (1) and the density loss (0) to give out the most reasonable three-dimensional density structure model and uncertainty of the coal mine goaf. The method overcomes the limitation of the traditional technology in complex goaf detection, and can give out a three-dimensional density structural model of the goaf of the coal mine and corresponding uncertainty.
The following specific embodiments of the present invention are provided, and it should be noted that the present invention is not limited to the following specific embodiments, and all equivalent changes made on the basis of the technical solutions of the present application fall within the protection scope of the present invention.
Example 1:
the embodiment provides a goaf scale measurement calculation method based on muon remote sensing data, which comprises the following steps:
step 1, laying a muon detector below a coal mine goaf model, and recording the incident position coordinates (x in ,y in ,z in ) And the exit position coordinates (x out ,y out ,z out ) The zenith angle θ and azimuth angle Φ of the track are calculated:
Figure 273912DEST_PATH_IMAGE004
Figure 283456DEST_PATH_IMAGE005
the method comprises the steps of carrying out a first treatment on the surface of the Wherein x=x in -x out ;y=y in -y out ;z=z in -z out The method comprises the steps of carrying out a first treatment on the surface of the Dividing angle intervals on zenith angle theta and azimuth angle phi and counting the number of cosmic ray muon in each angle interval to obtain a flux-direction distribution histogram parameter muoncounts (theta, phi) of the cosmic ray muon in each position under the coal mine goaf model to be detected and a flux-direction distribution histogram parameter muonopenky (theta, phi) of the cosmic ray muon of the open sky above the coal mine goaf model to be detected; thereby obtaining the residual flux ratio of the cosmic ray muon: ratio of Experiment (θ, φ) =muoncounts (θ, φ)/muonopensky (θ, φ), a histogram parameter ratio of the residual flux distribution of cosmic ray muon Experiment (θ, φ) can characterize the opacity parameter of the muon detector in various directions across the coal mine goaf model;
in the step 1, a coal mine goaf model is a natural coal mine goaf or a manually placed square block with different densities; the muon detector is a plastic scintillator array detector, a gas detector or a nuclear latex detector; the zenith angle theta ranges from 0 to 90 degrees, and the angle interval can be divided into 18 angle intervals according to 5-degree division; the azimuth angle phi ranges from 0 to 360 degrees, and the angle interval can be divided into 10 angle intervals according to 36-degree division, so that the whole detection space is divided into 18 x 10 = 180 angle intervals;
specifically, in this embodiment, observation data of the muon is first obtained, and the data may be obtained according to actual observation or may be obtained through analog calculation; the embodiment describes a data acquisition mode in the actual observation process, wherein the used equipment is a muon detector with position resolution, and the mode of acquiring data by muon equipment data is to record the concentrated incident position coordinates (x in ,y in ,z in ) Emission position coordinates (x out ,y out ,z out ) The determination of the position coordinates of the muon can be obtained through detection by various detectors (such as a plastic scintillator array detector, a gas detector, a nuclear latex detector and the like), and although the data acquisition mode is not unique, the finally-given data mainly comprises the position information of the muon; the method protects the data processing process after the position of the cosmic ray muon is acquired;
the experimental data can be obtained by carrying out experiments in natural coal mine goaf, or by manually setting squares with different densities to construct a coal mine goaf model for experiments; when the test starts, a cosmic ray position sensitive monitoring sensor is respectively arranged below the natural goaf model or the artificial goaf model; record the incidence of the received cosmic ray (x in ,y in ,z in ) And the exit position (x out ,y out ,z out ) The method comprises the steps of carrying out a first treatment on the surface of the Fig. 1 shows a schematic model of a simulated goaf constructed manually, fig. 2 shows a side view of a goaf model of a coal mine, wherein a gray area is a bedrock without abnormality, a white area is a goaf, and p1-p5 are 5 observation positions respectively.
The muon monitoring sensors at 5 observation positions record the incidence of the received cosmic rays respectivelyx in ,y in ,z in ) And the exit position (x out ,y out ,z out ) And respectively calculating the zenith angle theta and the azimuth angle phi of the track:
Figure 599031DEST_PATH_IMAGE006
Figure 668618DEST_PATH_IMAGE007
wherein x=x in -x out ;y=y in -y out ;z=z in -z out
Then dividing different intervals according to the angle range, respectively counting the number of muon in each angle interval (for example, the theta range is 0-90 degrees, the angle interval range can be divided into 18 intervals according to the 5-degree division, the phi range is 0-360 degrees, the angle interval range can be divided into 10 intervals according to the 36-degree division, then the whole detection space is divided into 18 x 10 = 180 angle intervals), and obtaining a distribution histogram parameter muoncounts (theta, phi) of the flux of the cosmic ray at each position under the model to be detected along with the direction. The histogram parameter muonopsesky (theta, phi) of the distribution of the cosmic ray flux of the open sky above the model to be measured along with the direction can be measured in the same way, and the residual flux ratio of the cosmic ray is further calculated: ratio Experiment (θ, φ) =muoncounts (θ, φ)/muonopsensky (θ, φ), a histogram parameter ratio of the cosmic ray residual flux distribution Experiment (θ, φ) is used to characterize the opacity parameter across the physical similarity model in various directions at the detector location.
Step 2, reconstructing a goaf model based on histogram parameters (theta, phi) of the residual flux distribution of the cosmic ray observed in the step 1, arranging the positions of muon detectors in a mode of the step 1, designating grid precision dimensions grid_x, grid_y and grid_z of the reconstructed goaf model according to the size of an observation area, dividing the observation area into three-dimensional grids, and forming a coal mine goaf three-dimensional grid lattice by using Nx Ny Nz grid solving areas; then, randomly lattice the three-dimensional grid of the goaf of the coal mineThe grid filled by the coal rock mass is set as 1, the grid of the goaf not filled is set as 0, the initial model form of the goaf of the coal mine is randomly generated, and the typical density of the grid of the 1 corresponding to the coal rock mass is 1.4g/cm 3 0 corresponds to an unfilled grid typical density of 0.5g/cm 3
As shown in fig. 3, the positions of the muon detectors P are arranged in the same manner as the test, P1-P5 are respectively 5 observation positions, and according to the reconstructed grid precision dimensions grid_x, grid_y and grid_z specified by the sizes of the observation regions, the three-dimensional grid is divided into a three-dimensional grid solving region (taking the initial form physical model grid of the coal mine goaf of fig. 3 as nx=27, ny=1 and nz=9 as an example) to form a three-dimensional grid lattice of the coal mine goaf. Then, randomly setting the grid of the three-dimensional grid lattice of the coal mine goaf filled by the coal rock mass as 1, and setting the unfilled grid of the unfilled goaf as 0; model initial morphology of randomly generated coal mine goaf (a model randomly generated is shown in FIG. 3), typical density of 1 corresponds to 1.4g/cm 3 0 corresponds to an unfilled grid typical density of 0.5g/cm 3
Step 3, model population gene coding and goaf model population t generation: the goaf model constructed according to the step 2 has Nx Ny x Nz grid density parameters which are respectively replaced by 1 and 0, and the gene length is Nx Ny x Nz; generating a population t with N individuals by randomly changing the arrangement sequence of 0 and 1, wherein the population is the goaf initial morphology model population; for example, fig. 2 shows an individual in a randomly generated population t.
Population gene coding: the biological DNA gene carries genetic information of organisms, but the coding units of the biological DNA gene are only 4 (adenine deoxyribonucleotide, guanine deoxyribonucleotide, cytosine deoxyribonucleotide and thymine deoxyribonucleotide), and various genetic information can be formed through different coded permutation and combination, and finally the biological DNA gene is expressed as diversity of the organisms. The code of the invention is similar to the above, and is a mode of corresponding density value by digital code, the typical density of the code 1 corresponding to the coal rock mass in the scheme is 1.4g ∈ -cm 3 The typical density of the grid filled with the goaf corresponding to the code 0 is 0.5g/cm 3
Goaf model population: after the coding mode is determined, different sequences with the length of Nx Ny Nz are generated through different sequences of 0 and 1 in a mode of step 2, coding information carried by one sequence represents a goaf density structure model, and the goaf density structure model is equivalent to an individual in the set of goaf density structures; the collection of these goaf model individuals is referred to as a goaf model population.
Step 4, creating Nx, ny and Nz geometric models with the size of grid_x, grid_y, grid_z in the Geant4 simulation program, performing density assignment according to the population genes in the step 3 to form a goaf physical model, and setting a muon detector at the same observation position as the step 1; using CRY cosmic ray generation sub-model cosmic ray sources to calculate muoncounts (theta, phi) received by muon detectors arranged at observation positions under goaf models to be detected, and then calculating ratio models (theta, phi) corresponding to each observation position;
specifically, in this embodiment, nx×ny×nz geometric models with dimensions of grid_x×grid_y×grid_z are created in the Geant4 simulation program, and density assignment is performed to form a goaf physical model according to "density genes" in step 3 (genes in biology are nucleotide sequences carrying genetic information; the expression mode of biology is referred to in the invention by reference to biology, to describe a length nx×ny×nz coding sequence generated by goaf model individuals according to the coding rules, because the coding values correspond to density information, they are called density genes), and a cosmic line detector is set at the same p1-p5 position as in step 1. Using CRY cosmic ray generation sub-model cosmic ray source to calculate the muoncounts (theta, phi) received by the cosmic ray detector set at the p1-p5 position under the goaf model to be measured, and then respectively calculating the corresponding ratio Model (θ,φ)。
Step 5, population adaptation value calculation: repeating the step 4 until the calculation of N model goaf individuals in the population is completed, and obtaining the orthometric of N groups of model residual flux distributionA graph parameter ratio model (θ, φ); the adaptation value alpha is defined as the histogram parameter ratio of the model residual flux distribution Model Histogram parameter ratio of (θ, φ) to experimental residual flux distribution Experiment Variance of (θ, φ):
Figure 194277DEST_PATH_IMAGE008
in the above expression, i denotes the number of bins of zenith angle θ, and j denotes the number of bins of azimuth angle Φ.
Step 6, excellent population individual hybridization variation: the adaptation values of N individuals in the step 5 are sorted from small to large, 2N/3 individuals are removed, the density chromosomes of the rest N/3 individuals are subjected to hybridization variation to generate new 2N/3 individuals and the original winning N/3 individuals to form a new generation population t+1 of N individuals;
step 7, repeating the steps 4-6, outputting the density genes of N/3 partial population individuals when the adaptation values of N/3 individuals in the population are smaller than the designed relative error threshold, and decoding the density genes as a coal mine goaf three-dimensional density structure model meeting the conditions;
and 8, the N/3 three-dimensional density structure models obtained in the step 7 are not unique, so that a bootstrap resampling method is continuously adopted to obtain the most reasonable coal mine goaf three-dimensional density structure model and corresponding uncertainty, and a density inversion model obtained through decoding calculation is shown in fig. 4. Specifically, the bootstrap resampling method comprises the following steps: the first step is that K resampling is carried out on N/3 population data to form a plurality of bootstrap data sets, each resampling is random sampling, and the number of samples is the same as the number of the population data set elements; the second step is that on the basis of an inversion frame, K bootstrap data sets are utilized to carry out K times of Geant4 modeling calculation to obtain an inversion result combination consisting of K density models, and the initial model, relevant inversion parameters and iteration times required by each calculation are the same; and thirdly, carrying out statistical analysis and evaluation on the calculation result combination, indicating the uncertainty of the model according to the difference condition among the bootstrap calculation models, and finally giving out the most reasonable coal mine goaf three-dimensional density structure model and the corresponding uncertainty.

Claims (6)

1.一种基于缪子遥感数据的采空区规模测量计算方法,其特征在于,包括以下步骤:1. A goaf scale measurement calculation method based on Muzi remote sensing data, is characterized in that, comprises the following steps: 步骤1,在煤矿采空区模型下方布设缪子探测器,记录宇宙线缪子在缪子探测器的入射位置坐标(xin,yin,zin)和出射位置坐标(xout,yout,zout),从而得到径迹的天顶角θ和方位角φ;对天顶角θ和方位角φ划分角度区间并统计每个角度区间的宇宙线缪子数量,得到待测煤矿采空区模型下每个位置处的宇宙线缪子的通量随方向分布直方图参数muoncounts(θ,φ),以及待测煤矿采空区模型上方的开阔天空的宇宙线缪子的通量随方向分布直方图参数muonopensky(θ,φ);从而得到宇宙线缪子剩余通量比ratio实验(θ,φ)=muoncounts(θ,φ)/muonopensky(θ,φ),ratio实验(θ,φ)能表征缪子探测器位置处各个方向上穿过煤矿采空区模型的不透明度参数;Step 1: Arrange muon detectors under the coal mine goaf model, and record the incident position coordinates (x in , y in , z in ) and exit position coordinates (x out , y out ) of cosmic ray muons in the muon detector , z out ), so as to obtain the zenith angle θ and azimuth φ of the track; divide the zenith angle θ and azimuth φ into angle intervals and count the number of cosmic ray muons in each angle interval, and obtain the goaf of the coal mine to be tested The flux of cosmic ray muons at each location under the area model varies with the direction distribution histogram parameter muoncounts(θ, φ), and the flux of cosmic ray muons in the open sky above the goaf model of the coal mine to be tested varies with the direction Distribution histogram parameter muonopensky(θ, φ); thus obtain cosmic ray muon residual flux ratio ratio experiment (θ, φ)=muoncounts(θ, φ)/muonopensky(θ, φ), ratio experiment (θ, φ) It can characterize the opacity parameters passing through the coal mine goaf model in all directions at the position of the muon detector; 步骤2,基于步骤1观测到的ratio实验(θ,φ)重构采空区模型,按照步骤1的方式布置缪子探测器的位置,并按照观测区域的大小指定重构的采空区模型网格精度尺寸grid_x,grid_y,grid_z,将观测区域划分三维网格共Nx×Ny×Nz个网格求解区域形成煤矿采空区三维网格点阵;然后,随机将煤矿采空区三维网格点阵被煤岩体填充的网格设置为1,未被填充的采空区的网格设置为0,随机生成煤矿采空区的模型初始形态;Step 2: Reconstruct the goaf model based on the ratio experiment (θ, φ) observed in step 1, arrange the position of the muon detector according to the method in step 1, and specify the reconstructed gob model according to the size of the observation area The grid precision size grid_x, grid_y, grid_z, the observation area is divided into three-dimensional grids, a total of Nx×Ny×Nz grid solution area to form a three-dimensional grid lattice of the coal mine goaf; then, the three-dimensional grid of the coal mine goaf is randomly divided into The grid of the lattice filled with coal and rock mass is set to 1, the grid of the unfilled goaf is set to 0, and the initial shape of the model of the coal mine goaf is randomly generated; 步骤3,模型种群基因编码及生成采空区模型种群t:按照步骤2构建的采空区模型具有Nx×Ny×Nz个网格密度参数,分别由1、0代替,基因长度为Nx×Ny×Nz;通过随机改变0、1的排列顺序生成具有N个个体的种群t,此种群即为采空区初始形态模型种群;Step 3, model population gene encoding and generation of goaf model population t: the goaf model constructed according to step 2 has Nx×Ny×Nz grid density parameters, which are replaced by 1 and 0 respectively, and the gene length is Nx×Ny ×Nz; by randomly changing the arrangement order of 0 and 1, a population t with N individuals is generated, which is the initial shape model population of the goaf; 步骤4,在Geant4模拟程序中创建Nx×Ny×Nz个尺寸为grid_x×grid_y×grid_z的几何模型,并按照步骤3中的种群基因进行密度赋值形成采空区物理模型,并在与步骤1相同的观测位置处设置缪子探测器;使用CRY宇宙线产生子模拟宇宙线源计算出在待测采空区模型下的观测位置处设置缪子探测器接收到的muoncounts(θ,φ),然后计算每个观测位置对应的ratio模型(θ,φ);Step 4, create Nx×Ny×Nz geometric models with the size of grid_x×grid_y×grid_z in the Geant4 simulation program, and perform density assignment according to the population gene in step 3 to form a physical model of the goaf, and perform the same as step 1 Set the muon detector at the observation position; use the CRY cosmic ray generator to simulate the cosmic ray source to calculate the muoncounts (θ, φ) received by the muon detector at the observation position under the goaf model to be tested, and then Calculate the ratio model (θ, φ) corresponding to each observation position; 步骤5:种群适应值计算:重复步骤4,直至完成种群中N个模型采空区个体的计算,得到N组模型剩余通量分布的直方图参数ratio模型(θ,φ);适应值α定义为模型剩余通量分布的直方图参数ratio模型(θ,φ)与实验剩余通量分布的直方图参数ratio实验(θ,φ)的方差:Step 5: Calculation of fitness value of the population: Repeat step 4 until the calculation of goaf individuals of N models in the population is completed, and the histogram parameter ratio model (θ, φ) of the residual flux distribution of N groups of models is obtained; the definition of fitness value α is the variance of the histogram parameter ratio model (θ, φ) of the model residual flux distribution and the histogram parameter ratio experiment (θ, φ) of the experimental residual flux distribution:
Figure 28434DEST_PATH_IMAGE001
Figure 28434DEST_PATH_IMAGE001
其中,i、j分别为天顶角θ和方位角φ的角度区间划分数量;Among them, i and j are the number of angular interval divisions of zenith angle θ and azimuth φ respectively; 步骤6,优秀种群个体杂交变异:步骤5中的N个个体的适应值按照从小到大排序,剔除掉后2N/3个体,将剩下的N/3个体的密度染色体进行杂交变异生成新的2N/3个体与原来优胜N/3个体组成新一代N个个体的种群t+1;Step 6, hybridization and variation of excellent population individuals: in step 5, the fitness values of the N individuals are sorted from small to large, and 2N/3 individuals are eliminated, and the density chromosomes of the remaining N/3 individuals are hybridized and mutated to generate new 2N/3 individuals and the original winning N/3 individuals form a population t+1 of N individuals in a new generation; 步骤7,重复步骤4~步骤6,当种群中N/3个个体的适应值小于所设计的相对误差阈值时,则输出这N/3部分种群个体的密度基因,解码作为符合条件的煤矿采空区三维密度结构模型;Step 7, repeat steps 4 to 6, when the fitness value of N/3 individuals in the population is less than the designed relative error threshold, then output the density genes of the N/3 population individuals and decode them as eligible coal mining Three-dimensional density structure model of void area; 步骤8,步骤7中得到的N/3个三维密度结构模型不具有唯一性,因此继续采用bootstrap重采样方法,得到最合理的煤矿采空区三维密度结构模型及相应的不确定度。Step 8, the N/3 3D density structure models obtained in step 7 are not unique, so continue to use the bootstrap resampling method to obtain the most reasonable 3D density structure model and corresponding uncertainty of the coal mine goaf.
2.如权利要求1所述的基于缪子遥感数据的采空区规模测量计算方法,其特征在于,所述步骤1中,煤矿采空区模型为天然的煤矿采空区或人工摆放的不同密度方块。2. the goaf scale measurement calculation method based on the remote sensing data of Muzi as claimed in claim 1, is characterized in that, in described step 1, the coal mine goaf model is a natural coal mine goaf or artificially placed Blocks of different densities. 3.如权利要求1所述的基于缪子遥感数据的采空区规模测量计算方法,其特征在于,所述缪子探测器为塑料闪烁体阵列探测器、气体探测器或核乳胶探测器。3. The goaf scale measurement and calculation method based on muon remote sensing data as claimed in claim 1, wherein the muon detector is a plastic scintillator array detector, a gas detector or a nuclear latex detector. 4.如权利要求1所述的基于缪子遥感数据的采空区规模测量计算方法,其特征在于,所 述步骤1中,天顶角θ和方位角φ:
Figure 428322DEST_PATH_IMAGE002
Figure 917072DEST_PATH_IMAGE003
;其中,x=xin-xout;y= yin-yout;z=zin-zout
4. the goaf scale measurement and calculation method based on Muzi remote sensing data as claimed in claim 1, is characterized in that, in described step 1, zenith angle θ and azimuth angle φ:
Figure 428322DEST_PATH_IMAGE002
;
Figure 917072DEST_PATH_IMAGE003
; among them, x=x in -x out ;y= y in -y out ;z=z in -z out .
5.如权利要求1所述的基于缪子遥感数据的采空区规模测量计算方法,其特征在于,所述步骤2中,1对应煤岩体的网格典型密度为1.4g/cm3,0对应采空区的网格典型密度为0.5g/cm35. The method for measuring and calculating the scale of gobs based on remote sensing data according to claim 1, characterized in that, in said step 2, the grid typical density corresponding to 1 is 1.4g/cm 3 , 0 corresponds to a typical grid density of 0.5g/cm 3 in the goaf. 6.如权利要求1所述的基于缪子遥感数据的采空区规模测量计算方法,其特征在于,所述步骤8中,采用bootstrap重采样方法包括:6. the goaf scale measurement and calculation method based on Muzi remote sensing data as claimed in claim 1, is characterized in that, in described step 8, adopts bootstrap resampling method to comprise: 第一步是对N/3种群数据进行K次重采样构成多个bootstrap数据集,每次重采样都是随机采样且采样个数与种群数据集元素个数相同;The first step is to resample the N/3 population data K times to form multiple bootstrap data sets. Each resampling is random sampling and the number of samples is the same as the number of elements in the population data set; 第二步是在反演框架的基础上利用K个bootstrap数据集进行K次Geant4建模计算得到由K个密度模型组成的反演结果组合,每次计算所需的初始模型、相关反演参数以及迭代次数都应相同;The second step is to use K bootstrap data sets to perform K Geant4 modeling calculations on the basis of the inversion framework to obtain a combination of inversion results consisting of K density models. The initial model and related inversion parameters required for each calculation and the number of iterations should be the same; 第三步则是对该计算结果组合进行统计分析和评价,根据bootstrap计算模型之间的差异情况指示模型不确定度,最终给出最合理的煤矿采空区三维密度结构模型及相应的不确定度。The third step is to conduct statistical analysis and evaluation on the combination of calculation results, and indicate the model uncertainty according to the difference between the bootstrap calculation models, and finally give the most reasonable three-dimensional density structure model of the coal mine goaf and the corresponding uncertainty Spend.
CN202310022850.0A 2023-01-08 2023-01-08 Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data Active CN115935120B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310022850.0A CN115935120B (en) 2023-01-08 2023-01-08 Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310022850.0A CN115935120B (en) 2023-01-08 2023-01-08 Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data

Publications (2)

Publication Number Publication Date
CN115935120A CN115935120A (en) 2023-04-07
CN115935120B true CN115935120B (en) 2023-06-13

Family

ID=86552561

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310022850.0A Active CN115935120B (en) 2023-01-08 2023-01-08 Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data

Country Status (1)

Country Link
CN (1) CN115935120B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108986208A (en) * 2018-07-11 2018-12-11 辽宁工程技术大学 A kind of reconstructing method of coal mine gob inbreak form
CN111789614A (en) * 2020-08-10 2020-10-20 上海联影医疗科技有限公司 An imaging system and method
CN115542410A (en) * 2022-09-01 2022-12-30 兰州大学 Target scanning and three-dimensional forward and backward modeling method based on cosmic ray muon

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI129611B (en) * 2019-08-23 2022-05-31 Muon Solutions Oy System and method for material density distribution survey based on cosmic muon detection

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108986208A (en) * 2018-07-11 2018-12-11 辽宁工程技术大学 A kind of reconstructing method of coal mine gob inbreak form
CN111789614A (en) * 2020-08-10 2020-10-20 上海联影医疗科技有限公司 An imaging system and method
CN115542410A (en) * 2022-09-01 2022-12-30 兰州大学 Target scanning and three-dimensional forward and backward modeling method based on cosmic ray muon

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
智宇;周静;陈雷;李沛玉;赵明锐;刘雯迪;贾世海;张昀昱;胡守扬;于伟翔;李笑梅.宇宙线缪子散射成像模拟与算法研究.原子能科学技术.2020,1-8. *

Also Published As

Publication number Publication date
CN115935120A (en) 2023-04-07

Similar Documents

Publication Publication Date Title
Yin et al. Simulation-based investigation on the accuracy of discrete fracture network (DFN) representation
CN108897066A (en) Carbonate rock crack density quantitative prediction method and device
CN113221228B (en) Hydropower station underground cave group surrounding rock mechanical parameter inversion method
CN113808265B (en) Construction method of modeling section of three-dimensional geological model of large-scale area
CN117665932A (en) Wave speed inversion method based on potential spatial features of three-dimensional seismic wave data of tunnel
CN118332902A (en) Three-dimensional numerical simulation method for seepage field evolution in fractured rock tunnels
CN115935120B (en) Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data
CN114254566A (en) A Neural Network Geological Exploration Inversion Method Based on 2.5-Dimensional Hybrid Spectral Element Method
CN118116506B (en) Underground lithology recognition method and system based on machine learning
Martinius et al. Uncertainty analysis of fluvial outcrop data for stochastic reservoir modelling
Lachérade et al. A sectorisation-based method for geostatistical modeling of pressuremeter test data: Application to the Grand Paris Express project (France)
CN118151236A (en) Seismic velocity deep learning inversion method and system based on data-space domain fusion
M. Tóth et al. Relationship between the geometric parameters of rock fractures, the size of percolation clusters and REV
CN115618697B (en) CAE simulation inversion method and server for multiple stratum parameters of shield tunnel
Libicki et al. Stochastic branching models of fault surfaces and estimated fractal dimensions
CN114609691B (en) Method, device, computer equipment and storage medium for determining geomorphic causal model
Zhao et al. Automatic identification of concealed dangerous rock blocks on high-steep slopes considering finite-sized discontinuity intersections
Penna et al. Seismic quantitative interpretation of diagenetic processes: A study in Brazilian pre-salt carbonate reservoirs
Parasyris et al. Near surface full waveform inversion via deep learning for subsurface imaging
US11422276B2 (en) Method and apparatus for estimating S-wave velocities by learning well logs
CN119004931B (en) Discrete element-based natural gas hydrate-containing sediment simulation method and system
Xu Identifying rock blocks based on hierarchical rock-mass structure model
CN117634235B (en) A hybrid Monte Carlo deep penetration calculation method for through-casing density logging
Abatchev Development of Specialized Nonlinear Inversion Algorithms, Basis Functions, Eikonal Solvers, and Their Integration for Use in Joint Seismic and Gravitational Tomographic Inversion
CN120233453A (en) Manganese ore deposit investigation effect evaluation method and system based on spread spectrum excitation method

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