CN109598694A - Fracturing fracture calculation method and system based on Radon transformation - Google Patents

Fracturing fracture calculation method and system based on Radon transformation Download PDF

Info

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
Application number
CN201710898148.5A
Other languages
Chinese (zh)
Inventor
杨心超
郭全仕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710898148.5A priority Critical patent/CN109598694A/en
Publication of CN109598694A publication Critical patent/CN109598694A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • G06T7/41Analysis of texture based on statistical description of texture
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth 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

Fracturing fracture calculation method and system based on Radon transformation
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.
CN201710898148.5A 2017-09-28 2017-09-28 Fracturing fracture calculation method and system based on Radon transformation Pending CN109598694A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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