CN104597504A - High-efficiency microseism seismic focus locating method - Google Patents
High-efficiency microseism seismic focus locating method Download PDFInfo
- Publication number
- CN104597504A CN104597504A CN201510047191.1A CN201510047191A CN104597504A CN 104597504 A CN104597504 A CN 104597504A CN 201510047191 A CN201510047191 A CN 201510047191A CN 104597504 A CN104597504 A CN 104597504A
- Authority
- CN
- China
- Prior art keywords
- wave
- source location
- microearthquake
- seismic
- filial generation
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a high-efficiency microseism seismic focus locating method which specifically includes the following steps: placing a plurality of detectors for P wave and S wave in a monitoring well nearby a fractured well in advance, collecting the P wave and the S wave generated by seism of the fractured well through the detectors after the P wave and the S wave reach the monitoring well through a stratum; sampling the detectors, converting simulation of the P wave and the S wave into digital signals capable of being recognized and displayed by a computer, recording the signals in binary files; utilizing seismic wave processing software to read the binary files to acquire the P wave data and the S wave data required by calculation of the location of the seismic focus, including parameters such as time and amplitude, and then conducting quick efficient calculation according to a microseism seismic focus locating method to determine the location of fracture in the fractured well. The method has the advantages of being high in efficiency, high in accuracy, small in errors and the like.
Description
Technical field
The present invention relates to a kind of high efficiency microearthquake seismic source location method, be applied to microearthquake seismic source location field.
Background technology
Oil/gas Well is because of the needs of increase yield, carry out waterfrac treatment, when hydraulic pressure is greater than reservoir stress, stratum can produce breaks, and then bring out microearthquake, produce seismic event, mainly P involves S ripple, existing microearthquake seismic source location method scheme be little pane on four dimensional spacetime territory search plain technology, this technology needs to travel through in space of thinking all searches plain region, and it is low to search plain efficiency, this technology also depends on the width of little pane simultaneously, search precision is not high, and this technology depends on the configuration of machine hardware and the width of little pane in addition, and therefore positioning error is bigger than normal.
Summary of the invention
For solving the defect of prior art, the invention discloses a kind of high efficiency microearthquake seismic source location method that efficiency is high, precision is high, error is little.
The invention discloses a kind of high efficiency microearthquake seismic source location method, concrete steps are as follows:
Step one: in the monitor well near fractured well, places multiple wave detector receiving P, S ripple in advance, P, S ripple that fractured well microearthquake produces, and is, after medium reaches monitor well, be detected device and collect with stratum;
Step 2: sample to wave detector, by P, S wave simulation signal, is converted into the digital signal of computing machine identifiable design, display, is recorded in binary file;
Step 3: use Seismic wave processing software initiatively to read these binary files, necessary P, S wave datum in source location is calculated to obtain, comprise the parameter such as the time of reaching, amplitude, then according to microearthquake seismic source location algorithm, fast, efficient calculation determine to break and occur in the position of fractured well.
Microearthquake seismic source location algorithm in a kind of high efficiency microearthquake seismic source location method of the present invention described in step 3 is specific as follows:
For objective function
Feasible zone X is nonempty set, and target is for finding
Wherein
;
Be described as follows:
1. initialization
Input evolution parameter: population scale N, crossover probability Pc, intersection factor F ∈ (0,1), evolutionary generation G, the bound of independent variable, in seismic source location, population scale is exactly separate the number of separating in set in initialized source location, general population scale is set to 20, crossover probability is set to 0.9, and the factor of intersecting is set to 0.5, and evolutionary generation is set to 1000 ~ 2000;
Random initial population:
Wherein each X(0) represent an initialized source location and separate (x, y, z, t) and comprise three dimensional space coordinate and origin time of earthquake four dimensions;
2. evaluate
Calculate
according to the parameter that source location is separated, push away to obtain the P ripple of waveform on each wave detector and the first arrival-time of S ripple according to this parameter is counter, according to amplitude of wave form value then, obtain an evaluation of estimate of separating about this position, evaluation of estimate is higher, then illustrate that position solution is the closer to real source location;
3. breeding restructuring
To each individuality in population
, stochastic generation three mutually different random integers
random integers
Wherein n represents the dimension that source location is separated, in the method, and n=4.
Each dimension in (x, y, z, t) is separated to source location, produces filial generation according to above-mentioned strategy, produce X(1 by X (0)), by that analogy follow-up;
4. select
To each filial generation, according to evaluation function
, the filial generation evaluation function as produced is more excellent, then replace the parent of its correspondence with filial generation, if filial generation evaluation function is more unexcellent, then directly inherit filial generation with parent;
If
meet end condition (reaching maximum evolutionary generation), with regard to exporting, there is optimum individual, if do not meet end condition, just return second step, continue evaluate and recombinate.
One of the present invention high efficiency microearthquake seismic source location method, adopt Basic Differential evolvement method, Basic Differential evolvement method is applicable to the Global Optimal Problem without constraint continuous variable, comprise linear programming, nonlinear programming, non-smooth blade, can also solve discrete variable problem for dividing evolvement method and retraining optimization problem after expansion, eliminate the region-wide of four-dimentional space and search element, location efficiency is high, without the discrete space of little pane, but continuous space search element, positioning precision is high, the discrete space of little pane and the error of time are abolished, positioning error is little.
Accompanying drawing explanation
Fig. 1 is the FB(flow block) of a kind of high efficiency microearthquake seismic source location method of the present invention.
Embodiment
As shown in Figure 1, the invention discloses a kind of high efficiency microearthquake seismic source location method, concrete steps are as follows:
Step one: in the monitor well near fractured well, places multiple wave detector receiving P, S ripple in advance, P, S ripple that fractured well microearthquake produces, and is, after medium reaches monitor well, be detected device and collect with stratum;
Step 2: sample to wave detector, by P, S wave simulation signal, is converted into the digital signal of computing machine identifiable design, display, is recorded in binary file;
Step 3: use Seismic wave processing software initiatively to read these binary files, necessary P, S wave datum in source location is calculated to obtain, comprise the parameter such as the time of reaching, amplitude, then according to microearthquake seismic source location algorithm, fast, efficient calculation determine to break and occur in the position of fractured well.
Microearthquake seismic source location algorithm in a kind of high efficiency microearthquake seismic source location method of the present invention described in step 3 is specific as follows:
For objective function
Feasible zone X is nonempty set, and target is for finding
Wherein
;
Be described as follows:
1. initialization
Input evolution parameter: population scale N, crossover probability Pc, intersection factor F ∈ (0,1), evolutionary generation G, the bound of independent variable, in seismic source location, population scale is exactly separate the number of separating in set in initialized source location, general population scale is set to 20, crossover probability is set to 0.9, and the factor of intersecting is set to 0.5, and evolutionary generation is set to 1000 ~ 2000;
Random initial population:
Wherein each X(0) represent an initialized source location and separate (x, y, z, t) and comprise three dimensional space coordinate and origin time of earthquake four dimensions;
2. evaluate
Calculate
according to the parameter that source location is separated, push away to obtain the P ripple of waveform on each wave detector and the first arrival-time of S ripple according to this parameter is counter, according to amplitude of wave form value then, obtain an evaluation of estimate of separating about this position, evaluation of estimate is higher, then illustrate that position solution is the closer to real source location;
3. breeding restructuring
To each individuality in population
, stochastic generation three mutually different random integers
random integers
Wherein n represents the dimension that source location is separated, in the method, and n=4.
Each dimension in (x, y, z, t) is separated to source location, produces filial generation according to above-mentioned strategy, produce X(1 by X (0)), by that analogy follow-up;
4. select
To each filial generation, according to evaluation function
, the filial generation evaluation function as produced is more excellent, then replace the parent of its correspondence with filial generation, if filial generation evaluation function is more unexcellent, then directly inherit filial generation with parent;
If
meet end condition (reaching maximum evolutionary generation), with regard to exporting, there is optimum individual, if do not meet end condition, just return second step, continue evaluate and recombinate.
Last it is noted that above embodiment only in order to illustrate the present invention and and unrestricted technical scheme described in the invention; Therefore, although this instructions with reference to each above-mentioned embodiment to present invention has been detailed description, those of ordinary skill in the art should be appreciated that and still can modify to the present invention or equivalent to replace; And all do not depart from technical scheme and the improvement thereof of the spirit and scope of the present invention, it all should be encompassed in right of the present invention.
Claims (2)
1. a high efficiency microearthquake seismic source location method, is characterized in that: concrete steps are as follows:
Step one: in the monitor well near fractured well, places multiple wave detector receiving P, S ripple in advance, P, S ripple that fractured well microearthquake produces, and is, after medium reaches monitor well, be detected device and collect with stratum;
Step 2: sample to wave detector, by P, S wave simulation signal, is converted into the digital signal of computing machine identifiable design, display, is recorded in binary file;
Step 3: use Seismic wave processing software initiatively to read these binary files, necessary P, S wave datum in source location is calculated to obtain, comprise the parameter such as the time of reaching, amplitude, then according to microearthquake seismic source location algorithm, fast, efficient calculation determine to break and occur in the position of fractured well.
2. one according to claim 1 high efficiency microearthquake seismic source location method, is characterized in that: the microearthquake seismic source location algorithm described in step 3 is specific as follows:
For objective function
Feasible zone X is nonempty set, and target is for finding
Wherein
;
Be described as follows:
(1) initialization
Input evolution parameter: population scale N, crossover probability Pc, intersection factor F ∈ (0,1), evolutionary generation G, the bound of independent variable, in seismic source location, population scale is exactly separate the number of separating in set in initialized source location, general population scale is set to 20, crossover probability is set to 0.9, and the factor of intersecting is set to 0.5, and evolutionary generation is set to 1000 ~ 2000;
Random initial population:
Wherein each X(0) represent an initialized source location and separate (x, y, z, t) and comprise three dimensional space coordinate and origin time of earthquake four dimensions;
(2) evaluate
Calculate
according to the parameter that source location is separated, push away to obtain the P ripple of waveform on each wave detector and the first arrival-time of S ripple according to this parameter is counter, according to amplitude of wave form value then, obtain an evaluation of estimate of separating about this position, evaluation of estimate is higher, then illustrate that position solution is the closer to real source location;
(3) breeding restructuring
To each individuality in population
, stochastic generation three mutually different random integers
random integers
Wherein n represents the dimension that source location is separated, in the method, and n=4.
Each dimension in (x, y, z, t) is separated to source location, produces filial generation according to above-mentioned strategy, produce X(1 by X (0)), by that analogy follow-up;
(4) select
To each filial generation, according to evaluation function
, the filial generation evaluation function as produced is more excellent, then replace the parent of its correspondence with filial generation, if filial generation evaluation function is more unexcellent, then directly inherit filial generation with parent;
If
meet end condition (reaching maximum evolutionary generation), with regard to exporting, there is optimum individual, if do not meet end condition, just return second step, continue evaluate and recombinate.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510047191.1A CN104597504A (en) | 2015-01-30 | 2015-01-30 | High-efficiency microseism seismic focus locating method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510047191.1A CN104597504A (en) | 2015-01-30 | 2015-01-30 | High-efficiency microseism seismic focus locating method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104597504A true CN104597504A (en) | 2015-05-06 |
Family
ID=53123409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510047191.1A Pending CN104597504A (en) | 2015-01-30 | 2015-01-30 | High-efficiency microseism seismic focus locating method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104597504A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104832160A (en) * | 2015-05-21 | 2015-08-12 | 河北煤炭科学研究院 | Downhole drilling trajectory microseism description method |
CN105842735A (en) * | 2016-05-20 | 2016-08-10 | 四川大学 | Complex-velocity-distribution regional rock micro-seismic seismic source positioning method |
CN106249295A (en) * | 2015-06-15 | 2016-12-21 | 中国石油化工股份有限公司 | A kind of borehole microseismic P, S ripple associating method for rapidly positioning and system |
CN110967751A (en) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101968549A (en) * | 2010-09-14 | 2011-02-09 | 北京安科兴业科技有限公司 | Small mine earthquake monitoring system and radio-detector arrangement method |
CN102495425A (en) * | 2011-11-14 | 2012-06-13 | 北京科技大学 | Energy-based method for automatically locating earthquake focus of microearthquake |
CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
-
2015
- 2015-01-30 CN CN201510047191.1A patent/CN104597504A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101968549A (en) * | 2010-09-14 | 2011-02-09 | 北京安科兴业科技有限公司 | Small mine earthquake monitoring system and radio-detector arrangement method |
CN102495425A (en) * | 2011-11-14 | 2012-06-13 | 北京科技大学 | Energy-based method for automatically locating earthquake focus of microearthquake |
CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
Non-Patent Citations (1)
Title |
---|
刘云彬等: "用微地震波法确定松辽盆地头台油田地应力场特征", 《大庆石油地质与开发》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104832160A (en) * | 2015-05-21 | 2015-08-12 | 河北煤炭科学研究院 | Downhole drilling trajectory microseism description method |
CN106249295A (en) * | 2015-06-15 | 2016-12-21 | 中国石油化工股份有限公司 | A kind of borehole microseismic P, S ripple associating method for rapidly positioning and system |
CN105842735A (en) * | 2016-05-20 | 2016-08-10 | 四川大学 | Complex-velocity-distribution regional rock micro-seismic seismic source positioning method |
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 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9880012B2 (en) | Hybrid road network and grid based spatial-temporal indexing under missing road links | |
Goswami | A brief introduction to nonlinear time series analysis and recurrence plots | |
US10229129B2 (en) | Method and apparatus for managing time series database | |
RU2594373C1 (en) | Identification of prevailing orientations of cracks | |
CN100456281C (en) | Data division apparatus, data division method | |
CN104597504A (en) | High-efficiency microseism seismic focus locating method | |
CN105224543A (en) | For the treatment of seasonal effect in time series method and apparatus | |
CN103714004A (en) | JVM online memory leak analysis method and system | |
CN105224660A (en) | A kind of disposal route of map point of interest POI data and device | |
CN105379186A (en) | Determining response similarity neighborhoods | |
KR20200052424A (en) | Method for segmenting time series data and apparatus thereof | |
CN103853752A (en) | Method and device for managing time series database | |
Luna-Romera et al. | An approach to silhouette and dunn clustering indices applied to big data in spark | |
CN103544325A (en) | Data processing method and device used for web page clicking distribution | |
Kim et al. | Applications of clustering and isolation forest techniques in real-time building energy-consumption data: Application to LEED certified buildings | |
CA3022158C (en) | Method and system for positioning seismic source in microseism monitoring | |
US20160163102A1 (en) | Computing device and method for processing point clouds | |
CN104598473A (en) | Information processing method and electronic device | |
Pouchard et al. | Prescriptive provenance for streaming analysis of workflows at scale | |
CN103678111A (en) | Evaluation method and evaluation device for similarity of source codes | |
Hadimlioglu et al. | Visualization of flooding using adaptive spatial resolution | |
US20160178590A1 (en) | System and method for predicting harmful materials | |
CN103914027A (en) | Method for monitoring and managing industrial on-site bottom layer device through expression analysis | |
Dréau et al. | Multi-element polynomial chaos expansion based on automatic discontinuity detection for nonlinear systems | |
CN104155692A (en) | Method and device for obtaining focus point of microseism |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150506 |