CN109598694A - Fracturing fracture calculation method and system based on Radon transformation - Google Patents
Fracturing fracture calculation method and system based on Radon transformation Download PDFInfo
- Publication number
- CN109598694A CN109598694A CN201710898148.5A CN201710898148A CN109598694A CN 109598694 A CN109598694 A CN 109598694A CN 201710898148 A CN201710898148 A CN 201710898148A CN 109598694 A CN109598694 A CN 109598694A
- Authority
- CN
- China
- Prior art keywords
- fracturing fracture
- bianry image
- transformation
- image
- converted
- 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
- 230000009466 transformation Effects 0.000 title claims abstract description 78
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 48
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 48
- 238000004364 calculation method Methods 0.000 title claims abstract description 14
- 230000004044 response Effects 0.000 claims description 14
- 238000000034 method Methods 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 6
- 230000007547 defect Effects 0.000 abstract description 5
- 238000004445 quantitative analysis Methods 0.000 abstract description 5
- 206010017076 Fracture Diseases 0.000 description 77
- 208000010392 Bone Fractures Diseases 0.000 description 76
- 238000011161 development Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Geometry (AREA)
- Probability & Statistics with Applications (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Quality & Reliability (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Disclose a kind of fracturing fracture calculation method and system based on Radon transformation.This method may include: that fracturing fracture data are converted to bianry image;Radon transformation is carried out to bianry image, obtains transformation bianry image;Based on transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.The present invention can be concisely and efficiently quantitative analysis fracturing fracture, effectively overcome the defect of conventional qualitative subjective analysis method.
Description
Technical field
The present invention relates to hydraulic fracturing microseismic field, more particularly, to a kind of based on Radon transformation
Fracturing fracture calculation method and system.
Background technique
In field of petroleum exploitation, the application of microseismic is mainly by generating in monitoring hydraulic fracturing process
Microseismic signals can monitor fracturing process, evaluate fracturing effect, and then well pattern is instructed to design and optimize engineering parameter.?
North America region, the oil reservoir that microseismic is widely used to hydraulically created fracture monitoring, high pressure injection operation generates are dynamic
The business scopes such as state monitoring, obtain the highly recognition of petroleum industrial circle, become abundant information, an essence during oil and gas development
Really, timely monitoring technology is that one of the important means of hydraulically created fracture real-time monitoring is carried out in unconventional development of resources.
Unconventional reservoir is implemented in transformation process reservoir by hydraulic fracturing, is forming one or more of major fractures
Meanwhile expand intrinsic fracture constantly and generating shearing slip with brittle rock, realize the communication to intrinsic fracture, rock bedding,
And secondary fracture is formed in the lateral pressure of major fracture, and continue branch on secondary fracture and form secondary secondary fracture.With
This analogizes, and forms intrinsic fracture and the interlaced fracture network of man-made fracture.The scale of fracture network and main development direction are
Evaluate the key index of fracturing transformation effect.
Fracturing fracture network modelling can be carried out according to the imaging of pressure break microseism focus, explanation personnel are typically based on the model
Pressure-break net qualitative interpretation is carried out, explanation results are affected by subjectivity and are not enough quantified.
Summary of the invention
The invention proposes a kind of fracturing fracture calculation methods and system based on Radon transformation, can be concisely and efficiently
Quantitative analysis fracturing fracture effectively overcomes the defect of conventional qualitative subjective analysis method.
According to an aspect of the invention, it is proposed that a kind of fracturing fracture calculation and analysis methods based on Radon transformation.It is described
Method may include: that fracturing fracture data are converted to bianry image;Radon transformation is carried out to the bianry image, is become
Change bianry image;Based on the transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
Preferably, fracturing fracture data are converted to bianry image includes: that fracturing fracture data are converted to fracturing fracture
Network;Image procossing is carried out to the fracturing fracture network, is converted by bianry image, obtains the bianry image.
Preferably, Radon transformation is carried out to the bianry image, obtaining transformation bianry image includes: with the binary map
The center of picture is polar origin, does line integral to image slices vegetarian refreshments along different angle, obtains the transformation bianry image.
It is preferably based on the transformation bianry image, the dimensional orientation and size for calculating fracturing fracture include: based on described
Bianry image is converted, judges the orientation and energy response of the fracturing fracture;According to the energy response, obtains the pressure break and split
The size of seam.
According to another aspect of the invention, it is proposed that a kind of fracturing fracture computing system based on Radon transformation, can wrap
Include: memory is stored with computer executable instructions;Processor, the computer that the processor is run in the memory can
It executes instruction, executes following steps: fracturing fracture data are converted into bianry image;Radon change is carried out to the bianry image
It changes, obtains transformation bianry image;Based on the transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
Preferably, fracturing fracture data are converted to bianry image includes: that fracturing fracture data are converted to fracturing fracture
Network;Image procossing is carried out to the fracturing fracture network, is converted by bianry image, obtains the bianry image.
Preferably, Radon transformation is carried out to the bianry image, obtaining transformation bianry image includes: with the binary map
The center of picture is polar origin, does line integral to image slices vegetarian refreshments along different angle, obtains the transformation bianry image.
It is preferably based on the transformation bianry image, the dimensional orientation and size for calculating fracturing fracture include: based on described
Bianry image is converted, judges the orientation and energy response of the fracturing fracture;According to the energy response, obtains the pressure break and split
The size of seam.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality
It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows the flow chart of the step of fracturing fracture calculation method according to the present invention based on Radon transformation.
Fig. 2 a and Fig. 2 b respectively illustrate x/y plane space according to an embodiment of the invention and Radon transformation parameter
The schematic diagram of the two-dimentional Radon shift theory in space.
Fig. 3 a and Fig. 3 b respectively illustrate bianry image according to an embodiment of the invention and convert bianry image
Schematic diagram.
Fig. 4 a and Fig. 4 b respectively illustrate fracturing fracture and it is according to an embodiment of the invention based on Radon transformation
The schematic diagram of fracturing fracture.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the flow chart of the step of fracturing fracture calculation method according to the present invention based on Radon transformation.
In this embodiment, the fracturing fracture calculation method according to the present invention based on Radon transformation may include: step
101, fracturing fracture data are converted into bianry image;Step 102, Radon transformation is carried out to bianry image, obtains transformation two-value
Image;Step 103, based on transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
In one example, fracturing fracture data are converted to bianry image includes: that fracturing fracture data are converted to pressure
Split fracture network figure;Image procossing is carried out to fracturing fracture network, is converted by bianry image, bianry image is obtained.
In one example, Radon transformation is carried out to bianry image, obtaining transformation bianry image includes: with bianry image
Center be polar origin, line integral done to image slices vegetarian refreshments along different angle, obtains transformation bianry image.
In one example, based on transformation bianry image, the dimensional orientation and size for calculating fracturing fracture include: based on change
Bianry image is changed, judges the orientation and energy response of fracturing fracture;According to energy response, the size of fracturing fracture is obtained.
Fig. 2 a and Fig. 2 b respectively illustrate x/y plane space according to an embodiment of the invention and Radon transformation parameter
The schematic diagram of the two-dimentional Radon shift theory in space.
Specifically, seam network planning mould is carried out to fracturing fracture network using Radon transformation algorithm and stitches net complicated dynamic behaviour
Theoretical basis is: under two-dimensional case, the principle of Radon transformation is in a plane along different straight line (straight line and origin
Distance is d, deflection α) line integral is done to f (x, y), obtained picture F (d, α) is exactly the Radon transformation of function f, that is,
It says, the transform value F (d, α) of each point (d, α) in the domain Radon has corresponded to some line integral value of original function, such as Fig. 2 a
With shown in Fig. 2 b.It is exactly the linear component contained inside bringing detection image that Radon, which converts one simply and directly applies,.Very
It is apparent that any straight line can all lead to extreme value of the Radon picture at the line correspondences (d, α), and in image straight line length
Degree is bigger, and the extreme value at (d, α) is also bigger.Therefore, it can use Radon transformation and form statistics carried out to pressure-break net.It will pressure
Crack data are converted to bianry image, and fracturing fracture data are converted to fracturing fracture network;To fracturing fracture network
Image procossing is carried out, is converted by bianry image, bianry image is obtained;Radon transformation is carried out to bianry image, with bianry image
Center be polar origin, line integral done to image slices vegetarian refreshments along different angle, obtains transformation bianry image;Based on transformation two
It is worth image, calculates the dimensional orientation and size of fracturing fracture.
This method can be concisely and efficiently quantitative analysis fracturing fracture, effectively overcome conventional qualitative subjective analysis method
Defect.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 3 a and Fig. 3 b respectively illustrate bianry image according to an embodiment of the invention and convert bianry image
Schematic diagram.
Fracturing fracture data are converted into bianry image, as shown in Figure 3a, fracturing fracture data are converted into fracturing fracture
Network;Image procossing is carried out to fracturing fracture network, is converted by bianry image, bianry image is obtained;To bianry image
Radon transformation is carried out to do line integral using the center of bianry image as polar origin to image slices vegetarian refreshments along different angle, obtain
Bianry image must be converted, as shown in Figure 3b, wherein abscissa is angle, θ, and ordinate d is the distance of deviation from origin;Based on transformation
Bianry image calculates the dimensional orientation and size of fracturing fracture, according to the transformation bianry image of Fig. 3 b, in Fig. 3 a development have 1. and
2. two cracks, orientation is respectively to deviate 45 degree and 135 degree of X-axis positive direction, and the energy response of crack 1. in fig 3b is big
In the energy response of crack 2., illustrate 2. the space geometry size of crack 1. is greater than crack.
Fig. 4 a and Fig. 4 b respectively illustrate fracturing fracture and it is according to an embodiment of the invention based on Radon transformation
The schematic diagram of fracturing fracture.
The ground micro-seismic monitoring record for having chosen a bite fractured well in southern china shale gas development zone, utilizes pressure break
Crack imaging technology obtains the fracturing fracture network model of the well, and as shown in fig. 4 a, white line is fracturing fracture, figure in figure
4b is transformed to Fig. 4 a progress Radon as a result, can be with the main development orientation of pressure-break net at Y-coordinate different in explanation figure 4a
With the relative size of seam network planning mould.
In conclusion the present invention can be concisely and efficiently quantitative analysis fracturing fracture, conventional qualitative subjective is effectively overcome
The defect of analysis method.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
According to an embodiment of the invention, providing a kind of fracturing fracture computing system based on Radon transformation, can wrap
Include: memory is stored with computer executable instructions;Processor, the computer executable instructions in processor run memory,
It executes following steps: fracturing fracture data is converted into bianry image;Radon transformation is carried out to bianry image, obtains transformation two
It is worth image;Based on transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
In one example, fracturing fracture data are converted to bianry image includes: that fracturing fracture data are converted to pressure
Split fracture network figure;Image procossing is carried out to fracturing fracture network, is converted by bianry image, bianry image is obtained.
In one example, Radon transformation is carried out to bianry image, obtaining transformation bianry image includes: with bianry image
Center be polar origin, line integral done to image slices vegetarian refreshments along different angle, obtains transformation bianry image.
In one example, based on transformation bianry image, the dimensional orientation and size for calculating fracturing fracture include: based on change
Bianry image is changed, judges the orientation and energy response of fracturing fracture;According to energy response, the size of fracturing fracture is obtained.
The present invention can be concisely and efficiently quantitative analysis fracturing fracture, effectively overcome conventional qualitative subjective analysis method
Defect.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.
Claims (8)
1. a kind of fracturing fracture calculation method based on Radon transformation, comprising:
Fracturing fracture data are converted into bianry image;
Radon transformation is carried out to the bianry image, obtains transformation bianry image;
Based on the transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
2. the fracturing fracture calculation method according to claim 1 based on Radon transformation, wherein by fracturing fracture data
Being converted to bianry image includes:
Fracturing fracture data are converted into fracturing fracture network;
Image procossing is carried out to the fracturing fracture network, is converted by bianry image, obtains the bianry image.
3. the fracturing fracture calculation method according to claim 1 based on Radon transformation, wherein to the bianry image
Radon transformation is carried out, obtaining transformation bianry image includes:
Using the center of the bianry image as polar origin, line integral done to image slices vegetarian refreshments along different angle, described in acquisition
Convert bianry image.
4. the fracturing fracture calculation method according to claim 1 based on Radon transformation, wherein based on the transformation two
It is worth image, the dimensional orientation and size for calculating fracturing fracture include:
Based on the transformation bianry image, the orientation and energy response of the fracturing fracture are judged;
According to the energy response, the size of the fracturing fracture is obtained.
5. a kind of fracturing fracture computing system based on Radon transformation, which is characterized in that the system includes:
Memory is stored with computer executable instructions;
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Fracturing fracture data are converted into bianry image;
Radon transformation is carried out to the bianry image, obtains transformation bianry image;
Based on the transformation bianry image, the dimensional orientation and size of fracturing fracture are calculated.
6. the fracturing fracture computing system according to claim 5 based on Radon transformation, wherein by fracturing fracture data
Being converted to bianry image includes:
Fracturing fracture data are converted into fracturing fracture network;
Image procossing is carried out to the fracturing fracture network, is converted by bianry image, obtains the bianry image.
7. the fracturing fracture computing system according to claim 5 based on Radon transformation, wherein to the bianry image
Radon transformation is carried out, obtaining transformation bianry image includes:
Using the center of the bianry image as polar origin, line integral done to image slices vegetarian refreshments along different angle, described in acquisition
Convert bianry image.
8. the fracturing fracture calculation method according to claim 5 based on Radon transformation, wherein based on the transformation two
It is worth image, the dimensional orientation and size for calculating fracturing fracture include:
Based on the transformation bianry image, the orientation and energy response of the fracturing fracture are judged;
According to the energy response, the size of the fracturing fracture is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710898148.5A CN109598694A (en) | 2017-09-28 | 2017-09-28 | Fracturing fracture calculation method and system based on Radon transformation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710898148.5A CN109598694A (en) | 2017-09-28 | 2017-09-28 | Fracturing fracture calculation method and system based on Radon transformation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109598694A true CN109598694A (en) | 2019-04-09 |
Family
ID=65955454
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710898148.5A Pending CN109598694A (en) | 2017-09-28 | 2017-09-28 | Fracturing fracture calculation method and system based on Radon transformation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109598694A (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297783A (en) * | 2014-06-30 | 2015-01-21 | 中国石油天然气集团公司 | Interpretation method and system for hydrofracture micro-seismic event |
CN105023013A (en) * | 2015-08-13 | 2015-11-04 | 西安电子科技大学 | Target detection method based on local standard deviation and Radon transformation |
-
2017
- 2017-09-28 CN CN201710898148.5A patent/CN109598694A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297783A (en) * | 2014-06-30 | 2015-01-21 | 中国石油天然气集团公司 | Interpretation method and system for hydrofracture micro-seismic event |
CN105023013A (en) * | 2015-08-13 | 2015-11-04 | 西安电子科技大学 | Target detection method based on local standard deviation and Radon transformation |
Non-Patent Citations (1)
Title |
---|
魏武: ""基于小波和Radon变换的桥梁裂缝检测"", 《计算机工程与设计》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Comprehensive evaluation of excavation-damaged zones in the deep underground caverns of the Houziyan hydropower station, Southwest China | |
Wu et al. | The pixel crack reconstruction method: From fracture image to crack geological model for fracture evolution simulation | |
US10662742B2 (en) | Method for hydrocarbon recovery with a fractal pattern and related apparatus | |
CN103775071A (en) | Method for measuring cranny evolvement of mining coal-rock mass | |
CN113687424B (en) | Carbonate fracture-cavity structure seismic characterization method based on deep learning | |
CN110672072A (en) | Information multivariate interpretation and characterization method for fractured rock mass structure | |
CN111429471A (en) | Geological disaster information management system and method | |
Lei et al. | Intelligent recognition of joints and fissures in tunnel faces using an improved mask region‐based convolutional neural network algorithm | |
CN109190218B (en) | Dynamic identification method for effective seam net of tight reservoir | |
Liu et al. | A digital twin-based framework for simulation and monitoring analysis of floating wind turbine structures | |
Zhang et al. | Deep learning‐based automatic detection of muck types for earth pressure balance shield tunneling in soft ground | |
CN116739183A (en) | Mine safety risk early warning prediction system | |
CN109298448B (en) | Prediction method and device for compact gas fracturing engineering dessert | |
KR101514623B1 (en) | Method and system for visualizing a discontinuity network of fractured rock masses in 3D | |
CN109598694A (en) | Fracturing fracture calculation method and system based on Radon transformation | |
CN117392452A (en) | Wall crack detection system and method based on image analysis | |
Youjun et al. | Analysis of pore scale fluid migration in a porous medium-application to coal rock seam | |
CN117347500A (en) | Rock fracture state identification method and related equipment | |
Han et al. | Research on a technical framework in smart construction based on big data | |
CN111881236A (en) | Stope visualization realization method and stope visualization realization device | |
CN113917562B (en) | Macro-microscopic structure representation and three-dimensional space construction method and device for deep-buried soft interlayer | |
CN110909446B (en) | Complex seam net reconstruction method based on event point and crack extension double-field coupling | |
Pollock et al. | 3D exploratory analysis of descriptive lithology records using regular expressions | |
CN115984286A (en) | Arch bridge point cloud segmentation method based on synthetic simulator and bridged neural network | |
CN111764963B (en) | Rock burst prediction method based on fast-RCNN |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190409 |
|
RJ01 | Rejection of invention patent application after publication |