CN106646662B - The prediction technique and device in Gas Outburst region - Google Patents
The prediction technique and device in Gas Outburst region Download PDFInfo
- 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
- tomography
- layer position
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting 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
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.
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 (2)
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 |
CN115079285B (en) * | 2022-07-25 | 2022-11-08 | 北京科技大学 | Dynamic prediction method for coal bed gas outburst danger visualization area |
Citations (5)
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 |
-
2016
- 2016-11-10 CN CN201610991422.9A patent/CN106646662B/en active Active
Patent Citations (5)
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)
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 | |
Malehmir et al. | Seismic methods in mineral exploration and mine planning: A general overview of past and present case histories and a look into the future | |
Paul et al. | Mapping of cleats and fractures as an indicator of in-situ stress orientation, Jharia coalfield, India | |
Malehmir et al. | Near-surface geophysical characterization of areas prone to natural hazards: a review of the current and perspective on the future | |
Ilg et al. | Normal faults and gas migration in an active plate boundary, southern Taranaki Basin, offshore New Zealand | |
Stabile et al. | A comprehensive approach for evaluating network performance in surface and borehole seismic monitoring | |
Peirce et al. | Constraints on crustal structure of adjacent OCCs and segment boundaries at 13° N on the Mid-Atlantic Ridge | |
CN112965136A (en) | Multi-stage advanced detection method for water-rich karst tunnel | |
Beeson et al. | The transtensional offshore portion of the northern San Andreas fault: Fault zone geometry, late Pleistocene to Holocene sediment deposition, shallow deformation patterns, and asymmetric basin growth | |
CN106646662B (en) | The prediction technique and device in Gas Outburst region | |
Qin et al. | Study of ground movement in a mining area with geological faults using FDM analysis and a stacking InSAR method | |
Maurer et al. | Geophysical characterization of slope instabilities | |
Green et al. | High-resolution geophysical techniques for improving hazard assessments of unstable rock slopes | |
CN114839671A (en) | Method for finely identifying coal measure stratum electrical interface by ground and ground roadway transient electromagnetic method | |
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 | |
Nkosi et al. | Physical property studies to elucidate the source of seismic reflectivity within the ICDP DSeis seismogenic zone: Klerksdorp goldfield, South Africa | |
Chang et al. | An analysis of the subsurface fault systems with audio-magnetotelluric surveys in the western Ilan Plain of NE Taiwan | |
Shinn et al. | 3D grid based screening process for large-scale CO2 geological storage in Gunsan Basin, Yellow Sea, Korea | |
Cantore et al. | Development of a site-conditions map for the Campania-Lucania region (southern Apennines, Italy) | |
Nkosi et al. | 3D seismic attributes for structural mapping and enhancement of deep gold mining: a case study from the West Wits Line goldfields, South Africa | |
Peterie et al. | Seismology and its applications in kansas | |
Rucker et al. | Reconnaissance evaluation of a potential future sinkhole using integrated simple surface geophysics and surface monitoring points | |
CN106772669A (en) | The detection method and device in Enriching Coalbed Methane region | |
Zhang et al. | Study of seismic diffractions caused by a fracture zone at In Salah carbon dioxide storage project | |
Brojerdi et al. | Reflection seismic imaging of the deeper structures at the Forsmark spent nuclear fuel repository site, central Sweden |
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 |