CN106646662B - The prediction technique and device in Gas Outburst region - Google Patents

The prediction technique and device in Gas Outburst region Download PDF

Info

Publication number
CN106646662B
CN106646662B CN201610991422.9A CN201610991422A CN106646662B CN 106646662 B CN106646662 B CN 106646662B CN 201610991422 A CN201610991422 A CN 201610991422A CN 106646662 B CN106646662 B CN 106646662B
Authority
CN
China
Prior art keywords
data
region
space
layer
tomography
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
CN201610991422.9A
Other languages
Chinese (zh)
Other versions
CN106646662A (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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201610991422.9A priority Critical patent/CN106646662B/en
Publication of CN106646662A publication Critical patent/CN106646662A/en
Application granted granted Critical
Publication of CN106646662B publication Critical patent/CN106646662B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

The present invention provides a kind of prediction techniques and device in Gas Outburst region, are related to the technical field that coal field gass protrudes hazard prediction, including:Obtain the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, the layer position data of destination layer position;According to seismic wave imaging data, earthquake prestack inversion data and layer position data determine the first subterranean zone and the second subterranean zone within the scope of pre-set space;Determine that the region that the first subterranean zone and the second subterranean zone overlap is region to be determined, wherein the possibility that geological disaster occurs for region to be determined is higher than other regions;It determines whether region to be determined is Gas Outburst region according to the log data in region to be determined, alleviates the technical problem for when predicting coal mine progress Gas Outburst region, causing default precision poor since prediction mode is more single in the prior art.

Description

The prediction technique and device in Gas Outburst region
Technical field
The present invention relates to coal field gass to protrude hazard prediction technical field, more particularly, to a kind of the pre- of Gas Outburst region Survey method and apparatus.
Background technology
With the increase of our national coal mine quantity, the incidence of gas burst accident is also gradually increasing in coal mine.Watt This protrusion accident refers to increase with seam mining depth, the increase of coal-mine gas content, under crustal stress effect so that weak Coal seam abrupt release have a large amount of gas and coal and caused by a kind of accident caused by geological disaster.Target, Gas Outburst Accident has become one of 3 big major impetus disasters in coal mine, our countries are exactly that gas burst accident number, Yi Ji occurs One of death toll country the most serious in gas burst accident.In this regard, related researcher just prevents gas burst accident A large amount of research has been carried out, for the present Research of current coal and gas prominent, can be mainly divided into a variety of research approaches, has been wrapped It includes:Forecasting research method based on seismic prospecting, the forecasting research method based on electromagnetic radiation, the prediction based on coalbed gas geology are ground Study carefully a variety of research methods such as method and other predictions.It is single-element that research element is used in above-mentioned a variety of research approaches, no Gas burst accident can comprehensively be analyzed from multiple angles.That is, what above-mentioned single prediction technique was reflected Outburst danger factor is relatively limited, causes to consider deficiency to the otherness of different coal field geology conditions in actual application, because This, there are some urgently improved places for above-mentioned a variety of research methods.
For when carrying out Gas Outburst region to coal mine and predicting, being led in the prior art since prediction mode is more single It causes to preset the poor technical problem of precision, currently no effective solution has been proposed.
Invention content
The purpose of the present invention is to provide a kind of prediction techniques and device in Gas Outburst region, to alleviate in the prior art When predicting coal mine progress Gas Outburst region, the technology for causing default precision poor since prediction mode is more single is asked Topic.
One side according to the ... of the embodiment of the present invention provides a kind of prediction technique in Gas Outburst region, including:It obtains Seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, the layer position data of destination layer position;According to described Seismic wave imaging data, the earthquake prestack inversion data and layer position data determine the first coal within the scope of the pre-set space Layer region and the second subterranean zone, wherein first subterranean zone is the region that coal seam is crushed, in second subterranean zone The hardness of coal body is less than or equal to default hardness;Determine the area that first subterranean zone and second subterranean zone overlap Domain is region to be determined, wherein the possibility that geological disaster occurs for the region to be determined is higher than other regions;It is waited for according to described Determine that the log data in region chooses Gas Outburst region in the region to be determined.
Further, existed according to the seismic imaging data, the earthquake prestack inversion data and layer position data Determine that the first subterranean zone includes within the scope of the pre-set space:It is true according to the seismic imaging data and layer position data The space distribution information in fixed first tomography space;The second tomography is determined according to the seismic imaging data and layer position data The spatial information in space, wherein the tomography scale in first tomography space is more than the tomography in second tomography space Scale;The spatial information of the space distribution information in first tomography space and second tomography space is carried out pair Than, and using the corresponding region of spatial information in second tomography space in target area as first coal Layer region, wherein the centre distance of the target area and first tomography space is less than or equal to pre-determined distance.
Further, the seismic imaging data include reflection wave imaging data, according to the seismic imaging data Determine that the space distribution information in the first tomography space includes with layer position data:According to the reflection wave imaging data, along The layer position data of the destination layer position calculate the relevant volume data in first tomography space, wherein the relevant volume data is Indicate the longitudinal component of the back wave and the similar degrees of data of cross stream component;The similarity value will be calculated and be determined as institute State the space distribution information in the first tomography space.
Further, the seismic imaging data include diffracted wave imaging data, according to the earthquake prestack inversion number Determine that the spatial information in the second tomography space includes according to layer position data:Along the layer position data of the destination layer position Calculate the RMS amplitude attribute of the diffracted wave imaging data;The root mean square amplitude attribute being calculated will be determined as institute State the spatial information in the second tomography space.
Further, existed according to the seismic imaging data, the earthquake prestack inversion data and layer position data Determine that the second subterranean zone includes within the scope of the pre-set space:The layer digit is extracted according to the earthquake prestack inversion data According to the density and Rockmass Shear modulus information of corresponding rock stratum;Pre-set density, Yi Jiti will be less than in the density extracted Less than the region of pre-set shear modulus information as second subterranean zone in the modulus of shearing information got.
Further, the log data includes apparent resistivity, gamma gamma, natural potential and acoustic speed, natural gal Agate, choosing Gas Outburst region in the region to be determined according to the log data in the region to be determined includes:Described Screening meets the data of preset condition in log data, wherein the preset condition is:The apparent resistivity is more than the first threshold Value, the gamma gamma are more than second threshold, and the natural potential is more than third threshold value, and the acoustic speed is less than the 4th threshold Value, the natural gamma are less than the 5th threshold value;It will be corresponding with the data for meeting preset condition in the region to be determined Region as the Gas Outburst region.
Further, the log data is the data that initial log data obtains later by data conversion, wherein institute It includes at least one of to state data conversion:Depth conversion, time conversion.
One side according to the ... of the embodiment of the present invention provides a kind of prediction meanss in Gas Outburst region, including:It obtains Unit, for obtaining the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, the layer digit of destination layer position According to;First determination unit, for according to the seismic imaging data, the earthquake prestack inversion data and layer position data The first subterranean zone and the second subterranean zone are determined within the scope of the pre-set space, wherein first subterranean zone is coal The broken region of layer, the hardness of coal body is less than or equal to default hardness in second subterranean zone;Second determination unit is used In determining that the region that first subterranean zone and second subterranean zone overlap is region to be determined, wherein described to wait for really Determine region and the possibility of geological disaster occurs higher than other regions;Selection unit, for the well logging according to the region to be determined Data choose Gas Outburst region in the region to be determined.
Further, first determination unit includes:First determining module, for according to the seismic imaging data The space distribution information in the first tomography space is determined with layer position data;Second determining module, for according to the seismic wave Imaging data and layer position data determine the spatial information in the second tomography space, wherein first tomography space Tomography scale is more than the tomography scale in second tomography space;Contrast module is used for the space in first tomography space The spatial information in distributed intelligence and second tomography space is compared, and by described second in target area The corresponding region of spatial information in tomography space is as first subterranean zone, wherein the target area with it is described The centre distance in the first tomography space is less than or equal to pre-determined distance.
Further, the seismic imaging data include reflection wave imaging data, and first determining module includes:The One computational submodule calculates described for according to the reflection wave imaging data along the layer position data of the destination layer position The relevant volume data in one tomography space, wherein the relevant volume data is to indicate the longitudinal component of the back wave and laterally divide The similar degrees of data of amount;First determination sub-module is determined as the first tomography sky for the similarity value will to be calculated Between space distribution information.
In the prediction technique in Gas Outburst region provided by the invention, the seismic wave within the scope of pre-set space is obtained first Then the layer position data of imaging data, earthquake prestack inversion data and destination layer position are determined default by above-mentioned three kinds of data Coal seam in spatial dimension is crushed region (that is, first subterranean zone) and the smaller region of seam hardness (that is, the second seam area Domain);Next, determining region to be determined according to the smaller region in the broken region in the coal seam determined and seam hardness;Finally lead to Cross test data judge this it is to be determined whether be Gas Outburst region.Using Gas Outburst region provided in an embodiment of the present invention Prediction technique achieved the purpose that accurately to predict Gas Outburst region, so alleviate in the prior art to coal mine into When row Gas Outburst region is predicted, cause to preset the poor technical problem of precision since prediction mode is more single.
Description of the drawings
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art Embodiment or attached drawing needed to be used in the description of the prior art are briefly described, it should be apparent that, in being described below Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor It puts, other drawings may also be obtained based on these drawings.
Fig. 1 is a kind of flow chart of the prediction technique in Gas Outburst region according to the ... of the embodiment of the present invention;
Fig. 2 is a kind of flow chart of the first subterranean zone of determination according to the ... of the embodiment of the present invention;
Fig. 3 is a kind of flow chart of the second subterranean zone of determination according to the ... of the embodiment of the present invention;
Fig. 4 is a kind of schematic diagram of the prediction meanss in Gas Outburst region according to the ... of the embodiment of the present invention.
Specific implementation mode
Technical scheme of the present invention is clearly and completely described below in conjunction with attached drawing, it is clear that described implementation Example is a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill The every other embodiment that personnel are obtained without making creative work, shall fall within the protection scope of the present invention.
In the description of the present invention, it should be noted that term "center", "upper", "lower", "left", "right", "vertical", The orientation or positional relationship of the instructions such as "horizontal", "inner", "outside" be based on the orientation or positional relationship shown in the drawings, merely to Convenient for the description present invention and simplify description, do not indicate or imply the indicated device or element must have a particular orientation, With specific azimuth configuration and operation, therefore it is not considered as limiting the invention.In addition, term " first ", " second ", " third " is used for description purposes only, and is not understood to indicate or imply relative importance.
In the description of the present invention, it should be noted that unless otherwise clearly defined and limited, term " installation ", " phase Even ", " connection " shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or be integrally connected;It can Can also be electrical connection to be mechanical connection;It can be directly connected, can also indirectly connected through an intermediary, Ke Yishi Connection inside two elements.For the ordinary skill in the art, above-mentioned term can be understood at this with concrete condition Concrete meaning in invention.
For when carrying out Gas Outburst region to coal mine and predicting, being led in the prior art since prediction mode is more single It causes to preset the poor technical problem of precision, inventor has carried out relevant grind from the angle of geologic origin to Gas Outburst region Study carefully.From geologic origin, Gas Outburst mostly occurs in coal body broken " deformation coal ", wherein if coal body destroys journey Degree is more serious, and the intensity of coal is lower, then the possibility that Gas Outburst occurs for the region is bigger, then the danger in the region Also bigger.Therefore, in embodiments of the present invention, for the geologic feature in Gas Outburst region, by a variety of seismic datas Attribute is analyzed, and then draws a circle to approve out coal body fault development band;Then, joined using the earthquake prestack inversion in a variety of seismic datas Number marks off coal body hypo-intense region;And using fault development band and coal body hypo-intense region overlapping region as Gas Outburst sensitivity Region;It is finally demarcated by log data parameter index, predicts Gas Outburst range.
Fig. 1 is a kind of flow chart of the prediction technique in Gas Outburst region according to the ... of the embodiment of the present invention, as shown in Figure 1, This method comprises the following steps:
Step S102 obtains the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, destination layer position Layer position data.
In embodiments of the present invention, the underground space range that pre-set space ranging from related technical personnel choose in advance, In, the horizontal spread and buried depth of the underground space range are that related technical personnel are pre-set.In embodiments of the present invention, it obtains The seismic wave got includes diffracted wave and back wave, and therefore, seismic imaging data include diffracted wave imaging data and back wave Imaging data;Earthquake prestack inversion data are the data to being obtained after seismic data process using prestack inversion algorithm;Layer position Data are the layer bit space location coordinate information of destination layer position, wherein since rock stratum has multiple layers of position, in the present invention In embodiment, in order to simplify calculating process, destination layer position is chosen in multiple layers of position of rock stratum in advance, wherein the destination layer position For the maximum layer position of possibility of tomography occurs in multiple layers of position.
Step S104, according to seismic wave imaging data, earthquake prestack inversion data and layer position data are in pre-set space range The first subterranean zone of interior determination and the second subterranean zone, wherein the first subterranean zone is the region that coal seam is crushed, the second seam area The hardness of coal body is less than or equal to default hardness in domain.
In embodiments of the present invention, get seismic imaging data, earthquake prestack inversion data and layer position data it Afterwards, so that it may to determine the first subterranean zone and the second subterranean zone within the scope of pre-set space according to above-mentioned three kinds of data.Specifically Ground, the first subterranean zone are the region that coal seam is crushed, and the second subterranean zone is the smaller region of hardness, that is, hardness be less than or Person is equal to the region of default hardness.After determining the first subterranean zone and the second subterranean zone, so that it may with according to the first coal seam Region and the second subterranean zone determine Gas Outburst sensitizing range (that is, region to be determined in step S106), and Gas Outburst is quick Sensillary area domain shows that the possibility of region generation Gas Outburst is higher.In embodiments of the present invention, the first subterranean zone and second Subterranean zone may include the region overlapped, can not also include the region overlapped.
Therefore, in embodiments of the present invention, by the smaller region in the broken region in comprehensive analysis coal seam and seam hardness, and The mode that Gas Outburst sensitizing range is determined according to analysis result determines Gas Outburst sensitizing range by various analyses Domain avoids in turn and determines precision of prediction caused by Gas Outburst sensitizing range only on one side in the prior art Poor problem.Further, in embodiments of the present invention, when determining the first subterranean zone and the second subterranean zone, equally A plurality of types of data (that is, seismic imaging data, earthquake prestack inversion data, layer position data of destination layer position) are used, Can equally avoid in the prior art the data only with one side predicted caused by Gas Outburst sensitizing range to determine The poor problem of precision.
Step S106 determines that the region that the first subterranean zone and the second subterranean zone overlap is region to be determined, wherein wait for Determine that the possibility of geological disaster occurs for region higher than other regions.
By foregoing description it is found that in the present invention is implemented, passing through the first subterranean zone of comprehensive analysis and the second seam area Domain determines Gas Outburst sensitizing range (that is, region to be determined).Due to true according to the first subterranean zone and the second subterranean zone There are many kinds of the modes for determining Gas Outburst sensitizing range, in embodiments of the present invention, preferably the first subterranean zone and the second coal seam The region overlapped in region is region to be determined.That is, if some region, which is coal seam, is crushed region, and the coal seam is broken The hardness in coal seam is smaller in broken region, then the possibility that Gas Outburst occurs for the region is larger, at this point, by the region (that is, weight Close region) it is used as Gas Outburst sensitizing range.
Step S108 chooses Gas Outburst region according to the log data in region to be determined in region to be determined.
In embodiments of the present invention, after determining region to be determined, it can not directly determine that the region is Gas Outburst Region, it is also necessary to the region further be tested, specifically, log data may be used to region to be determined into traveling Whether one pacing tries, to include Gas Outburst region in determination region to be determined.
Log data is that can be obtained the best geologic data of resolution ratio highest, continuity, log data so far In include abundant geological information, different test datas records the history of this Geological Evolution to varying degrees, never homonymy Face reflect landing surface mineralizing evolution condition and influence factor, as sea level variability, palaeoenvironment, ancient geoaraply, Palaeoclimatic Information and Its situation of change etc..Therefore, in embodiments of the present invention, region to be determined is carried out using log data further determining energy Enough improve the precision of prediction in Gas Outburst region.
When the raising of the precision of prediction in Gas Outburst region can save the prediction of related technical personnel to a certain extent Between, to improve the manufacturing schedule of coal mine, simultaneously, additionally it is possible to ensure the personal safety of coal mine working personnel.
In the prediction technique in Gas Outburst region provided by the invention, the seismic wave within the scope of pre-set space is obtained first Then the layer position data of imaging data, earthquake prestack inversion data and destination layer position are determined default by above-mentioned three kinds of data Coal seam in spatial dimension is crushed region (that is, first subterranean zone) and the smaller region of seam hardness (that is, the second seam area Domain);Next, determining region to be determined according to the smaller region in the broken region in the coal seam determined and seam hardness;Finally lead to Cross test data judge this it is to be determined whether be Gas Outburst region.Using Gas Outburst region provided in an embodiment of the present invention Prediction technique achieved the purpose that accurately to predict Gas Outburst region, so alleviate in the prior art to coal mine into When row Gas Outburst region is predicted, cause to preset the poor technical problem of precision since prediction mode is more single.
Fig. 2 is a kind of flow chart of the first subterranean zone of determination according to the ... of the embodiment of the present invention, as shown in Fig. 2, specifically, According to seismic wave imaging data, earthquake prestack inversion data and layer position data determine the first seam area within the scope of pre-set space Domain includes the following steps S201 to step S203:
Step S201 determines the space distribution information in the first tomography space according to seismic wave imaging data and layer position data;
Step S202 determines the spatial information in the second tomography space according to seismic wave imaging data and layer position data, Wherein, the tomography scale in the first tomography space is more than the tomography scale in the second tomography space;
Step S203 carries out the spatial information in the space distribution information in the first tomography space and the second tomography space Comparison, and using the corresponding region of spatial information in the second tomography space in target area as the first seam area Domain, wherein the centre distance of target area and the first tomography space is less than or equal to pre-determined distance.
In embodiments of the present invention, when determining that coal seam is crushed region, it is possible, firstly, to according to getting in step S102 Seismic imaging data and the layer position data of destination layer position determine the space distribution information in large scale tomography space, wherein big ruler It is above-mentioned first tomography space to spend tomography space.Then, determine that small scale is disconnected according to earthquake prestack inversion data and layer position data The spatial information of sheaf space (that is, second tomography space), wherein the tomography scale in large scale tomography space is more than small scale The tomography scale in tomography space.
The space distribution information and small scale tomography space of determining large scale tomography space spatial information it Afterwards, so that it may the space distribution information in the first tomography space and the spatial information in the second tomography space to be compared, and The spatial information in the second tomography space in target area is corresponding as the first subterranean zone.Of the invention real It applies in example, above-mentioned target area is the peripheral region of space distribution information, it is understood that be that target area is that the first tomography is empty Between peripheral region, it is also understood that centre distance between target area and the first tomography space be less than or equal to it is default Distance, wherein the size of pre-determined distance can be set according to actual needs.
In embodiments of the present invention, in the space for the space distribution information and the second tomography space for knowing the first tomography space After spread information, so that it may the corresponding region of spatial information of the first tomography spatial peripheral will be in as the first coal seam Region.So far, the broken region (that is, first subterranean zone) in coal seam has determined that, by above-mentioned determining method it is found that In the embodiment of the present invention, during determining the first subterranean zone, using 2 aspects (that is, space distribution information and space exhibition Cloth information) be determined, compared with the existing technology in determine the broken region in coal seam using on one side by the way of, Neng Gougeng Add and accurately determines the broken region in coal seam.After accurately determining coal output layer and being crushed region, can also more accurately it determine Go out Gas Outburst region.
In the optional embodiment of the present invention, if seismic imaging data include reflection wave imaging data, In this case, the space distribution information in the first tomography space is determined according to seismic wave imaging data and layer position data, including it is as follows Step:
Step S2011 calculates the first tomography space according to reflection wave imaging data along the layer position data of destination layer position Relevant volume data, wherein relevant volume data is the similar degrees of data of the longitudinal component and cross stream component that indicate back wave;
Similarity value will be calculated be determined as the space distribution information in the first tomography space in step S2012.
In embodiments of the present invention, can be according to the reflection wave imaging data in seismic wave imaging data, and utilize relevant Technology seeks the relevant volume data in the first tomography space along layer position, that is, seeks the phase about large scale tomography space distribution rule Dry attribute.Specifically, it is that back wave is calculated in layer plane by coherent technique when seeking coherence properties using coherent technique Local waveform similarity on vertical and horizontal, obtains the estimated value of seismic data correlation, in turn, similar by what is be calculated Property value is determined as space distribution information.
Include the feelings of diffracted wave imaging data in seismic imaging data in another optional embodiment of the present invention Under condition, determine that the spatial information in the second tomography space includes following step according to earthquake prestack inversion data and layer position data Suddenly:
Step S2021 calculates the RMS amplitude attribute of diffracted wave imaging data along the layer position data of destination layer position;
The root mean square amplitude attribute being calculated is determined as the spatial information in the second tomography space by step S2022.
It in embodiments of the present invention, can be with seismic imaging number when determining the spatial information in the second tomography space Diffracted wave imaging data in, and RMS amplitude attribute is sought along layer position data, wherein RMS amplitude attribute refers to The average value of Amplitude-squared extracts square root again.Above-mentioned RMS amplitude attribute is calculated to be also understood that as along the diffraction of destination layer position Face calculates the RMS amplitude attribute of diffracted wave.Wherein, the RMS amplitude attribute is for detecting strong shock preparation value, and is somebody's turn to do RMS amplitude attribute corresponds to small scale and is broken spatial information.
By foregoing description it is found that in embodiments of the present invention, according to reflection wave imaging data, and utilizing coherent technique edge It layer position data and seeks coherence properties about large scale tomography space distribution rule;And according to diffracted wave imaging data, edge It layer position data and seeks RMS amplitude attribute, with the spatial information in determination small scale tomography space, wherein root mean square shakes Width attribute refers to that the average value of Amplitude-squared extracts square root again;Next, comparing large scale tomography space distribution information and small scale It is broken spatial information, draws a circle to approve out the small scale region of fracture being distributed around large scale tomography, and using the region as coal The broken region of layer.By the calculating process of above-mentioned data it is found that during determining that coal seam is crushed region, seismic wave is used The layer position data of imaging data and destination layer position, that is, the present invention obtains coal seam fracture area by handling multidimensional data Domain carries out the determination that coal seam is crushed region relative to only with seismic imaging data, or carries out coal only with layer position data The mode of the determination in the broken region of layer improves the determination precision that coal seam is crushed region.
It should be noted that when calculating RMS amplitude attribute, window when can be the setting of layer plane in advance, then, Calculate RMS amplitude attribute during according to it is pre-set when window carry out RMS amplitude attribute calculating.
Fig. 3 is a kind of flow chart of the second subterranean zone of determination according to the ... of the embodiment of the present invention, as shown in figure 3, specifically, According to seismic wave imaging data, earthquake prestack inversion data and layer position data determine the second seam area within the scope of pre-set space Domain includes the following steps S301 and step S302:
Step S301, according to the density and Rockmass Shear modulus of data corresponding rock stratum in inverting data extraction layer position before earthquake Information;
Step S302 will be less than pre-set density, and be less than in the modulus of shearing information extracted in the density extracted The region of pre-set shear modulus information is as the second subterranean zone.
In embodiments of the present invention, after determining that coal seam is crushed region, it is also necessary to determine the area that hardness is smaller in coal seam Domain determines the second subterranean zone.During determining the second subterranean zone, the earthquake prestack inversion data of use and layer position Data in turn, determine the second subterranean zone that is, the determination process is handled multidimensional data.
Since earthquake prestack inversion data include density and modulus of shearing inverting information, determining the second coal seam When region, first, the density and modulus of shearing information of the corresponding rock stratum in destination layer position are extracted by earthquake prestack inversion data. After extracting stratum density and modulus of shearing information, the density extracted can be compared with pre-set density, and will The modulus of shearing information extracted is compared with pre-set shear modulus information, determines that density is less than pre-set density, and shearing Modulus information be less than pre-set shear modulus information corresponding to region as the second subterranean zone, i.e., it is smaller as seam hardness Region.
In the optional embodiment of the present invention, if log data includes apparent resistivity, gamma gamma, naturally electricity Position and acoustic speed, natural gamma, then choosing Gas Outburst in region to be determined according to the log data in region to be determined Region includes the following steps:
Step S1 meets the data of preset condition in log data screening, wherein preset condition is that apparent resistivity is more than the One threshold value, gamma gamma are more than second threshold, and natural potential is more than third threshold value, and acoustic speed is less than the 4th threshold value, natural gal Agate is less than the 5th threshold value;
Step S2, using region corresponding with the data for meeting preset condition in region to be determined as Gas Outburst area Domain.
In embodiments of the present invention, determine region to be determined (that is, Gas Outburst sensitizing range in above-mentioned steps S106 Domain) after, it needs to analyze the log data in region to be determined.Test data is met to the area corresponding to the data of preset condition Domain is as Gas Outburst region.
In embodiments of the present invention, log data be chosen for apparent resistivity, gamma gamma, natural potential and acoustic speed, Natural gamma is respectively provided with pair for above-mentioned log data in advance before screening log data meets the data of preset condition The threshold value answered.For example, be apparent resistivity setting threshold value be first threshold, be gamma gamma be arranged threshold value be second threshold, For natural potential setting threshold value be third threshold value, be acoustic speed setting threshold value be the 4th threshold value, be natural gamma setting Threshold value be the 5th threshold value.
During screening, mainly screening apparent resistivity and be more than first threshold, gamma gamma is more than second threshold, from Right current potential is more than third threshold value, and acoustic speed is less than the 4th threshold value and natural gamma is less than the data of the 4th threshold value.If sieve It selects apparent resistivity and is more than first threshold, gamma gamma is more than second threshold, and natural potential is more than third threshold value, and acoustic speed is small In the 4th threshold value, natural gamma is less than the data of the 5th threshold value, then corresponding region is made in region to be determined by the data For Gas Outburst region.That is, the region for meeting " three high two is low " in above-mentioned log data is Gas Outburst region, In, " three high two is low " refers to high apparent resistivity, high gamma gamma, high natural potential and wave velocity, low natural gamma in a low voice.
By the above-mentioned description to log data it is found that log data can reflect the geological structure origin cause of formation in the region, because This chooses Gas Outburst region by log data in region to be determined, is to region to be determined in embodiments of the present invention The geological structure origin cause of formation analyzed, in turn, Gas Outburst region is determined according to analysis result.
The present invention an optional embodiment in, above-mentioned log data be initial log data by data conversion it The data obtained afterwards, wherein transformation includes at least one of:Depth conversion, time conversion.
It should be noted that during the format of test data has very much, in embodiments of the present invention, the log data of selection is The log data that initial log data obtains after time conversion and depth conversion, after time change and depth conversion The log data arrived and seismic imaging data, earthquake prestack inversion data and layer position data match.
The prediction technique in Gas Outburst region provided in an embodiment of the present invention, according to Gas Outburst position on physical mechanism The construction and rock physical property of coal body;On geophysical prospecting method, by using seismic attributes analysis, prestack inversion number According to and log data, consider various geophysics indexs, in practical applications not so as to avoid Individual forecast method Foot place.The prediction technique in Gas Outburst region provided in an embodiment of the present invention can make early warning to Gas Outburst region, to keep away Exempt from that correspondingly geological disaster occurs during coal mining, for safety of coal mines, efficiently production provides safeguard.
The embodiment of the present invention additionally provides a kind of prediction meanss in Gas Outburst region, the prediction dress in the Gas Outburst region The prediction technique for being mainly used for executing the Gas Outburst region that the above of the embodiment of the present invention is provided is set, below to the present invention The prediction meanss in the Gas Outburst region that embodiment provides do specific introduction.
Fig. 4 is a kind of schematic diagram of the prediction in Gas Outburst region according to the ... of the embodiment of the present invention, as shown in figure 4, this watt The prediction meanss of this outburst area include mainly acquiring unit 41, the first determination unit 43, the second determination unit 45 and choose single Member 47, wherein:
Acquiring unit 41, for obtaining the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, mesh Mark the layer position data of layer position.
In embodiments of the present invention, the underground space range that pre-set space ranging from related technical personnel choose in advance, In, the horizontal spread and buried depth of the underground space range are that related technical personnel are pre-set.In embodiments of the present invention, ground Seismic wave includes diffracted wave and back wave, and therefore, seismic imaging data include diffracted wave imaging data and reflection wave imaging data; Earthquake prestack inversion data are the data to being obtained after seismic data process using prestack inversion algorithm;Layer position data are target The layer bit space location coordinate information of layer position, wherein since rock stratum has multiple layers of position, in embodiments of the present invention, it is Simplified calculating process chooses destination layer position in the multiple layers of position of rock stratum in advance, wherein the destination layer position is in multiple layers of position The maximum layer position of possibility of tomography occurs.
First determination unit 43, for according to seismic wave imaging data, earthquake prestack inversion data and layer position data to be pre- If determining the first subterranean zone and the second subterranean zone in spatial dimension, wherein the first subterranean zone is the region that coal seam is crushed, The hardness of coal body is less than or equal to default hardness in second subterranean zone.
In embodiments of the present invention, get seismic imaging data, earthquake prestack inversion data and layer position data it Afterwards, so that it may to determine the first subterranean zone and the second subterranean zone within the scope of pre-set space according to above-mentioned three kinds of data.Specifically Ground, the first subterranean zone are the region that coal seam is crushed, and the second subterranean zone is the smaller region of hardness, that is, hardness be less than or Person is equal to the region of default hardness.After determining the first subterranean zone and the second subterranean zone, so that it may with according to the first coal seam Region and the second subterranean zone determine Gas Outburst sensitizing range (that is, region to be determined in step S106), and Gas Outburst is quick Sensillary area domain shows that the possibility of region generation Gas Outburst is higher.In embodiments of the present invention, the first subterranean zone and second Subterranean zone may include the region overlapped, can not also include the region overlapped.
Therefore, in embodiments of the present invention, by the smaller region in the broken region in comprehensive analysis coal seam and seam hardness, and The mode that Gas Outburst sensitizing range is determined according to analysis result determines Gas Outburst sensitizing range by various analyses Domain avoids in turn and determines precision of prediction caused by Gas Outburst sensitizing range only on one side in the prior art Poor problem.Further, in embodiments of the present invention, when determining the first subterranean zone and the second subterranean zone, equally A plurality of types of data (that is, seismic imaging data, earthquake prestack inversion data, layer position data of destination layer position) are used, Can equally avoid in the prior art the data only with one side predicted caused by Gas Outburst sensitizing range to determine The poor problem of precision.
Second determination unit 45, for determining that the region that the first subterranean zone and the second subterranean zone overlap is area to be determined Domain, wherein the possibility that geological disaster occurs for region to be determined is higher than other regions.
By foregoing description it is found that in the present invention is implemented, passing through the first subterranean zone of comprehensive analysis and the second seam area Domain determines Gas Outburst sensitizing range (that is, region to be determined).Due to true according to the first subterranean zone and the second subterranean zone Determine Gas Outburst sensitizing range mode have it is many in, in embodiments of the present invention, preferably the first subterranean zone and the second coal seam The region overlapped in region is region to be determined.That is, if some region, which is coal seam, is crushed region, and the coal seam is broken The hardness in coal seam is smaller in broken region, then the possibility that Gas Outburst occurs for the region is larger, at this point, by the region (that is, weight Close region) it is used as Gas Outburst sensitizing range.
Selection unit 47, for choosing Gas Outburst area in region to be determined according to the log data in region to be determined Domain.
In embodiments of the present invention, after determining region to be determined, it can not directly determine that the region is Gas Outburst Region, it is also necessary to the region further be tested, specifically, log data may be used to region to be determined into traveling Whether one pacing tries, be Gas Outburst region with determination region to be determined.
Log data is that can be obtained the best geologic data of resolution ratio highest, continuity, log data so far In include abundant geological information, different test datas records the history of this Geological Evolution to varying degrees, never homonymy Face reflect landing surface mineralizing evolution condition and influence factor, as sea level variability, palaeoenvironment, ancient geoaraply, Palaeoclimatic Information and Its situation of change etc..Therefore, in embodiments of the present invention, region to be determined is carried out using log data further determining energy Enough improve the precision of prediction in Gas Outburst region.
When the raising of the precision of prediction in Gas Outburst region can save the prediction of related technical personnel to a certain extent Between, to improve the manufacturing schedule of coal mine, simultaneously, additionally it is possible to ensure the personal safety of coal mine working personnel.
In the prediction technique in Gas Outburst region provided by the invention, the seismic wave within the scope of pre-set space is obtained first Then the layer position data of imaging data, earthquake prestack inversion data and destination layer position are determined default by above-mentioned three kinds of data Coal seam in spatial dimension is crushed region (that is, first subterranean zone) and the smaller region of seam hardness (that is, the second seam area Domain);Next, determining region to be determined according to the smaller region in the broken region in the coal seam determined and seam hardness;Finally lead to Cross test data judge this it is to be determined whether be Gas Outburst region.Using Gas Outburst region provided in an embodiment of the present invention Prediction technique achieved the purpose that accurately to predict Gas Outburst region, so alleviate in the prior art to coal mine into When row Gas Outburst region is predicted, cause to preset the poor technical problem of precision since prediction mode is more single.
Optionally, the first determination unit includes:First determining module, for according to seismic wave imaging data and layer position data Determine the space distribution information in the first tomography space;Second determining module, for according to seismic wave imaging data and layer position data Determine the spatial information in the second tomography space, wherein the tomography scale in the first tomography space is more than the second tomography space Tomography scale;Contrast module, for believing the space distribution information in the first tomography space and the spatial in the second tomography space Breath is compared, and using the corresponding region of spatial information in the second tomography space in target area as the first coal Layer region, wherein the centre distance of target area and the first tomography space is less than or equal to pre-determined distance.
Optionally, seismic imaging data include reflection wave imaging data, and the first determining module includes:First calculates submodule Block, for according to reflection wave imaging data, the relevant volume data in the first tomography space to be calculated along the layer position data of destination layer position, Wherein, be concerned with the similar degrees of data that volume data is the longitudinal component and cross stream component that indicate back wave;First determination sub-module is used In the space distribution information that similarity value will be calculated be determined as the first tomography space.
Optionally, seismic imaging data include diffracted wave imaging data, and the first determining module further includes:Second calculates son Module, the RMS amplitude attribute for calculating diffracted wave imaging data along the layer position data of destination layer position;Second determines son Module, the root mean square amplitude attribute for will be calculated are determined as the spatial information in the second tomography space.
Optionally, the first determination unit further includes:Extraction module, for according to earthquake prestack inversion data extraction layer digit According to the density and Rockmass Shear modulus information of corresponding rock stratum;Third determining module, it is pre- for will be less than in the density extracted If being less than the region of pre-set shear modulus information in density, and the modulus of shearing information extracted as the second subterranean zone.
Optionally, log data includes apparent resistivity, gamma gamma, natural potential and acoustic speed, natural gamma, is chosen Unit includes:Screening module, the data for meeting preset condition in log data screening, wherein preset condition is apparent resistance Rate is more than first threshold, and gamma gamma is more than second threshold, and natural potential is more than third threshold value, and acoustic speed is less than the 4th threshold Value, natural gamma are less than the 5th threshold value;4th determining module, the data for log data screening to be met to preset condition correspond to Region as Gas Outburst region.
Optionally, log data is the data that initial log data obtains later by data conversion, wherein data conversion Including at least one of:Depth conversion, time conversion.
Finally it should be noted that:The above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations;To the greatest extent Present invention has been described in detail with reference to the aforementioned embodiments for pipe, it will be understood by those of ordinary skill in the art that:Its according to So can with technical scheme described in the above embodiments is modified, either to which part or all technical features into Row equivalent replacement;And these modifications or replacements, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution The range of scheme.

Claims (8)

1. a kind of prediction technique in Gas Outburst region, which is characterized in that including:
Obtain the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, the layer position data of destination layer position;
According to the seismic imaging data, the earthquake prestack inversion data and layer position data are in the pre-set space model Enclose the first subterranean zone of interior determination and the second subterranean zone, wherein first subterranean zone is the region that coal seam is crushed, described The hardness of coal body is less than or equal to default hardness in second subterranean zone;
Determine that the region that first subterranean zone and second subterranean zone overlap is region to be determined, wherein described to wait for Determine that the possibility of geological disaster occurs for region higher than other regions;
Gas Outburst region is chosen in the region to be determined according to the log data in the region to be determined;
Wherein, according to the seismic imaging data, the earthquake prestack inversion data and layer position data described default Determine that the first subterranean zone includes in spatial dimension:
The space distribution information in the first tomography space is determined according to the seismic imaging data and layer position data;
The spatial information in the second tomography space is determined according to the seismic imaging data and layer position data, wherein The tomography scale in first tomography space is more than the tomography scale in second tomography space;
The spatial information of the space distribution information in first tomography space and second tomography space is compared, And using the corresponding region of spatial information in second tomography space in target area as first coal seam Region, wherein the centre distance of the target area and first tomography space is less than or equal to pre-determined distance.
2. prediction technique according to claim 1, which is characterized in that the seismic imaging data include reflection wave imaging Data determine that the space distribution information in the first tomography space includes according to the seismic imaging data and layer position data:
According to the reflection wave imaging data, the phase in first tomography space is calculated along the layer position data of the destination layer position Stem body data, wherein the relevant volume data is the similar degrees of data of the longitudinal component and cross stream component that indicate the back wave;
The similarity value will be calculated it is determined as the space distribution information in first tomography space.
3. prediction technique according to claim 1, which is characterized in that the seismic imaging data include diffracted wave imaging Data determine that the spatial information in the second tomography space includes according to the seismic imaging data and layer position data:
The RMS amplitude attribute of the diffracted wave imaging data is calculated along the layer position data of the destination layer position;
The root mean square amplitude attribute being calculated is determined as to the spatial information in second tomography space.
4. prediction technique according to claim 1, which is characterized in that according to the seismic imaging data, the earthquake Prestack inversion data and layer position data determine that the second subterranean zone includes within the scope of the pre-set space:
The density and Rockmass Shear modulus letter of the corresponding rock stratum of layer position data are extracted according to the earthquake prestack inversion data Breath;
Pre-set density will be less than in the density extracted, and is less than default cut in the modulus of shearing information extracted The region of shear modulu information is as second subterranean zone.
5. prediction technique according to any one of claim 1 to 4, which is characterized in that the log data includes regarding electricity Resistance rate, gamma gamma, natural potential and acoustic speed, natural gamma, according to the log data in the region to be determined described Gas Outburst region is chosen in region to be determined includes:
Screening meets the data of preset condition in the log data, wherein the preset condition is:The apparent resistivity is big In first threshold, the gamma gamma is more than second threshold, and the natural potential is more than third threshold value, and the acoustic speed is less than 4th threshold value, the natural gamma are less than the 5th threshold value;
Using region corresponding with the data for meeting preset condition in the region to be determined as the Gas Outburst area Domain.
6. prediction technique according to any one of claim 1 to 4, which is characterized in that the log data is initial surveys The data that well data obtain later by data conversion, wherein the data conversion includes at least one of:Depth conversion, Time converts.
7. a kind of prediction meanss in Gas Outburst region, which is characterized in that including:
Acquiring unit, for obtaining the seismic imaging data within the scope of pre-set space, earthquake prestack inversion data, destination layer position Layer position data;
First determination unit, for according to the seismic imaging data, the earthquake prestack inversion data and the layer digit According to determining first subterranean zone and the second subterranean zone within the scope of the pre-set space, wherein first subterranean zone is The broken region in coal seam, the hardness of coal body is less than or equal to default hardness in second subterranean zone;
Second determination unit, for determining that the region that first subterranean zone and second subterranean zone overlap is to be determined Region, wherein the possibility that geological disaster occurs for the region to be determined is higher than other regions;
Selection unit, for choosing Gas Outburst area in the region to be determined according to the log data in the region to be determined Domain;
Wherein, first determination unit includes:
First determining module, the sky for determining the first tomography space according to the seismic imaging data and layer position data Between distributed intelligence;
Second determining module, the sky for determining the second tomography space according to the seismic imaging data and layer position data Between spread information, wherein the tomography scale in first tomography space be more than second tomography space tomography scale;
Contrast module is used for the spatial of the space distribution information in first tomography space and second tomography space Information is compared, and using the corresponding region of spatial information in second tomography space in target area as First subterranean zone, wherein the centre distance in the target area and first tomography space is less than or equal to pre- If distance.
8. prediction meanss according to claim 7, which is characterized in that the seismic imaging data include reflection wave imaging Data, first determining module include:
First computational submodule, for according to the reflection wave imaging data, being calculated along the layer position data of the destination layer position The relevant volume data in first tomography space, wherein the relevant volume data be indicate the back wave longitudinal component and The similar degrees of data of cross stream component;
First determination sub-module is determined as the spatial distribution in first tomography space for the similarity value will to be calculated Information.
CN201610991422.9A 2016-11-10 2016-11-10 The prediction technique and device in Gas Outburst region Active CN106646662B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610991422.9A CN106646662B (en) 2016-11-10 2016-11-10 The prediction technique and device in Gas Outburst region

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610991422.9A CN106646662B (en) 2016-11-10 2016-11-10 The prediction technique and device in Gas Outburst region

Publications (2)

Publication Number Publication Date
CN106646662A CN106646662A (en) 2017-05-10
CN106646662B true CN106646662B (en) 2018-09-18

Family

ID=58806118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610991422.9A Active CN106646662B (en) 2016-11-10 2016-11-10 The prediction technique and device in Gas Outburst region

Country Status (1)

Country Link
CN (1) CN106646662B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112360562A (en) * 2020-11-11 2021-02-12 中国矿业大学 Method for finely judging and identifying potential of coal and gas outburst dangerous area of driving face

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6618678B1 (en) * 2000-05-26 2003-09-09 Jason Geosystems B.V. Method of joint analysis and interpretation of the subsurface from multiple seismic derived layer property data sets
CN102359924A (en) * 2011-09-19 2012-02-22 中国地质大学(北京) Detection method for coal petrography intensity based on multi-wave seismic data
CN104850897A (en) * 2015-02-25 2015-08-19 中国矿业大学 Prediction method for coal and gas outburst based on seismic information
CN105607121A (en) * 2016-02-02 2016-05-25 中国矿业大学(北京) Coal collapse column identification method and apparatus
CN106019371A (en) * 2016-05-13 2016-10-12 中国矿业大学 Outburst coal seam roadway minor fault advanced qualitative forecast method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6618678B1 (en) * 2000-05-26 2003-09-09 Jason Geosystems B.V. Method of joint analysis and interpretation of the subsurface from multiple seismic derived layer property data sets
CN102359924A (en) * 2011-09-19 2012-02-22 中国地质大学(北京) Detection method for coal petrography intensity based on multi-wave seismic data
CN104850897A (en) * 2015-02-25 2015-08-19 中国矿业大学 Prediction method for coal and gas outburst based on seismic information
CN105607121A (en) * 2016-02-02 2016-05-25 中国矿业大学(北京) Coal collapse column identification method and apparatus
CN106019371A (en) * 2016-05-13 2016-10-12 中国矿业大学 Outburst coal seam roadway minor fault advanced qualitative forecast method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"利用地震AVO反演预测煤与瓦斯突出区";胡朝元 等;《天然气地球科学》;20110831;第22卷(第4期);第728-732页 *

Also Published As

Publication number Publication date
CN106646662A (en) 2017-05-10

Similar Documents

Publication Publication Date Title
Dezayes et al. Structure of the low permeable naturally fractured geothermal reservoir at Soultz
Paul et al. Mapping of cleats and fractures as an indicator of in-situ stress orientation, Jharia coalfield, India
Ilg et al. Normal faults and gas migration in an active plate boundary, southern Taranaki Basin, offshore New Zealand
Malehmir et al. Near-surface geophysical characterization of areas prone to natural hazards: a review of the current and perspective on the future
Stabile et al. A comprehensive approach for evaluating network performance in surface and borehole seismic monitoring
Cianfarra et al. Quantification of fracturing within fault damage zones affecting Late Proterozoic carbonates in Svalbard
Peirce et al. Constraints on crustal structure of adjacent OCCs and segment boundaries at 13° N on the Mid-Atlantic Ridge
CN106646662B (en) The prediction technique and device in Gas Outburst region
Erteleva et al. Methods for assessing the seismic hazard of stable continental areas using combined paleoseismological and geophysical data
Rossland et al. Permian–holocene tectonostratigraphic evolution of the mandal high, central graben, North Sea
Maurer et al. Geophysical characterization of slope instabilities
Martinsen et al. Outcrops revitalized: Tools, techniques and applications
Almalki et al. Efficiency of seismic attributes in detecting near-surface cavities
Isakbek et al. Geophysical investigation of the landslide-prone slope downstream from the Rogun Dam construction site (Tajikistan)
Cakir et al. Shallow seismic site characterizations at 25 ANSS/PNSN stations and compilation of site-specific data for the entire strongmotion network in Washington and Oregon
Chung et al. Seismic site classifications for the St. Louis urban area
Cantore et al. Development of a site-conditions map for the Campania-Lucania region (southern Apennines, Italy)
Carvalho et al. A geophysical study of the Carcavai fault zone, Portugal
Peterie et al. Seismology and its applications in kansas
BLACKWELL et al. PAUL MORGAN
Rucker et al. Reconnaissance evaluation of a potential future sinkhole using integrated simple surface geophysics and surface monitoring points
Brojerdi et al. Reflection seismic imaging of the deeper structures at the Forsmark spent nuclear fuel repository site, central Sweden
Cowell Geophysical Investigation of a Geologic Unconformity at Seward Park-Seattle, WA
Kuehn et al. Automatic determination of full moment tensor solutions from P-wave first motion amplitudes
Gorbatikov et al. Refinement of the deep structure and tectonic movement kinematics in the area of the 2001 salsk earthquake based on new geophysical data

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