CN104133246A - Microseism event scanning positioning method and device - Google Patents
Microseism event scanning positioning method and device Download PDFInfo
- Publication number
- CN104133246A CN104133246A CN201410373504.8A CN201410373504A CN104133246A CN 104133246 A CN104133246 A CN 104133246A CN 201410373504 A CN201410373504 A CN 201410373504A CN 104133246 A CN104133246 A CN 104133246A
- Authority
- CN
- China
- Prior art keywords
- event
- grid
- data
- microearthquake
- correction
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a microseism event scanning positioning method and device. The microseism event scanning positioning method includes the steps of collecting simple-component or three-component microseism data; obtaining a check shot, and generating static correction values of geophone points according to the check shot; carrying out time difference correction on the microseism data with the first-break time of the check shot, and respectively carrying out stacking to generate stacked tracks; judging whether the number of arrays, with peak amplitudes larger than a preset amplitude threshold value, in arrays with the preset number reaches a set number or not; if the number of the arrays, with the peak amplitudes larger than the preset amplitude threshold value, in the arrays with the preset number reaches the set number, defining the time when the number of the arrays reaches the set number to be arrival time of one effective event; generating first-break polarity pij by picking up take-off time or extreme value time of the arrays; carrying out static correction on the effective events, multiplying pre-stack amplitudes of the effective events obtained after static correction by the corresponding first-break polarity pij according to event numbers and array numbers corresponding to the effective events to generate corrected data; carrying out seismic focus positioning on the corrected data with a multi-scale mesh energy scanning method.
Description
Technical field
The invention relates to seismic exploration technique, especially in regard to the data processing technique of seismic prospecting, concretely, is about a kind of microearthquake event scanning localization method and device.
Background technology
In unconventional oil gas waterfrac treatment engineering or geothermal utilization engineering, by microearthquake, monitor the acoustie emission event producing rock burst can be obtained time, to the location of these events, can disclose the information such as space distribution that subsurface rock breaks, for the parameter optimization of the engineerings such as waterfrac treatment and recruitment evaluation etc. provide the most directly, help.
In prior art, there are a lot of micro seismic monitoring methods, when not having alternative well to carry out monitoring in well, it is comparatively conventional that monitoring method is arranged on ground, the data that obtain for this method cause first arrival to pick up difficulty because signal to noise ratio (S/N ratio) is lower while carrying out microseism location, major part has adopted energy scan localization method, but general all methods based on blind scanning, the phase place that there is no to consider the radiation waveform that causes due to different focal mechanism changes, therefore for walking sliding situation of breaking, can cause the counteracting of positive-negative phase, thereby reduce the precision of microseism location.
Summary of the invention
The invention provides a kind of microearthquake event scanning localization method and device, to obtain the accurate positioning result of event on the basis in polarity correction and static correction by multiple dimensioned grid energy scan.
To achieve these goals, the invention provides a kind of microearthquake event scanning localization method, described microearthquake event scanning localization method comprises:
By many, ground, arrange and gather simple component or three-component microearthquake data;
Obtain check shot, and according to described check shot, generate the static correction value of each geophone station;
Utilize the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generate stack road;
Judge whether the number that peak amplitude in the arrangement of described predetermined number is greater than the arrangement of a default amplitude threshold reaches a setting number;
Timing definition while if so, the number of described arrangement being reached to described setting number is the time of arrival of a validity event;
Take-off time or the extreme value time of by picking up each, arranging, generate first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event;
Described validity event is carried out to static correction;
The indication using prestack seismic amplitude of the described validity event after static correction is multiplied by corresponding first arrival polarity p by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data;
After utilizing multiple dimensioned grid energy scan method to described correction, data are carried out focus location.
In one embodiment, described in obtain check shot, and according to described check shot, generate the static correction value of each geophone station, comprising:
Utilize perforation or primacord fuse to obtain check shot;
Described check shot is picked up to described first break time, described first break time and initial velocity model are just being drilled to the first break time obtaining and subtracting each other, generate the static correction value of each geophone station.
In one embodiment, described microearthquake event scanning localization method also comprises: the described amplitude threshold that described stack road is set.
In one embodiment, described utilize multiple dimensioned grid energy scan method to described correction after data carry out focus location, comprising:
Subsurface source the first distribution range be take to the first yardstick subdivision as a plurality of the first grids;
While obtaining walking of the first grid described in each by ray-tracing scheme;
During according to the walking of described the first grid, data after described correction are carried out to TEC time error correction stack, generate the first superposition value;
Described first superposition value of each the first grid relatively, grid element center corresponding to described first superposition value of maximum of usining is as Primary Location point;
On described Primary Location point basis, divide the second distribution range, selecting the second yardstick is a plurality of the second grids by described the second distribution range subdivision;
While obtaining walking of the second grid described in each by ray-tracing scheme;
During according to the walking of described the second grid, data after described correction are carried out to TEC time error correction stack, generate the second superposition value;
Relatively described second superposition value of each the second grid, selects grid element center corresponding to maximum described the second superposition value as source point.
In one embodiment, described the first distribution range is greater than described the second distribution range, and described the first yardstick is greater than described the second yardstick.
To achieve these goals, the present invention also provides a kind of microearthquake event scanning locating device, and described microearthquake event scanning locating device comprises:
Data acquisition unit, for being arranged and gathered simple component or three-component microearthquake data by many, ground;
Static correction value generation unit, for obtaining check shot, and generates the static correction value of each geophone station according to described check shot;
Stack road generation unit, for utilizing the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generates stack road;
Judging unit, for judging that the arrangement peak amplitude of described predetermined number is greater than a number of presetting the arrangement of amplitude threshold and whether reaches a setting number;
Time setting unit, the timing definition when the number of described arrangement is reached to described setting number is the time of arrival of a validity event;
First arrival polarity generation unit, for take-off time or the extreme value time of arranging by picking up each, generates first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event;
Static correction unit, for carrying out static correction to described validity event;
Data generating unit after proofreading and correct, for being multiplied by corresponding first arrival polarity p by the indication using prestack seismic amplitude of the described validity event after static correction by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data;
Focus positioning unit, carries out focus location for data after utilizing multiple dimensioned grid energy scan method to described correction.
In one embodiment, described static correction value generation unit comprises:
Check shot acquisition module, for utilizing perforation or primacord fuse to obtain check shot;
Static correction value generation module, for described check shot is picked up to described first break time, is just drilling by described first break time and initial velocity model the first break time obtaining and is subtracting each other, and generates the static correction value of each geophone station.
In one embodiment, described microearthquake event scanning locating device also comprises: amplitude threshold setting unit, and for the described amplitude threshold in described stack road is set.
In one embodiment, described focus positioning unit comprises:
The first mesh generation module, for take subsurface source the first distribution range the first yardstick subdivision as a plurality of the first grids;
The first acquisition module while walking, when obtaining the walking of the first grid described in each by ray-tracing scheme;
The first superposition value generation module, carries out TEC time error correction stack to data after described correction when according to the walking of described the first grid, and generates the first superposition value;
Primary Location point selection module, for comparing described first superposition value of each the first grid, selects grid element center corresponding to maximum described the first superposition value as Primary Location point;
The second mesh generation module, for divide the second distribution range on described Primary Location point basis, second yardstick of take is a plurality of the second grids by described the second distribution range subdivision;
The second acquisition module while walking, when obtaining the walking of the second grid described in each by ray-tracing scheme;
The second superposition value generation module, carries out TEC time error correction stack to data after described correction when according to the walking of described the second grid, and generates the second superposition value;
Source point is selected module, for comparing described second superposition value of each the second grid, selects grid element center corresponding to maximum described the second superposition value as source point.
In one embodiment, described the first distribution range is greater than described the second distribution range, and described the first yardstick is greater than described the second yardstick.
In microearthquake event scanning localization method of the present invention, by check shot, calculate static correction, make bearing calibration simpler, strong adaptability.Utilize check shot first break time to carry out TEC time error correction and detect validity event, fast and easy by arranging stack.By first break picking polarity after arranging stack, be used for polarity correction, the positive-negative polarity energy offset while having avoided energy scan.Multiple dimensioned grid energy scan method has been avoided solving and has been absorbed in local better solution, and calculated amount is relatively little, and precision is high.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, to the accompanying drawing of required use in embodiment or description of the Prior Art be briefly described below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skills, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is the microearthquake event scanning localization method process flow diagram of the embodiment of the present invention;
Fig. 2 is that the embodiment of the present invention generates the method flow diagram of the static correction value of each geophone station according to check shot;
Fig. 3 is the method flow diagram that after the embodiment of the present invention utilizes multiple dimensioned grid energy scan method to described correction, data are carried out focus location;
Fig. 4 is the structured flowchart of the microearthquake event scanning locating device of the embodiment of the present invention;
Fig. 5 is the structured flowchart of the static correction value generation unit 402 of the embodiment of the present invention;
Fig. 6 is the structured flowchart of the focus positioning unit 409 of the embodiment of the present invention;
Fig. 7 is that the ground monitoring of embodiment of the present invention check shot records schematic diagram;
Fig. 8 is embodiment of the present invention static correction value curve synoptic diagram;
Fig. 9 is that the ground monitoring of embodiment of the present invention perforation big gun records result schematic diagram after static correction;
Figure 10 is that 20 records of embodiment of the present invention ground monitoring are arranged the rear result schematic diagram of stack by 12;
Figure 11 is embodiment of the present invention actual ground monitoring microearthquake event typical record schematic diagram;
Figure 12 is result schematic diagram after the static correction of embodiment of the present invention actual ground monitoring microearthquake event typical record;
Figure 13 is result schematic diagram after embodiment of the present invention actual ground monitoring microearthquake event typical record static correction+polarity correction;
Figure 14 is the positioning result vertical view of embodiment of the present invention ground micro seismic monitoring record.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, rather than whole embodiment.Embodiment based in the present invention, those of ordinary skills, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.
As shown in Figure 1, the invention provides a kind of microearthquake event scanning localization method, described microearthquake event scanning localization method comprises:
Step 101: arrange and gather simple component or three-component microearthquake data by many, ground;
Step 102: obtain check shot, and generate the static correction value of each geophone station according to described check shot;
Step 103: utilize the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generate stack road;
Step 104: judge whether the number that peak amplitude in the arrangement of described predetermined number is greater than the arrangement of a default amplitude threshold reaches a setting number;
Step 105: if the timing definition when number of described arrangement is reached to described setting number is the time of arrival of a validity event;
Step 106: take-off time or the extreme value time of arranging by picking up each, generate first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event;
Step 107: described validity event is carried out to static correction;
Step 108: the indication using prestack seismic amplitude of the described validity event after static correction is multiplied by corresponding first arrival polarity p by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data;
Step 109: after utilizing multiple dimensioned grid energy scan method to described correction, data are carried out focus location.
Flow process is as shown in Figure 1 known, first the present invention utilizes the first break time of check shot that the microearthquake data microearthquake data that gather are carried out TEC time error correction and superposeed respectively by the arrangement of predetermined number, when the number that peak amplitude is greater than the arrangement of a default amplitude threshold in the arrangement of the predetermined number of judgement reaches a setting number, the timing definition that the number of arrangement is reached while setting number is the time of arrival of a validity event.Then, take-off time or the extreme value time of by picking up each, arranging, generate first arrival polarity p
ij, validity event is carried out to static correction.Finally, the indication using prestack seismic amplitude of the validity event after static correction is multiplied by corresponding first arrival polarity p by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data, utilize multiple dimensioned grid energy scan method to carry out focus location to proofreading and correct rear data.The present invention is used for polarity correction, the positive-negative polarity energy offset while having avoided energy scan by first break picking polarity after arranging stack; Utilize multiple dimensioned grid energy scan method, avoided solving and be absorbed in local better solution, and calculated amount is relatively little, precision is high.
In one embodiment, as shown in Figure 2, in step 102, obtain check shot, and according to check shot, generate the static correction value of each geophone station, comprising:
Step 201: utilize perforation or primacord fuse to obtain check shot.Obtain the method for check shot with multiple, the invention is not restricted to utilize perforation or primacord fuse to obtain the method for check shot.
Step 202: described check shot is picked up to described first break time, described first break time and initial velocity model are just being drilled to the first break time obtaining and subtracting each other, generate the static correction value of each geophone station.
In one embodiment, before step 104, this microearthquake event scanning localization method also comprises: the amplitude threshold in described stack road is set, i.e. default amplitude threshold in step 104.
During concrete enforcement, as shown in Figure 3, step 109 can comprise:
Step 301: subsurface source the first distribution range be take to the first yardstick subdivision as a plurality of the first grids;
Step 302: while obtaining walking of the first grid described in each by ray-tracing scheme;
Step 303: during according to the walking of described the first grid, data after described correction are carried out to TEC time error correction stack, generate the first superposition value;
Step 304: described first superposition value of each the first grid relatively, grid element center corresponding to described first superposition value of maximum of usining is as Primary Location point;
Step 305: divide the second distribution range on described Primary Location point basis, selecting the second yardstick is a plurality of the second grids by described the second distribution range subdivision;
Step 306: while obtaining walking of the second grid described in each by ray-tracing scheme;
Step 307: during according to the walking of described the second grid, data after described correction are carried out to TEC time error correction stack, generate the second superposition value;
Step 308: compare described second superposition value of each the second grid, select grid element center corresponding to maximum described the second superposition value as source point.
In Fig. 3, to proofreading and correct rear data, carry out focus when location, adopted multiple dimensioned grid energy scan method, avoided solving and be absorbed in local better solution, calculated amount is relatively little, and precision is high.
In one embodiment, the first distribution range shown in Fig. 3 is greater than described the second distribution range, and the first yardstick is greater than described the second yardstick.
As shown in Figure 4, the present invention also provides a kind of microearthquake event scanning locating device, described microearthquake event scanning locating device comprises: data acquisition unit 401, static correction value generation unit 402, stack road generation unit 403, judging unit 404, time setting unit 405, first arrival polarity generation unit 406, static correction unit 407, data generating unit 408 and focus positioning unit 409 after proofreading and correct.
Data acquisition unit 401 can be arranged and be gathered simple component or three-component microearthquake data by many, ground.
Static correction value generation unit 402 can obtain check shot, and according to described check shot, generates the static correction value of each geophone station.
Stack road generation unit 403 can utilize the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generates stack road.
Judging unit 404 can be for judging whether the number that peak amplitude in the arrangement of described predetermined number is greater than the arrangement of a default amplitude threshold reaches a setting number.
Time setting unit 405 can the timing definition when the number of described arrangement is reached to described setting number be the time of arrival of a validity event.
First arrival polarity generation unit 406 can, for take-off time or the extreme value time of arranging by picking up each, generate first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event.
Static correction unit 407 can be for carrying out static correction to described validity event.
After proofreading and correct, data generating unit 408 can be for being multiplied by corresponding first arrival polarity p by the indication using prestack seismic amplitude of the described validity event after static correction by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data.
Focus positioning unit 409 can carry out focus location for data after utilizing multiple dimensioned grid energy scan method to described correction.
Structured flowchart is as shown in Figure 4 known, microearthquake event scanning locating device of the present invention, first utilize the first break time of check shot that the microearthquake data microearthquake data that gather are carried out TEC time error correction and superposeed respectively by the arrangement of predetermined number, when the number that peak amplitude is greater than the arrangement of a default amplitude threshold in the arrangement of the predetermined number of judgement reaches a setting number, the timing definition that the number of arrangement is reached while setting number is the time of arrival of a validity event.Then, take-off time or the extreme value time of by picking up each, arranging, generate first arrival polarity p
ij, validity event is carried out to static correction.Finally, the indication using prestack seismic amplitude of the validity event after static correction is multiplied by corresponding first arrival polarity p by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data, utilize multiple dimensioned grid energy scan method to carry out focus location to proofreading and correct rear data.Microearthquake event scanning locating device of the present invention is used for polarity correction, the positive-negative polarity energy offset while having avoided energy scan by first break picking polarity after arranging stack; Utilize multiple dimensioned grid energy scan method, avoided solving and be absorbed in local better solution, and calculated amount is relatively little, precision is high.
In one embodiment, as shown in Figure 5, static correction value generation unit 402 comprises: check shot acquisition module 501 and static correction value generation module 502.
Check shot acquisition module 501 can be for utilizing perforation or primacord fuse to obtain check shot.Static correction value generation module 502 can be for described check shot is picked up to described first break time, described first break time and initial velocity model just drilled to the first break time obtaining and subtracting each other, and generates the static correction value of each geophone station.
In one embodiment, described microearthquake event scanning locating device also comprises: amplitude threshold setting unit, and for the described amplitude threshold in described stack road is set.
In one embodiment, as shown in Figure 6, focus positioning unit 409 comprises: the first mesh generation module 601, the first 602, the first superposition value generation modules 603 of acquisition module while walking, Primary Location point selection module 604, the second mesh generation module 605, the second 606, the second superposition value generation modules 607 of acquisition module while walking, source point is selected module 608.
The first mesh generation module 601 is for take subsurface source the first distribution range the first yardstick subdivision as a plurality of the first grids.
The first acquisition module 602 while walking is when obtaining the walking of the first grid described in each by ray-tracing scheme.
The first superposition value generation module 603 carries out TEC time error correction stack to data after described correction when according to the walking of described the first grid, and generates the first superposition value.
Primary Location point selection module 604, for comparing described first superposition value of each the first grid, selects grid element center corresponding to maximum described the first superposition value as Primary Location point.
The second mesh generation module 605 for dividing the second distribution range on described Primary Location point basis, and second yardstick of take is a plurality of the second grids by described the second distribution range subdivision.
The second acquisition module 606 while walking is when obtaining the walking of the second grid described in each by ray-tracing scheme.
The second superposition value generation module 607 carries out TEC time error correction stack to data after described correction when according to the walking of described the second grid, and generates the second superposition value.
Source point selects module 608 for comparing described second superposition value of each the second grid, selects grid element center corresponding to maximum described the second superposition value as source point.
In one embodiment, the first distribution range is greater than the second distribution range, and the first yardstick is greater than described the second yardstick.
In order to understand better the present invention, below in conjunction with concrete example, describe the present invention in detail.
The present invention implements at Chinese middle and south block, the radial survey line of having disposed 12 directions for fractured well gathers, be total to image data 500GB, by microearthquake event scanning localization method of the present invention, process and obtained rationally positioning result reliably, feasibility and the practicality of method has been described.Details are as follows for specific implementation step:
1) by 12, ground, arrange 1771 geophone stations and receive, gather simple component microearthquake data.
2) utilize perforation to obtain check shot, the ground monitoring record of this check shot as shown in Figure 7.By by this check shot first break picking, just drilling with the initial velocity model of having set up the first break time obtaining and subtracting each other, obtain the static correction value dt of each geophone station, correcting energy curve is as shown in Figure 8.The first most hyperbola distribution of check shot after static correction, as shown in Figure 9.
3) to all microearthquake data, utilize check shot first break time carry out TEC time error correction and superpose respectively by 12 arrangements, form stack road, Figure 10 is that 20 records of ground monitoring are arranged the rear result schematic diagram of stack by 12.
4) stack road is arranged to amplitude threshold A0=2000, when the peak amplitude of 8 arrangements is greater than threshold value, be made as to the time of arrival of a validity event the corresponding time.
5) because focal mechanism is different, during different arrangement stack, its first arrival phase place can be different, and as shown in figure 11, because the rear signal to noise ratio (S/N ratio) of stack is higher, road number obviously reduces actual ground monitoring microearthquake event typical record.By picking up, can obtain first arrival polarity p
ij, pij represents j polarity of arranging of i event, p here
ij=1 representing positive polarity, p
ij=-1 representing negative polarity, p
ij=0 represents the invalid data that polarity cannot judge.
6) all validity events are carried out to static correction, after the static correction of actual ground monitoring microearthquake event typical record, result as shown in figure 12.
7) on the basis of above-mentioned static correction, the indication using prestack seismic amplitude of all validity events is multiplied by corresponding polarity p by its corresponding event sequence number and arrangement sequence number
ij, after actual ground monitoring microearthquake event typical record static correction+polarity correction, result as shown in figure 13.
8) utilize multiple dimensioned grid energy scan method to carry out focus location.First underground possible focus distribution range is carried out to subdivision (10 meters) with large scale, while asking for walking of each grid by ray tracing, by calculate walk time carry out TEC time error correction stack, the superposition value that each grid is obtained relatively obtains grid element center corresponding to maximal value and is Primary Location point, on this anchor point basis, its less neighborhood is carried out to the more subdivision of small scale (1 meter), scan again stack and obtain grid element center point corresponding to maximal value as source point, finally obtain the positioning result that precision is higher (as shown in figure 14).
Through after above processing procedure, obtained more satisfactory positioning result, microearthquake event scanning localization method of the present invention has taken into full account static correction and first arrival polarity and has implemented corresponding correction, by multiple dimensioned energy scan, has obtained microseismic event positioning result accurately.
In microearthquake event scanning localization method of the present invention, by check shot, calculate static correction, make bearing calibration simpler, strong adaptability.Utilize check shot first break time to carry out TEC time error correction and detect validity event, fast and easy by arranging stack.By first break picking polarity after arranging stack, be used for polarity correction, the positive-negative polarity energy offset while having avoided energy scan.Multiple dimensioned grid energy scan method has been avoided solving and has been absorbed in local better solution, and calculated amount is relatively little, and precision is high.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt complete hardware implementation example, implement software example or in conjunction with the form of the embodiment of software and hardware aspect completely.And the present invention can adopt the form that wherein includes the upper computer program of implementing of computer-usable storage medium (including but not limited to magnetic disk memory, CD-ROM, optical memory etc.) of computer usable program code one or more.
The present invention is with reference to describing according to process flow diagram and/or the block scheme of the method for the embodiment of the present invention, equipment (system) and computer program.Should understand can be in computer program instructions realization flow figure and/or block scheme each flow process and/or the flow process in square frame and process flow diagram and/or block scheme and/or the combination of square frame.Can provide these computer program instructions to the processor of multi-purpose computer, special purpose computer, Embedded Processor or other programmable data processing device to produce a machine, the instruction of carrying out by the processor of computing machine or other programmable data processing device is produced for realizing the device in the function of flow process of process flow diagram or a plurality of flow process and/or square frame of block scheme or a plurality of square frame appointments.
These computer program instructions also can be stored in energy vectoring computer or the computer-readable memory of other programmable data processing device with ad hoc fashion work, the instruction that makes to be stored in this computer-readable memory produces the manufacture that comprises command device, and this command device is realized the function of appointment in flow process of process flow diagram or a plurality of flow process and/or square frame of block scheme or a plurality of square frame.
These computer program instructions also can be loaded in computing machine or other programmable data processing device, make to carry out sequence of operations step to produce computer implemented processing on computing machine or other programmable devices, thereby the instruction of carrying out is provided for realizing the step of the function of appointment in flow process of process flow diagram or a plurality of flow process and/or square frame of block scheme or a plurality of square frame on computing machine or other programmable devices.
In the present invention, applied specific embodiment principle of the present invention and embodiment are set forth, the explanation of above embodiment is just for helping to understand method of the present invention and core concept thereof; , for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention meanwhile.
Claims (10)
1. a microearthquake event scanning localization method, is characterized in that, described microearthquake event scanning localization method comprises:
By many, ground, arrange and gather simple component or three-component microearthquake data;
Obtain check shot, and according to described check shot, generate the static correction value of each geophone station;
Utilize the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generate stack road;
Judge whether the number that peak amplitude in the arrangement of described predetermined number is greater than the arrangement of a default amplitude threshold reaches a setting number;
Timing definition while if so, the number of described arrangement being reached to described setting number is the time of arrival of a validity event;
Take-off time or the extreme value time of by picking up each, arranging, generate first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event;
Described validity event is carried out to static correction;
The indication using prestack seismic amplitude of the described validity event after static correction is multiplied by corresponding first arrival polarity p by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data;
After utilizing multiple dimensioned grid energy scan method to described correction, data are carried out focus location.
2. microearthquake event scanning localization method according to claim 1, is characterized in that, described in obtain check shot, and according to described check shot, generate the static correction value of each geophone station, comprising:
Utilize perforation or primacord fuse to obtain check shot;
Described check shot is picked up to described first break time, described first break time and initial velocity model are just being drilled to the first break time obtaining and subtracting each other, generate the static correction value of each geophone station.
3. microearthquake event scanning localization method according to claim 1, is characterized in that, described microearthquake event scanning localization method also comprises: the described amplitude threshold that described stack road is set.
4. microearthquake event scanning localization method according to claim 1, is characterized in that, described utilize multiple dimensioned grid energy scan method to described correction after data carry out focus location, comprising:
Subsurface source the first distribution range be take to the first yardstick subdivision as a plurality of the first grids;
While obtaining walking of the first grid described in each by ray-tracing scheme;
During according to the walking of described the first grid, data after described correction are carried out to TEC time error correction stack, generate the first superposition value;
Described first superposition value of each the first grid relatively, grid element center corresponding to described first superposition value of maximum of usining is as Primary Location point;
On described Primary Location point basis, divide the second distribution range, selecting the second yardstick is a plurality of the second grids by described the second distribution range subdivision;
While obtaining walking of the second grid described in each by ray-tracing scheme;
During according to the walking of described the second grid, data after described correction are carried out to TEC time error correction stack, generate the second superposition value;
Relatively described second superposition value of each the second grid, selects grid element center corresponding to maximum described the second superposition value as source point.
5. microearthquake event scanning localization method according to claim 4, is characterized in that, described the first distribution range is greater than described the second distribution range, and described the first yardstick is greater than described the second yardstick.
6. a microearthquake event scanning locating device, is characterized in that, described microearthquake event scanning locating device comprises:
Data acquisition unit, for being arranged and gathered simple component or three-component microearthquake data by many, ground;
Static correction value generation unit, for obtaining check shot, and generates the static correction value of each geophone station according to described check shot;
Stack road generation unit, for utilizing the first break time of described check shot described microearthquake data are carried out to TEC time error correction and superpose respectively by the arrangement of predetermined number, generates stack road;
Judging unit, for judging that the arrangement peak amplitude of described predetermined number is greater than a number of presetting the arrangement of amplitude threshold and whether reaches a setting number;
Time setting unit, the timing definition when the number of described arrangement is reached to described setting number is the time of arrival of a validity event;
First arrival polarity generation unit, for take-off time or the extreme value time of arranging by picking up each, generates first arrival polarity p
ij, wherein, p
ijrepresent j polarity of arranging of i event;
Static correction unit, for carrying out static correction to described validity event;
Data generating unit after proofreading and correct, for being multiplied by corresponding first arrival polarity p by the indication using prestack seismic amplitude of the described validity event after static correction by its corresponding event sequence number and arrangement sequence number
ij, generate and proofread and correct rear data;
Focus positioning unit, carries out focus location for data after utilizing multiple dimensioned grid energy scan method to described correction.
7. microearthquake event scanning locating device according to claim 6, is characterized in that, described static correction value generation unit comprises:
Check shot acquisition module, for utilizing perforation or primacord fuse to obtain check shot;
Static correction value generation module, for described check shot is picked up to described first break time, is just drilling by described first break time and initial velocity model the first break time obtaining and is subtracting each other, and generates the static correction value of each geophone station.
8. microearthquake event scanning locating device according to claim 6, is characterized in that, described microearthquake event scanning locating device also comprises: amplitude threshold setting unit, and for the described amplitude threshold in described stack road is set.
9. microearthquake event scanning locating device according to claim 6, is characterized in that, described focus positioning unit comprises:
The first mesh generation module, for take subsurface source the first distribution range the first yardstick subdivision as a plurality of the first grids;
The first acquisition module while walking, when obtaining the walking of the first grid described in each by ray-tracing scheme;
The first superposition value generation module, carries out TEC time error correction stack to data after described correction when according to the walking of described the first grid, and generates the first superposition value;
Primary Location point selection module, for comparing described first superposition value of each the first grid, selects grid element center corresponding to maximum described the first superposition value as Primary Location point;
The second mesh generation module, for divide the second distribution range on described Primary Location point basis, second yardstick of take is a plurality of the second grids by described the second distribution range subdivision;
The second acquisition module while walking, when obtaining the walking of the second grid described in each by ray-tracing scheme;
The second superposition value generation module, carries out TEC time error correction stack to data after described correction when according to the walking of described the second grid, and generates the second superposition value;
Source point is selected module, for comparing described second superposition value of each the second grid, selects grid element center corresponding to maximum described the second superposition value as source point.
10. microearthquake event scanning locating device according to claim 9, is characterized in that, described the first distribution range is greater than described the second distribution range, and described the first yardstick is greater than described the second yardstick.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410373504.8A CN104133246B (en) | 2014-07-31 | 2014-07-31 | A kind of micro-seismic event Scan orientation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410373504.8A CN104133246B (en) | 2014-07-31 | 2014-07-31 | A kind of micro-seismic event Scan orientation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104133246A true CN104133246A (en) | 2014-11-05 |
| CN104133246B CN104133246B (en) | 2016-11-02 |
Family
ID=51805989
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410373504.8A Active CN104133246B (en) | 2014-07-31 | 2014-07-31 | A kind of micro-seismic event Scan orientation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104133246B (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104360384A (en) * | 2014-11-13 | 2015-02-18 | 中国石油天然气集团公司 | Microseism event positioning method and device based on automatic scanning of longitudinal and transverse wave energy |
| CN105403918A (en) * | 2015-12-09 | 2016-03-16 | 中国科学院地质与地球物理研究所 | Three-component microseism data effective event identification method and system |
| CN105954795A (en) * | 2016-04-25 | 2016-09-21 | 吉林大学 | Grid successive dissection method used for microseismic positioning |
| CN106168679A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | The processing method of seismic acquisition record polarity |
| CN106199520A (en) * | 2015-05-30 | 2016-12-07 | 中国石油化工股份有限公司 | Microseism seismic source location method and apparatus |
| CN106249297A (en) * | 2015-06-08 | 2016-12-21 | 中国石油化工股份有限公司 | Fracturing microseism seismic source location method and system based on Signal estimation |
| CN107367754A (en) * | 2016-05-11 | 2017-11-21 | 中国石油化工股份有限公司 | Microseism first arrival recognition methods and device based on three-component polarization gradient |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
| CN106199520B (en) * | 2015-05-30 | 2018-08-31 | 中国石油化工股份有限公司 | Microseism seismic source location method and apparatus |
| CN109188516A (en) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | The micro-seismic event localization method of the domain Radon energy scan superposition |
| CN110907991A (en) * | 2019-12-11 | 2020-03-24 | 重庆大学 | Seismic source positioning method and system based on data field potential value and readable storage medium |
| CN110967751A (en) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium |
| CN110967762A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Microseism stratum velocity calibration method and system |
| CN112213768A (en) * | 2020-09-25 | 2021-01-12 | 南方科技大学 | Ground micro-seismic positioning method and system combining seismic source mechanism inversion |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
| US6032104A (en) * | 1997-02-20 | 2000-02-29 | Elf Exploration Production | Method of processing seismic signals so as to reduce the effects of dip in the analysis of moveouts |
| CN102841373A (en) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Microseism positioning method based on azimuth angle constraint |
| CN102928873A (en) * | 2012-10-30 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for positioning ground micro-seismic based on four-dimensional energy focusing |
| CN103064111A (en) * | 2012-12-12 | 2013-04-24 | 中国石油天然气集团公司 | Micro seismic event recognition method based on morphological filtering |
| CN103076634A (en) * | 2013-01-06 | 2013-05-01 | 刘建中 | Method for carrying out earthquake source mechanism analysis by utilizing micro earthquakes |
| CN103399300A (en) * | 2013-07-31 | 2013-11-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Wave packet superposition microseism ground location method |
| CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
-
2014
- 2014-07-31 CN CN201410373504.8A patent/CN104133246B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
| US6032104A (en) * | 1997-02-20 | 2000-02-29 | Elf Exploration Production | Method of processing seismic signals so as to reduce the effects of dip in the analysis of moveouts |
| CN102841373A (en) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Microseism positioning method based on azimuth angle constraint |
| CN102928873A (en) * | 2012-10-30 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for positioning ground micro-seismic based on four-dimensional energy focusing |
| CN103064111A (en) * | 2012-12-12 | 2013-04-24 | 中国石油天然气集团公司 | Micro seismic event recognition method based on morphological filtering |
| CN103076634A (en) * | 2013-01-06 | 2013-05-01 | 刘建中 | Method for carrying out earthquake source mechanism analysis by utilizing micro earthquakes |
| CN103399300A (en) * | 2013-07-31 | 2013-11-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Wave packet superposition microseism ground location method |
| CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
Non-Patent Citations (1)
| Title |
|---|
| 张栋，韩文功，贾静敏等: "横波微地震测井资料处理方法探索及初步应用", 《石油仪器》, vol. 27, no. 2, 30 April 2013 (2013-04-30), pages 62 - 65 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104360384A (en) * | 2014-11-13 | 2015-02-18 | 中国石油天然气集团公司 | Microseism event positioning method and device based on automatic scanning of longitudinal and transverse wave energy |
| CN104360384B (en) * | 2014-11-13 | 2017-03-08 | 中国石油天然气集团公司 | Micro-seismic event localization method based on wave energy automatically scanning in length and breadth and device |
| CN106168679B (en) * | 2015-05-18 | 2018-03-09 | 中国石油化工股份有限公司 | Seismic acquisition records the processing method of polarity |
| CN106168679A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | The processing method of seismic acquisition record polarity |
| CN106199520B (en) * | 2015-05-30 | 2018-08-31 | 中国石油化工股份有限公司 | Microseism seismic source location method and apparatus |
| CN106199520A (en) * | 2015-05-30 | 2016-12-07 | 中国石油化工股份有限公司 | Microseism seismic source location method and apparatus |
| CN106249297A (en) * | 2015-06-08 | 2016-12-21 | 中国石油化工股份有限公司 | Fracturing microseism seismic source location method and system based on Signal estimation |
| CN106249297B (en) * | 2015-06-08 | 2018-04-06 | 中国石油化工股份有限公司 | Hydraulic fracturing microseism seismic source location method and system based on signal estimation |
| CN105403918B (en) * | 2015-12-09 | 2018-06-22 | 中国科学院地质与地球物理研究所 | A kind of three-component microseism data validity event recognition methods and system |
| CN105403918A (en) * | 2015-12-09 | 2016-03-16 | 中国科学院地质与地球物理研究所 | Three-component microseism data effective event identification method and system |
| CN105954795A (en) * | 2016-04-25 | 2016-09-21 | 吉林大学 | Grid successive dissection method used for microseismic positioning |
| CN107367754A (en) * | 2016-05-11 | 2017-11-21 | 中国石油化工股份有限公司 | Microseism first arrival recognition methods and device based on three-component polarization gradient |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
| CN108445533B (en) * | 2018-03-08 | 2019-06-07 | 青岛海洋地质研究所 | The method of long offset building phytal zone rate pattern based on OBS |
| CN110967751A (en) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium |
| CN110967751B (en) * | 2018-09-29 | 2022-03-08 | 中国石油化工股份有限公司 | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium |
| CN110967762A (en) * | 2018-09-30 | 2020-04-07 | 中国石油化工股份有限公司 | Microseism stratum velocity calibration method and system |
| CN109188516B (en) * | 2018-10-31 | 2021-07-20 | 中国石油化工股份有限公司 | Microseism event positioning method for Radon domain energy scanning and stacking |
| CN109188516A (en) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | The micro-seismic event localization method of the domain Radon energy scan superposition |
| CN110907991A (en) * | 2019-12-11 | 2020-03-24 | 重庆大学 | Seismic source positioning method and system based on data field potential value and readable storage medium |
| CN110907991B (en) * | 2019-12-11 | 2021-03-16 | 重庆大学 | Seismic source positioning method and system based on data field potential value and readable storage medium |
| CN112213768A (en) * | 2020-09-25 | 2021-01-12 | 南方科技大学 | Ground micro-seismic positioning method and system combining seismic source mechanism inversion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104133246B (en) | 2016-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104133246A (en) | Microseism event scanning positioning method and device | |
| CN105807316B (en) | Ground observation microseism velocity model corrections method based on amplitude superposition | |
| CN103105624A (en) | Longitudinal and transversal wave time difference positioning method based on base data technology | |
| Zhou et al. | Machine‐learning‐facilitated earthquake and anthropogenic source detections near the Weiyuan Shale Gas Blocks, Sichuan, China | |
| CN104280775B (en) | Microseism monitoring and positioning method based on full-waveform vector offset superposition | |
| CN103336297B (en) | Micro rupture vector scan method | |
| CN103984024B (en) | A kind of three-component geophone horizontal component data auto-correction method | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN107132578B (en) | A kind of microseism ground monitoring velocity model corrections algorithm | |
| CN106680869A (en) | Microearthquake event detecting and positioning method and apparatus | |
| CN103869363B (en) | Microseism localization method and device | |
| CN105301639A (en) | Speed field inversion method and device based on VSP double-weight travel time tomography | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN102928873A (en) | Method for positioning ground micro-seismic based on four-dimensional energy focusing | |
| CN101915938A (en) | Offset imaging method and device for converted waves | |
| CN110516650A (en) | A kind of blind space seismic source location system of shallow-layer based on shock sensor array | |
| CN104765064A (en) | Microseism interference imaging method | |
| CN107703540B (en) | A kind of microseism positioning and chromatography imaging method | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| CN103105622A (en) | Homomorphous wave time difference positioning method based on data base technology | |
| CN105549077A (en) | Micro-earthquake epicenter positioning method calculated based on multilevel multi-scale grid similarity coefficient | |
| Rowan et al. | The three‐dimensional structure of Kilauea Volcano, Hawaii, from travel time tomography | |
| CN104614762B (en) | Loose sandstone gas reservoir boundary determining method and device | |
| CN103630938A (en) | Imaging system and imaging method for well earthquake using hammer head of down-hole hammer as focus | |
| CN105717538B (en) | Undulating surface seismic data migration datum plane conversion method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |