CN105403917A - Crack detection method and device - Google Patents
Crack detection method and device Download PDFInfo
- Publication number
- CN105403917A CN105403917A CN201510726803.XA CN201510726803A CN105403917A CN 105403917 A CN105403917 A CN 105403917A CN 201510726803 A CN201510726803 A CN 201510726803A CN 105403917 A CN105403917 A CN 105403917A
- Authority
- CN
- China
- Prior art keywords
- earthquake data
- dimensional earthquake
- sector
- dimensional
- stress field
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 7
- 239000000203 mixture Substances 0.000 claims description 14
- 241001269238 Data Species 0.000 claims description 8
- 238000013481 data capture Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 13
- 238000003384 imaging method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 239000002245 particle Substances 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 239000010430 carbonatite Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical group 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a crack detection method and device. The method comprises obtaining multi-dimensional seismic data of a target work area, determining a principal stress field direction of the target work area according to analysis processing of the multi-dimensional seismic data, arranging the multi-dimensional seismic data to be sub-azimuth multi-dimensional seismic data according to orientation angles and the principle stress field direction, and performing ellipse fitting according to seismic parameter distribution of a sub-azimuth multi-dimensional data set to obtain a growth direction and growth strength of a crack. Through utilization of the technical scheme provided in an embodiment of the invention, accuracy of determination of the growth direction and strength of the crack can be improved.
Description
Technical field
The present invention relates to seismogeology exploratory development interpretation technique field, particularly relate to a kind of method and device of Crack Detection.
Background technology
In the exploration and exploitation of oil and natural gas, the seam hole type hydrocarbon-bearing pool that carbonate fracture and hole are combined to form becomes the major domain of Tarim Oilfield carbonate exploration exploitation.Along with the lasting propelling of seam hole, oil field from the exploration stage to the development phase, the prediction in crack has become the important step of carbonatite Efficient Development.The Accurate Prediction in crack can have larger control action in exploitation links such as evading drilling engineering risk, raising reservoir Drilling ratio, Carbonate acidizing pressure break reservoir reconstruction.
Crack detection method of the prior art mainly comprises shear wave splitting, microearthquake and compressional wave FRACTURE PREDICTION technology.Because the acquisition cost of shear wave is higher, the problems such as the monitoring range of microearthquake is limited, Crack Detection technology comparatively conventional is at present compressional wave FRACTURE PREDICTION technology.Compressional wave FRACTURE PREDICTION technology usually prestack can be divided into predict according to the type of seismic data and poststack predicts two classes.Poststack compressional wave FRACTURE PREDICTION technology mainly through horizontal uncontinuities in geometric attribute Study of Seismic road such as relevant, curvature, according to the large fracture feature that described uncontinuity research causes due to tomography, little fracture etc.Prestack compressional wave FRACTURE PREDICTION technology can carry out FRACTURE PREDICTION by Prestack seismic data orientative feature, the main process of prestack compressional wave azimuthal anisotropy crack prediction method is as follows: calculate different azimuth attribute, ellipse fitting is carried out to different azimuth attribute, transverse can reflection crack orientation, and oval short major axis ratio can reflection crack density.
Crack prediction method of the prior art does not consider that Crack Detection result is subject to the impact of geology and seismic condition.Because carbonatite nonuniformity is stronger, differing greatly between adjacent seismic trace, adopting method of the prior art to carry out, that FRACTURE PREDICTION may cause predicting the outcome is inaccurate.
Summary of the invention
The object of the present invention is to provide a kind of method and apparatus of Crack Detection, the accuracy of FRACTURE PREDICTION can be improved.
Method and the device of a kind of Crack Detection provided by the invention are achieved in that
A method for Crack Detection, described method comprises:
Obtain the multi-dimensional earthquake data in target work area;
According to the analyzing and processing to described multi-dimensional earthquake data, determine the principal stress field direction in described target work area;
According to described principal stress field direction, described multi-dimensional earthquake data are arranged composition orientation multi-dimensional earthquake data according to position angle;
Ellipse fitting is carried out in seismologic parameter distribution according to described point of orientation multi-dimensional earthquake data set, obtains growth direction and the growth strength in crack.
Optionally, in one embodiment of the invention, described multi-dimensional earthquake data comprise based on longitudinal survey line, horizontal survey line, degree of depth survey line, offset distance, azimuthal five dimension geological datas.
Optionally, in one embodiment of the invention, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
The distribution of the sector divided in adjustment point orientation multi-dimensional earthquake data, until described principal stress field direction is arranged in complete sector.
Optionally, in one embodiment of the invention, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
According to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the degree of covering of the sector of different directions matches in described point of orientation multi-dimensional earthquake data.
Optionally, in one embodiment of the invention, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
The number of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the signal to noise ratio (S/N ratio) of geological data and the degree of covering of described sector match in described sector.
Optionally, in one embodiment of the invention, described seismologic parameter comprise following in any one: amplitude, frequency, when walking.
A device for Crack Detection, described device comprises:
Multi-dimensional earthquake data capture unit, for obtaining the multi-dimensional earthquake data in target work area;
Principal stress field direction-determining unit, for according to the analyzing and processing to described multi-dimensional earthquake data, determines the principal stress field direction in described target work area;
Divide orientation multi-dimensional earthquake data ordering unit, for described multi-dimensional earthquake data being arranged composition orientation multi-dimensional earthquake data according to position angle according to described principal stress field direction;
Ellipse fitting unit, carries out ellipse fitting for the seismologic parameter distribution according to described point of orientation multi-dimensional earthquake data set, obtains growth direction and the growth strength in crack.
Optionally, in one embodiment of the invention, described multi-dimensional earthquake data comprise based on longitudinal survey line, horizontal survey line, degree of depth survey line, offset distance, azimuthal five dimension geological datas.
Optionally, in one embodiment of the invention, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector distribution adjustment unit, for adjusting the distribution of the sector divided in point orientation multi-dimensional earthquake data, until described principal stress field direction is arranged in complete sector.
Optionally, in one embodiment of the invention, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector-size adjustment unit, for according to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the degree of covering of the sector of different directions matches in described point of orientation multi-dimensional earthquake data.
Optionally, in one embodiment of the invention, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector number adjustment unit, for adjusting the number of the sector divided in point orientation multi-dimensional earthquake data, until the signal to noise ratio (S/N ratio) of geological data and the degree of covering of described sector match in described sector.
Optionally, in one embodiment of the invention, described seismologic parameter comprise following in any one: amplitude, frequency, when walking.
The technical scheme of the method and apparatus of a kind of Crack Detection of the present invention considers the anisotropy of geologic condition, according to principal stress field direction and point orientation multi-dimensional earthquake data in the multi-dimensional earthquake data acquisition target work area in target work area, macroscopic fracture direction can be obtained from principal stress field direction, improve the accuracy determining micro crack.Ellipse fitting is carried out in seismologic parameter distribution according to point orientation multi-dimensional earthquake data, can obtain seismologic parameter value, improve the accuracy of ellipse fitting, improve and determine the growth direction in crack and the accuracy of growth strength from multiple orientation.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the method flow diagram of a kind of embodiment of crack prediction method provided by the invention;
Fig. 2 is oil well distribution and the principal stress field directional diagram of study area, oil field in Tarim Basin;
Fig. 3 is the front and back comparison diagram that orientation multi-dimensional earthquake data sector distribution is distinguished in adjustment provided by the invention research;
Fig. 4 is comparison diagram study area good Li Tage group data being carried out to ellipse fitting by different azimuth number;
Fig. 5 is comparison diagram hawk mountain, study area group data being carried out to ellipse fitting by different azimuth number;
Fig. 6 is the modular structure schematic diagram of a kind of embodiment of crack detecting device provided by the invention.
Embodiment
Technical scheme in the present invention is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present invention, 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, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, should belong to the scope of protection of the invention.
The present invention is based on the anisotropy of geologic condition in formation at target locations, take stress field as foundation, carries out Crack Detection from preferred seismologic parameter data to carbonatite formation at target locations.The invention provides a kind of method and device of Crack Detection, Fig. 1 is the method flow diagram of a kind of embodiment of crack detection method provided by the invention, and as shown in Figure 1, described method can comprise:
Step S1: the multi-dimensional earthquake data obtaining target work area.
In seismic data process, often two-dimension earthquake data or 3-d seismic data set are studied, but in some application scenarios, two dimension or 3D seismic data can not meet the demand of research.The geological data only having single orientation to generate may comprise the geophysical structure of accuracy deficiency.Such as, for the stratum with effect of anisotropy, can show the geology imaging obtained according to different position angles may be significantly different.With anticipate accuracy to geophysical structure such as tomography, little vertical displacement and the fracture of secondary earthquake scale (such as, the crack that measurement result is less than tens of meters, described crack may lower than the resolution of the detection of receiver or other detection instrument) carry out imaging and may need to carry out imaging along each position angle.The geological data obtained from multiple position angle can be called wide-azimuth geological data, and wide-azimuth geological data can have higher value for the imaging of carbonatite structure (such as Tarim Oilfield).Based on above-mentioned theory analysis, the present invention can carry out Crack Detection based on multi-dimensional earthquake data to described target work area.In one embodiment of the invention, described multi-dimensional earthquake data comprise based on longitudinal survey line, horizontal survey line, degree of depth survey line, offset distance, azimuthal five dimension geological datas, and described five dimension geological datas can be the 3D seismic datas based on multiple offset distance wide-azimuth.Utilize described multi-dimensional earthquake data instead of normally used two dimension, 3D seismic data carries out imaging to geophysical structure and can improve the accuracy of image and provide the extraneous information about geologic structure.Such as, wide-azimuth geological data can comprise different azimuthal multiple image traces, so from target work area described in multiple angle analysis the structure in crack.Similarly, multiple offset distance geological data progressively can change the offset distance between emitting antenna and receiving antenna and the multi-fold data gathered from formation at target locations in same position.Multiple offset distance geological data is comparatively large relative to the quantity of information of single offset distance geological data, by the signal to noise ratio (S/N ratio) that can improve geological data after described multiple offset distance geological data superposition.
Step S2: according to the analyzing and processing to described multi-dimensional earthquake data, determines the principal stress field direction in described target work area.
In stratal configuration, the spatial composing relations of reservoir and tight zone can determine longitudinal development models in crack.Usually, the genesis analysis in crack is mainly subject to the control of stress field and stratal configuration.When terrestrial stress increases slowly and rock stratum pressure is major principal stress, comparatively loose and the rock stratum intergranular pore that oozing property of hole is higher is grown better, grain compressive strength is little, and mainly with the form of Particle Breakage release terrestrial stress, crack shows as Particle Breakage crack more.Hole between reservoir ooze densification, rock compressive strength large, particle not easily breaks, and along with the release of adjacent reservoirs section terrestrial stress, the terrestrial stress in fine and close reservoir is more concentrated, finally with the form of rock breakdown release stress, can produce macroscopic fracture.Described Particle Breakage crack is often distributed in the adjacent domain of macroscopic fracture.Therefore, under the condition that mechanical properties of rock is close, the growth direction in crack and development degree and principal stress field direction are closely connected.
In one embodiment of the invention, on the basis understanding described target work area structural setting, according to the analysis to described multi-dimensional earthquake data, the principal stress field direction in described target work area can be determined.By determining principal stress field direction, at least can determine the direction of some macroscopic fractures, after determining the direction of macroscopic fracture, for detection micro crack (such as Particle Breakage splits), can detection time be saved, increase work efficiency.
Step S3: described multi-dimensional earthquake data are arranged composition orientation multi-dimensional earthquake data according to position angle according to described principal stress field direction.
Described multi-dimensional earthquake data in the present invention can be the geological datas based on wide-azimuth, in one embodiment of the invention, ellipse fitting method can be adopted to carry out matching to described multi-dimensional earthquake data, and then obtain growth direction and the growth strength in crack.So, first by according to described principal stress field direction, described multi-dimensional earthquake data can be arranged composition orientation multi-dimensional earthquake data according to position angle.Described point of orientation multi-dimensional earthquake data can be that described multi-dimensional earthquake data are divided into several sector regions according to 360 degree of position angles.
Particularly, in one embodiment of the invention, the distribution of the sector divided in point orientation multi-dimensional earthquake data can be adjusted, make described principal stress field direction be arranged in complete sector.According to above-mentioned analysis, the direction of principal stress field can represent the direction of macroscopic fracture, when dividing the sector of point orientation multi-dimensional earthquake data, situation about being just in time divided in principal stress field direction on the separatrix of two adjacent sectors easily occurs.That is, described macroscopic fracture is divided in two sectors, has a negative impact to follow-up process.At this time, can adjust the distribution of sector, thus make described principal stress field direction be arranged in complete sector.Such as, Fig. 3 is the front and back comparison diagram of adjustment provided by the invention point orientation multi-dimensional earthquake data sector distribution.As shown in Figure 3, the left figure of Fig. 3 is point orientation multi-dimensional earthquake data before the distribution of adjustment sector, as described in Figure, on the separatrix of just in time displacement two adjacent sectors, two principal stress field directions, after adjustment, shown in figure as right in Fig. 3, two principal stress field directions are arranged in complete sector.
In another embodiment of the present invention, can according to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, makes the degree of covering of the sector of different directions in described point orientation multi-dimensional earthquake data match.As shown in Figure 3, in some sectors of point orientation multi-dimensional earthquake data, indicate the degree of covering of geological data, usually, degree of covering is more, and the energy of sector is larger.In follow-up ellipse fitting process, can analyze under the condition of each sector balancing energy, determine the direction and intensity in crack according to the difference of earthquake parameter value.Particularly, by the size of adjustment sector, the degree of covering of the sector of different directions in described point orientation multi-dimensional earthquake data can be matched, ensures the equilibrium of described sector energy.
In another embodiment of the present invention, the number of the sector divided in point orientation multi-dimensional earthquake data can be adjusted, the degree of covering of the signal to noise ratio (S/N ratio) of geological data in described sector and described sector is matched.The resolution of signal to noise ratio (S/N ratio) and imaging has substantial connection, and the sector that signal to noise ratio (S/N ratio) is high generally needs less degree of covering just can imaging, and the lower sector of signal to noise ratio (S/N ratio) generally needs more degree of covering just can imaging.The number of sector can be adjusted, the degree of covering of the signal to noise ratio (S/N ratio) of geological data in described sector and described sector is matched.Particularly, when the signal to noise ratio (S/N ratio) of geological data is higher, can sector be increased, thus reduce the energy of single sector; When the signal to noise ratio (S/N ratio) of geological data is on the low side, can sector be reduced, thus increase the energy of single sector.
Step S4: the parameter distribution according to described point of orientation multi-dimensional earthquake data set carries out ellipse fitting, obtains growth direction and the growth strength in crack.
By described multi-dimensional earthquake data according to after arrangement composition orientation, position angle multi-dimensional earthquake data, the parameter value of seismologic parameter can be extracted from sector, carry out ellipse fitting according to described parameter value.Described seismologic parameter can comprise following in any one: amplitude, frequency, when walking.Particularly, can extract from each sector of described point of orientation multi-dimensional earthquake data above-mentioned in the parameter value of any one parameter and position angle corresponding to described parameter value.Such as, point orientation multi-dimensional earthquake data being divided into 4 sectors from one extract amplitude, an amplitude and position angle corresponding to described amplitude can be chosen from described 4 sectors respectively, according to 4 amplitudes and position angle corresponding to described amplitude, fit to an elliptic equation.The major axis of described elliptic equation can represent the growth direction in crack, and the minor axis of described elliptic equation can represent the growth strength in crack.
Below by a concrete application scenarios explanation, Fig. 2 is oil well distribution and the principal stress field directional diagram of study area, oil field in Tarim Basin, as shown in Figure 2, two groups of major developmental fractures of described study area can comprise direction, 60 degree, east by north and direction, 120 degree, east by north, can determine that the principal stress field direction of described study area is direction, 60 degree, east by north and direction, 120 degree, east by north.Fig. 3 is the front and back comparison diagram to described study area adjustment point orientation multi-dimensional earthquake data sector distribution, does not repeat them here.Fig. 4 is the comparison diagram described study area good Li Tage group data being carried out to ellipse fitting by different azimuth number, as shown in Figure 4, the crack pattern obtained after carrying out ellipse fitting after described good Li Tage group data are divided into 18 directions is more close to principal stress field, can determine for good Li Tage group data, the direction divided is more, and result is more accurate.Fig. 5 is comparison diagram hawk mountain, described study area group data being carried out to ellipse fitting by different azimuth number, because hawk mountain group intrinsic signal is more weak, orientation can keep the signal to noise ratio (S/N ratio) of data less, as shown in Figure 5, the crack pattern obtained after carrying out ellipse fitting after therefore described hawk mountain group data being divided into 6 directions is more close to principal stress field.
The present invention also provides a kind of device of Crack Detection on the other hand, Fig. 6 is the modular structure schematic diagram of a kind of embodiment of crack detecting device provided by the invention, by reference to the accompanying drawings 6, this device 60 can comprise: multi-dimensional earthquake data capture unit 61, principal stress field direction-determining unit 62, point orientation multi-dimensional earthquake data ordering unit 63, ellipse fitting unit 64, wherein
Multi-dimensional earthquake data capture unit 61, for obtaining the multi-dimensional earthquake data in target work area.
Principal stress field direction-determining unit 62, for according to the analyzing and processing to described multi-dimensional earthquake data, determines the principal stress field direction in described target work area.
Divide orientation multi-dimensional earthquake data ordering unit 63, for described multi-dimensional earthquake data being arranged composition orientation multi-dimensional earthquake data according to position angle according to described principal stress field direction.
Ellipse fitting unit 64, carries out ellipse fitting for the seismologic parameter distribution according to described point of orientation multi-dimensional earthquake data set, obtains growth direction and the growth strength in crack.
Described point of orientation multi-dimensional earthquake data ordering unit can comprise: sector distribution adjustment unit, for adjusting the distribution of the sector divided in point orientation multi-dimensional earthquake data, makes described principal stress field direction be arranged in complete sector.
Described point of orientation multi-dimensional earthquake data ordering unit can comprise: sector-size adjustment unit, for according to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, makes the degree of covering of the sector of different directions in described point orientation multi-dimensional earthquake data match.
Described point of orientation multi-dimensional earthquake data ordering unit can comprise: sector number adjustment unit, for adjusting the number of the sector divided in point orientation multi-dimensional earthquake data, the degree of covering of the signal to noise ratio (S/N ratio) of geological data in described sector and described sector is matched.
As can be seen here, the technical scheme of the method and apparatus of a kind of Crack Detection of the present invention considers the anisotropy of geologic condition, according to principal stress field direction and point orientation multi-dimensional earthquake data in the multi-dimensional earthquake data acquisition target work area in target work area, macroscopic fracture direction can be obtained from principal stress field direction, improve the accuracy determining micro crack.Ellipse fitting is carried out in seismologic parameter distribution according to point orientation multi-dimensional earthquake data, can obtain seismologic parameter value, improve the accuracy of ellipse fitting, and then the growth direction in crack and the accuracy of growth strength is determined in raising from multiple orientation.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiments, between each embodiment identical similar part mutually see.Especially, for system embodiment, because it is substantially similar to embodiment of the method, so description is fairly simple, relevant part illustrates see the part of embodiment of the method.
Those skilled in the art can also recognize the various illustrative components, blocks (illustrativelogicalblock) that the embodiment of the present invention is listed, unit, and step can pass through electronic hardware, computer software, or both combinations realize.For the replaceability (interchangeability) of clear displaying hardware and software, above-mentioned various illustrative components (illustrativecomponents), unit and step have universally described their function.Such function is the designing requirement realizing depending on specific application and whole system by hardware or software.Those skilled in the art for often kind of specifically application, can use the function described in the realization of various method, but this realization can should not be understood to the scope exceeding embodiment of the present invention protection.
Various illustrative logical block described in the embodiment of the present invention, or unit can pass through general processor, digital signal processor, special IC (ASIC), field programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or the design of above-mentioned any combination realizes or operates described function.General processor can be microprocessor, and alternatively, this general processor also can be any traditional processor, controller, microcontroller or state machine.Processor also can be realized by the combination of calculation element, such as digital signal processor and microprocessor, multi-microprocessor, and a Digital Signal Processor Core combined by one or more microprocessor, or other similar configuration any realizes.
The software module that method described in the embodiment of the present invention or the step of algorithm directly can embed hardware, processor performs or the combination of both.Software module can be stored in the storage medium of other arbitrary form in RAM storer, flash memory, ROM storer, eprom memory, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area.Exemplarily, storage medium can be connected with processor, with make processor can from storage medium reading information, and write information can be deposited to storage medium.Alternatively, storage medium can also be integrated in processor.Processor and storage medium can be arranged in ASIC, and ASIC can be arranged in user terminal.Alternatively, processor and storage medium also can be arranged in the different parts in user terminal.
In one or more exemplary design, the above-mentioned functions described by the embodiment of the present invention can realize in the combination in any of hardware, software, firmware or this three.If realized in software, these functions can store on the medium with computer-readable, or are transmitted on the medium of computer-readable with one or more instruction or code form.Computer readable medium comprises computer storage medium and is convenient to make to allow computer program transfer to the telecommunication media in other place from a place.Storage medium can be that any general or special computer can the useable medium of access.Such as, such computer readable media can include but not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other anyly may be used for carrying or store the medium that can be read the program code of form with instruction or data structure and other by general or special computer or general or special processor.In addition, any connection can be properly termed computer readable medium, such as, if software is by a concentric cable, fiber optic cables, twisted-pair feeder, Digital Subscriber Line (DSL) or being also comprised in defined computer readable medium with wireless way for transmittings such as such as infrared, wireless and microwaves from a web-site, server or other remote resource.Described video disc (disk) and disk (disc) comprise Zip disk, radium-shine dish, CD, DVD, floppy disk and Blu-ray Disc, and disk is usually with magnetic duplication data, and video disc carries out optical reproduction data with laser usually.Above-mentioned combination also can be included in computer readable medium.
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.
Claims (12)
1. a method for Crack Detection, is characterized in that, described method comprises:
Obtain the multi-dimensional earthquake data in target work area;
According to the analyzing and processing to described multi-dimensional earthquake data, determine the principal stress field direction in described target work area;
According to described principal stress field direction, described multi-dimensional earthquake data are arranged composition orientation multi-dimensional earthquake data according to position angle;
Ellipse fitting is carried out in seismologic parameter distribution according to described point of orientation multi-dimensional earthquake data set, obtains growth direction and the growth strength in crack.
2. the method for a kind of Crack Detection according to claim 1, is characterized in that, described multi-dimensional earthquake data comprise based on longitudinal survey line, horizontal survey line, degree of depth survey line, offset distance, azimuthal five dimension geological datas.
3. the method for a kind of Crack Detection according to claim 1, is characterized in that, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
The distribution of the sector divided in adjustment point orientation multi-dimensional earthquake data, until described principal stress field direction is arranged in complete sector.
4. the method for a kind of Crack Detection according to claim 1, is characterized in that, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
According to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the degree of covering of the sector of different directions matches in described point of orientation multi-dimensional earthquake data.
5. the method for a kind of Crack Detection according to claim 1, is characterized in that, described according to described principal stress field direction by described multi-dimensional earthquake data according to position angle arrangement composition orientation multi-dimensional earthquake data, comprising:
The number of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the signal to noise ratio (S/N ratio) of geological data and the degree of covering of described sector match in described sector.
6. the method for a kind of Crack Detection according to claim 1, is characterized in that, described seismologic parameter comprise following in any one: amplitude, frequency, when walking.
7. a device for Crack Detection, is characterized in that, described device comprises:
Multi-dimensional earthquake data capture unit, for obtaining the multi-dimensional earthquake data in target work area;
Principal stress field direction-determining unit, for according to the analyzing and processing to described multi-dimensional earthquake data, determines the principal stress field direction in described target work area;
Divide orientation multi-dimensional earthquake data ordering unit, for described multi-dimensional earthquake data being arranged composition orientation multi-dimensional earthquake data according to position angle according to described principal stress field direction;
Ellipse fitting unit, carries out ellipse fitting for the seismologic parameter distribution according to described point of orientation multi-dimensional earthquake data set, obtains growth direction and the growth strength in crack.
8. the device of a kind of Crack Detection according to claim 7, is characterized in that, described multi-dimensional earthquake data comprise based on longitudinal survey line, horizontal survey line, degree of depth survey line, offset distance, azimuthal five dimension geological datas.
9. the device of a kind of Crack Detection according to claim 7, is characterized in that, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector distribution adjustment unit, for adjusting the distribution of the sector divided in point orientation multi-dimensional earthquake data, until described principal stress field direction is arranged in complete sector.
10. the device of a kind of Crack Detection according to claim 7, is characterized in that, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector-size adjustment unit, for according to the described degree of covering of multi-dimensional earthquake data and the distribution situation of offset distance, the size of the sector divided in adjustment point orientation multi-dimensional earthquake data, until the degree of covering of the sector of different directions matches in described point of orientation multi-dimensional earthquake data.
The device of 11. a kind of Crack Detection according to claim 7, is characterized in that, described point of orientation multi-dimensional earthquake data ordering unit, comprising:
Sector number adjustment unit, for adjusting the number of the sector divided in point orientation multi-dimensional earthquake data, until the signal to noise ratio (S/N ratio) of geological data and the degree of covering of described sector match in described sector.
The device of 12. a kind of Crack Detection according to claim 7, is characterized in that, described seismologic parameter comprise following in any one: amplitude, frequency, when walking.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510726803.XA CN105403917B (en) | 2015-10-30 | 2015-10-30 | A kind of method and device of Crack Detection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510726803.XA CN105403917B (en) | 2015-10-30 | 2015-10-30 | A kind of method and device of Crack Detection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105403917A true CN105403917A (en) | 2016-03-16 |
| CN105403917B CN105403917B (en) | 2018-10-16 |
Family
ID=55469514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510726803.XA Active CN105403917B (en) | 2015-10-30 | 2015-10-30 | A kind of method and device of Crack Detection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105403917B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105954797A (en) * | 2016-04-19 | 2016-09-21 | 中国石油天然气集团公司 | Fracture identification method and fracture identification device of seismic data |
| CN107045143A (en) * | 2017-02-22 | 2017-08-15 | 中国石油天然气股份有限公司 | A kind of method and device of Prediction of fracture |
| CN109581485A (en) * | 2018-12-04 | 2019-04-05 | 成都捷科思石油天然气技术发展有限公司 | A method of carrying out automatic slit detection directly on pre-stack depth migration seismic data |
| CN111505715A (en) * | 2020-04-26 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack parameters based on central incidence angle of depth domain |
| CN112394408A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Anisotropic medium crack prediction method and device |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030018436A1 (en) * | 2001-07-20 | 2003-01-23 | Stark Tracy Joseph | System for multi-dimensional data analysis |
| CN103076623A (en) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | Crack detection method based on prestack coherence |
-
2015
- 2015-10-30 CN CN201510726803.XA patent/CN105403917B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030018436A1 (en) * | 2001-07-20 | 2003-01-23 | Stark Tracy Joseph | System for multi-dimensional data analysis |
| CN103076623A (en) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | Crack detection method based on prestack coherence |
Non-Patent Citations (4)
| Title |
|---|
| 曲良超 等: "西峰油田白马中区长8段储层裂缝发育特点及水淹预测", 《石油地质与工程》 * |
| 李相文 等: "GeoEast-EasyTrack在缝洞型碳酸盐岩研究中的应用", 《内蒙古石油化工》 * |
| 王洪求 等: "不同地震属性的方位各向异性分析及裂缝预测", 《石油地球物理勘探》 * |
| 赵伦 等: "复杂碳酸盐岩储集层裂缝发育特征及形成机制—以哈萨克斯坦让纳若尔油田为例", 《石油勘探与开发》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105954797A (en) * | 2016-04-19 | 2016-09-21 | 中国石油天然气集团公司 | Fracture identification method and fracture identification device of seismic data |
| CN105954797B (en) * | 2016-04-19 | 2018-05-04 | 中国石油天然气集团公司 | The fault recognition method and apparatus of seismic data |
| CN107045143A (en) * | 2017-02-22 | 2017-08-15 | 中国石油天然气股份有限公司 | A kind of method and device of Prediction of fracture |
| CN107045143B (en) * | 2017-02-22 | 2018-12-28 | 中国石油天然气股份有限公司 | A kind of method and device of Prediction of fracture |
| CN109581485A (en) * | 2018-12-04 | 2019-04-05 | 成都捷科思石油天然气技术发展有限公司 | A method of carrying out automatic slit detection directly on pre-stack depth migration seismic data |
| CN111505715A (en) * | 2020-04-26 | 2020-08-07 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack parameters based on central incidence angle of depth domain |
| CN112394408A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Anisotropic medium crack prediction method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105403917B (en) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104330828B (en) | Dessert reservoir forecasting method and forecasting device | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| CN103109207B (en) | For detecting the method for subsurface seismic events in vertical transverse isotropic medium | |
| US20170205531A1 (en) | Geological modeling workflow | |
| US10215868B2 (en) | Automated horizon auto-picking on multiple volumes | |
| US10571584B2 (en) | Global inversion based estimation of anisotropy parameters for orthorhombic media | |
| CN104793247B (en) | Near-surface velocity modeling method and model building device | |
| US20130100770A1 (en) | Method for Imaging Microseismic Events Using an Azimuthally-dependent Focal Mechanism | |
| CN105403917A (en) | Crack detection method and device | |
| CN107045143B (en) | A kind of method and device of Prediction of fracture | |
| CN105759312A (en) | Reservoir characteristic method well-seismic calibration method | |
| CN105277976A (en) | Earthquake forward modelling method based on rock outcrop radar detection | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| Grechka et al. | Microseismic at HFTS2: A story of three stimulated wells | |
| CN105301638B (en) | A kind of method and apparatus in extraction base of weathering face | |
| CN108646292A (en) | Fracture spacing prediction technique, device and computer storage media | |
| CN109469477B (en) | Method and device for predicting extension direction of artificial crack | |
| Roy et al. | Integrated characterization of hydraulic fracture treatments in the Barnett Shale: The Stocker geophysical experiment | |
| CN104267434B (en) | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device | |
| CN109298448A (en) | A kind of prediction technique and device of tight gas fracturing engineering dessert | |
| CN104614762A (en) | Loose sandstone gas reservoir boundary determination method and device | |
| Lanzarone et al. | The value of VSP data through early phases of field appraisal and development: A modeling and acquisition case study in the Gulf of Mexico | |
| US11474272B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
| Bellatreche et al. | Detailed seismic lithofluid distribution using Bayesian stochastic inversion for a thinly bedded reservoir: A case study over Huntington UK Central North Sea | |
| Le Calvez et al. | Combining surface and downhole microseismic monitoring: a multi-well completion in Oman |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |