CN102508310B - Detection method for porosity distribution of upper formation of fire district of coal field - Google Patents
Detection method for porosity distribution of upper formation of fire district of coal field Download PDFInfo
- Publication number
- CN102508310B CN102508310B CN 201110335168 CN201110335168A CN102508310B CN 102508310 B CN102508310 B CN 102508310B CN 201110335168 CN201110335168 CN 201110335168 CN 201110335168 A CN201110335168 A CN 201110335168A CN 102508310 B CN102508310 B CN 102508310B
- Authority
- CN
- China
- Prior art keywords
- porosity
- factor
- survey line
- fire district
- borings
- 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.)
- Expired - Fee Related
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 14
- 239000003245 coal Substances 0.000 title abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims description 8
- 230000001939 inductive effect Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 239000004575 stone Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000011065 in-situ storage Methods 0.000 abstract 1
- 230000001902 propagating effect Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 23
- 238000002485 combustion reaction Methods 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Disclosed is a detection method for porosity distribution of an upper formation of a fire district of a coal field. A plurality of measuring lines are arranged on the ground surface of the upper portion of the fire district of the coal field, a plurality of measuring points are disposed on each measuring line at intervals, holes are constructed and drilled on the measuring points on each measuring line one by one, then a transmitting antenna and a receiving antenna are respectively disposed in each two adjacent drilled holes on the same measuring line one by one, drilled hole radar multi-offset-distance cross-hole transmission detection is sequentially carried out for all the drilled holes on each measuring line, and the time span of electromagnetic waves propagating from the transmitting antennas to the receiving antennas under conditions that offset distances are different are acquired; and data of the acquired propagation time span of the electromagnetic waves are processed by the aid of a computer, mosaic and combination of porosity data on a two-dimensional section between each two adjacent drilled holes on each measuring line is realized, and a contour map of porosity distribution of the formation on the section of each measuring line is drawn. Accordingly, porosity distribution of the upper formation of the fire district of the coal field is obtained, and the condition of the porosity distribution of the loose and broken formation in an upper space of the fire district of the coal field is quickly and effectively detected in an in-situ state.
Description
Technical field
The present invention relates to a kind of detection method, especially a kind of detection method that is applicable to that fire district, coalfield upper formation factor of porosity distributes.
Background technology
It is the air feed of fire district, research coalfield combustion system and the basis of air-out passage and the migration rule of firing-fighting medium in the scall of Huo Qu top that the factor of porosity of fire district, coalfield upper formation distributes.But also there is not to survey under a kind of state in position the method that the scall factor of porosity distributes at present.The measurement of formation porosity is used more in oil-gas geology, mainly contains well logging method (as acoustic logging, density logging and nuclear magnetic resonance log), seismic inversion method.But the scope that acoustic logging, density logging and nuclear magnetic resonance log method are surveyed is very limited, can only on vertically, estimate stratum pores'growth situation in the position that drilling well is arranged, and can not show factor of porosity situation of change in the horizontal, therefore be not suitable for the detection of fire district, the coalfield upper formation factor of porosity that alters a great deal on the space.The seismic inversion method has obtained more application in the petroleum and gas geology and exploration at home and abroad in recent years, and this method is set up sound impedance or speed and reservoir porosity relation, thereby asked for factor of porosity by the rock physics analysis.The seismic inversion method is primarily aimed at the darker stratum of buried depth, but the resolution for the more shallow and looser ravel area stratum, Huo Qu top of buried depth is relatively poor, and this method need be carried out shot firing operation in the survey district, can cause further subsiding of Huo Qu subsidence area, even cause casualties.The seismic inversion method also has application in conjunction with the method for well-log information joint inversion formation porosity, this method is by extracting seismic properties, set up the relation of seismic properties and each well logging factor of porosity, utilize nonlinear neural network and multiple linear ASSOCIATE STATISTICS etc., the factor of porosity variation that obtains the whole district thereby the factor of porosity that quantitatively calculates no wellblock distributes.But because the seismic properties of extracting and factor of porosity are difficult to directly establish the relational expression of definite meaning, so bigger to the error of calculation of formation porosity.As seen through the above analysis, well logging method and seismic inversion method all can not be used for changing on the space exploration factor of porosity distribution of big and loose and broken formation.
Summary of the invention
The purpose of this invention is to provide a kind of method that fire district, coalfield upper formation factor of porosity distributes of surveying, to be implemented in the distribution of surveying fire district, coalfield upper formation factor of porosity under the original position state.
The detection method that fire district, coalfield of the present invention upper formation factor of porosity distributes comprises the steps:
A, arrange many surveys line at the ground surface on Huo Qu top, coalfield, on between each survey line every a plurality of measuring points are set, one by one to the measuring point construction drill on the every survey line, on same survey line, in adjacent two borings emitting antenna and receiving antenna are set respectively one by one then, successively the radar multiple offset distance of holing of all borings on the every survey line is striden the hole transmission and surveyed by emitting antenna and receiving antenna, be captured in the travel-time of electromagnetic wave from the emitting antenna to the receiving antenna under the different offset distance situations;
B, with electromagnetic wave travel-time of collecting by Computer Processing, obtain on the every survey line slowness of each pixel on the two dimensional cross-section between adjacent two borings
, according to formula:
Calculate on the every survey line bulk dielectric constant of each pixel on the two dimensional cross-section between adjacent two borings
In the formula:
Be electromagnetic wave velocity of propagation in a vacuum;
The not burnt rock stone sample of c, fire district, collection coalfield upper formation is with the specific inductive capacity of decameter measurement rock sample relative vacuum
D, utilize formula:
Calculate the factor of porosity of each pixel on the two dimensional cross-section between two whenever adjacent borings
By computing machine with on each bar survey line between adjacent two borings the factor of porosity data on the two dimensional cross-section splice merging respectively, and make the isogram that every formation porosity on the survey line section distributes, thereby draw the distribution of fire district, coalfield upper formation factor of porosity, the more big then stratum of factor of porosity is more loose, and the more little then stratum of factor of porosity is more closely knit; In the formula:
Be the specific inductive capacity of air,
It is a factor relevant with direction of an electric field.
Described ground surface on Huo Qu top, coalfield arranges that the spacing of many surveys line is 50m; The described spacing that a plurality of measuring points is set at interval at each survey line is 40 ~ 50m.
Beneficial effect: the present invention adopts under the state boring radar that fire district upper formation is carried out multiple offset distance in position and strides the hole transmission measurement, is captured in the travel-time of electromagnetic wave from the emitting antenna to the receiving antenna under the different offset distance situations; By computing machine the gained data are handled, obtained the slowness of each pixel on the two dimensional cross-section between adjacent two borings
And bulk dielectric constant
Gather Huo Qu top and do not burn the rock sample that becomes the stratum, and measure the specific inductive capacity of its relative vacuum
Utilize formula:
Calculate the factor of porosity of each pixel on the two dimensional cross-section between the two adjacent borings
By computing machine with on each bar survey line between the adjacent boring factor of porosity data on the two dimensional cross-section splice merging respectively, and make the isogram that every formation porosity on the survey line section distributes, thereby detect the distribution of Huo Qu top, coalfield scall factor of porosity, the more big representative of factor of porosity stratum is more loose, and the more little stratum of then representing of factor of porosity is more closely knit.The inventive method is easy, data are accurate, compare with logging methods such as acoustic logging, density logging and nuclear magnetic resonance logs, has overcome the shortcoming that well logging method can not show factor of porosity situation of change in the horizontal; Compare with the seismic inversion method, overcome the seismic inversion method and only be applicable to and survey the darker formation porosity of buried depth, and for the relatively poor shortcoming of resolution on the more shallow stratum of buried depth; In addition, this method need not to set up the relation of seismic properties and each well factor of porosity, has therefore avoided earthquake in conjunction with the joint inversion method of the well-log information problem bigger to the error of calculation of formation porosity; Simultaneously, this method need not dangerous operation such as to blow out at detecting area, has avoided the generation of accidents such as casualties.Therefore, the invention solves the defective of existing method of testing, realized fast, effectively surveying under the state in position the factor of porosity that changes big and loose and broken formation on the space on Huo Qu top, coalfield and distributed, have applicability widely.
Embodiment
Embodiment 1: western certain the fire district, coalfield of China, fire district area 100898m
2The equal gross thickness 23m in coal seam, 26 ° at inclination angle, because of coal combustion (fiery district center temperature reaches 1200 ℃), cause face of land outcrop of coal seam and upper formation bulk settling, formed the loose crushing band, be the air feed of understanding fire district, coalfield combustion system, distribution and the migration rule of firing-fighting medium in the scall of Huo Qu top of air-out passage, utilize the boring radar that multiple offset distance has been carried out in the distribution of this fire district upper formation factor of porosity and stride hole transmission detection, detection steps is as follows:
A, arrange many surveys line at the ground surface on Huo Qu top, coalfield, the spacing of many surveys line is 50m; At interval a plurality of measuring points are set on each survey line, the spacing of a plurality of measuring points is 40 ~ 50m; One by one to the measuring point construction drill on the every survey line, on same survey line, in adjacent two borings emitting antenna and receiving antenna are set respectively one by one then, successively the radar multiple offset distance of holing of all borings on the every survey line is striden the hole transmission and surveyed by emitting antenna and receiving antenna, be captured in the travel-time of electromagnetic wave from the emitting antenna to the receiving antenna under the different offset distance situations;
Survey the dipole antenna that used boring radar antenna adopts the 60MHz frequency.During detection, it is motionless at first emitting antenna to be fixed on the aperture, and receiving antenna moves down from the position, aperture, and every 3m once sampling is at the bottom of the hole; Afterwards emitting antenna is moved down 3m, receiving antenna is mobile every 3m once sampling from bottom to top; Carry out so repeatedly, at the bottom of emitting antenna moves to the hole, finish boring radar multiple offset distance and stride hole transmission detection.
B, adopt method of conjugate gradient to carry out iterative by computing machine electromagnetic wave travel-time of collecting, obtain on the every survey line slowness of each pixel on the two dimensional cross-section between adjacent two borings
, according to formula:
Calculate on the every survey line bulk dielectric constant of each pixel on the two dimensional cross-section between adjacent two borings
In the formula:
Be electromagnetic wave velocity of propagation in a vacuum;
The not burnt rock stone sample of c, fire district, collection coalfield upper formation is with the specific inductive capacity of decameter measurement rock sample relative vacuum
D, utilize formula:
Calculate the factor of porosity of each pixel on the two dimensional cross-section between two whenever adjacent borings
By computing machine with on each bar survey line between the adjacent boring factor of porosity data on the two dimensional cross-section splice merging respectively, and make the isogram that every formation porosity on the survey line section distributes, thereby draw the distribution of fire district, coalfield upper formation factor of porosity, the more big then stratum of factor of porosity is more loose, and the more little then stratum of factor of porosity is more closely knit; In the formula:
Be the specific inductive capacity of air,
It is a factor relevant with direction of an electric field.
Get the factor of porosity distribution situation of coalfield fire district upper formation clear by above-mentioned detection, for the air feed of studying fire district combustion system, air-out passage and in the fire extinguishing work progress migration rule of firing-fighting medium in fire district upper formation the basis is provided.
Claims (1)
1. the detection method of coalfield fire district upper formation factor of porosity distribution is characterized in that comprising the steps:
A, arrange many surveys line at the ground surface on Huo Qu top, coalfield, the spacing of many surveys line is 50m; At interval a plurality of measuring points are set on each survey line, the spacing of a plurality of measuring points is 40 ~ 50m; One by one to the measuring point construction drill on the every survey line, on same survey line, in adjacent two borings emitting antenna and receiving antenna are set respectively one by one then, successively the radar multiple offset distance of holing of all borings on the every survey line is striden the hole transmission and surveyed by emitting antenna and receiving antenna, be captured in the travel-time of electromagnetic wave from the emitting antenna to the receiving antenna under the different offset distance situations;
B, with electromagnetic wave travel-time of collecting by Computer Processing, obtain on the every survey line slowness of each pixel on the two dimensional cross-section between adjacent two borings
, according to formula:
Calculate on the every survey line bulk dielectric constant of each pixel on the two dimensional cross-section between adjacent two borings
In the formula:
Be electromagnetic wave velocity of propagation in a vacuum;
The not burnt rock stone sample of c, fire district, collection coalfield upper formation is with the specific inductive capacity of decameter measurement rock sample relative vacuum
D, utilize formula:
Calculate the factor of porosity of each pixel on the two dimensional cross-section between two whenever adjacent borings
By computing machine with on each bar survey line between adjacent two borings the factor of porosity data on the two dimensional cross-section splice merging respectively, and make the isogram that every formation porosity on the survey line section distributes, thereby draw the distribution of fire district, coalfield upper formation factor of porosity, the more big then stratum of factor of porosity is more loose, and the more little then stratum of factor of porosity is more closely knit; In the formula:
Be the specific inductive capacity of air,
It is a factor relevant with direction of an electric field.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110335168 CN102508310B (en) | 2011-10-31 | 2011-10-31 | Detection method for porosity distribution of upper formation of fire district of coal field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110335168 CN102508310B (en) | 2011-10-31 | 2011-10-31 | Detection method for porosity distribution of upper formation of fire district of coal field |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508310A CN102508310A (en) | 2012-06-20 |
CN102508310B true CN102508310B (en) | 2013-08-28 |
Family
ID=46220415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110335168 Expired - Fee Related CN102508310B (en) | 2011-10-31 | 2011-10-31 | Detection method for porosity distribution of upper formation of fire district of coal field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508310B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005706A (en) * | 2015-07-31 | 2015-10-28 | 陕西煤业化工技术研究院有限责任公司 | Method for estimating boundary of burning area of single coal seam |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105628574B (en) * | 2014-10-29 | 2019-04-09 | 中国石油天然气股份有限公司 | The measurement method of rock sample porosity |
CN105352996B (en) * | 2015-10-20 | 2018-08-10 | 中国矿业大学 | A kind of model test method of test underground coal-field fire overlying strata temperature change |
CN105588870B (en) * | 2016-01-04 | 2018-07-20 | 中国矿业大学 | A kind of coal-field fire detection device and method based on sulfur isotope component ratio |
CN109187643A (en) * | 2018-06-25 | 2019-01-11 | 国网福建省电力有限公司宁德供电公司 | A kind of transmission tower column foot gap measurement device based on wide-angle mensuration |
CN111596372B (en) * | 2019-12-22 | 2022-05-17 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Water-containing porosity and mineralization degree inversion method based on electromagnetic wave measurement system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (en) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | Transient electromagnetical method for reservoir pore space anisotropy |
CN102052070A (en) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | Complex lithology clastic rock stratum constituent logging processing method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8452538B2 (en) * | 2009-08-27 | 2013-05-28 | Conocophillips Company | Petrophysical evaluation of subterranean formations |
US8614573B2 (en) * | 2009-09-23 | 2013-12-24 | Schlumberger Technology Corporation | Estimating porosity and fluid volume |
-
2011
- 2011-10-31 CN CN 201110335168 patent/CN102508310B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (en) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | Transient electromagnetical method for reservoir pore space anisotropy |
CN102052070A (en) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | Complex lithology clastic rock stratum constituent logging processing method |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005706A (en) * | 2015-07-31 | 2015-10-28 | 陕西煤业化工技术研究院有限责任公司 | Method for estimating boundary of burning area of single coal seam |
CN105005706B (en) * | 2015-07-31 | 2018-06-29 | 陕西煤业化工技术研究院有限责任公司 | A kind of single coal bed burns the evaluation method on area boundary |
Also Published As
Publication number | Publication date |
---|---|
CN102508310A (en) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508310B (en) | Detection method for porosity distribution of upper formation of fire district of coal field | |
CN102426384B (en) | Method for detecting underground goaf and karst distribution | |
CN108957521B (en) | Long-distance three-dimensional advanced geological prediction method for tunnel | |
CN107526101A (en) | A kind of collection for obtaining earthquake reflected wave and processing method | |
CN102288990B (en) | Ultra-shallow three-dimensional seismic reflection wave exploration method | |
CN102866417A (en) | Device and method for seismic cross hole computed tomography (CT) detection and tomography of underground cave | |
CN108152854A (en) | A kind of lossless detection method and its application based on fine motion power spectral density | |
CN202837558U (en) | Underground karst cave earthquake cross-hole CT (computer tomography) detection and tomographic imaging device | |
Orfanos et al. | Analysis of different geophysical methods in the detection of an underground opening at a controlled test site | |
CN106772554A (en) | A kind of multi-channel transient surface wave exploration method under MODEL OVER COMPLEX TOPOGRAPHY | |
CN103424769A (en) | Combined multi-wave seismic exploration method in gob | |
CN108490485A (en) | Double track tunnel Microseismic monitoring system and its positioning accuracy appraisal procedure | |
CN102590874A (en) | Method for detecting ground surface crack of upland coal-mining subsidence paddy field | |
CN105093314B (en) | A kind of method for determining microseism focus | |
CN106371137B (en) | Seismic detection device for shallow abnormal body and three-dimensional observation method thereof | |
CN111812707A (en) | Method for detecting boulders in subway line site based on micro-motion and surface wave exploration | |
Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
CN112882125A (en) | Tunnel-landslide mass comprehensive detection and stability analysis method | |
CN108614307B (en) | Inverted arch quality detection method and system by using comprehensive geophysical prospecting technology | |
Long et al. | Locating geothermal resources using seismic exploration in Xian county, China | |
Araffa et al. | Integration of geophysical techniques to detect geotechnical hazards: a case study in Mokattam, Cairo, Egypt | |
CN113050085A (en) | Advanced geological prediction method | |
CN116774276A (en) | Quantum measurement-based seismic wave acquisition method and system | |
Isakova et al. | GPR for mapping fractures for the extraction of facing granite from a quarry: A case study from Republic of Karelia | |
CN104375187A (en) | Method for designing cannon well depths based on microlog and non-seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130828 |