CN105223611B - Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set - Google Patents
Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set Download PDFInfo
- Publication number
- CN105223611B CN105223611B CN201410317648.1A CN201410317648A CN105223611B CN 105223611 B CN105223611 B CN 105223611B CN 201410317648 A CN201410317648 A CN 201410317648A CN 105223611 B CN105223611 B CN 105223611B
- Authority
- CN
- China
- Prior art keywords
- pressure break
- area
- point set
- point
- plane
- 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.)
- Active
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of method for determining pressure break reservoir leading edge SRF based on micro-seismic monitoring three-dimensional point set, calculating pressure break reservoir area SRA and pressure break reservoir volume SRV.This method first projects to whole three-dimensional data points on horizontal plane, and the Alpha shapes of plane point set are determined using the Alpha shape methods in computational geometry, and it is exactly SRA to calculate the area that SRF is surrounded.Then microseism three-dimensional point set is layered in vertical direction, according to the plane distribution of case point in layer, the SRF and SRA of each layer of calculating, the area SRA of individual layer is multiplied by floor height and obtains the volume SRV of individual layer, and all layers of volume summation obtains total SRV.This method determines SRF shapes, the area and volume of SRA and SRV quantitative description crack areas of crack area according to the case point three-dimensional coordinate data of micro-seismic monitoring, for the size of production capacity after quantitative estimation pressure break, improves the accuracy of evaluation of Fracturing Effect on Compact Sandstone.
Description
Technical field
The present invention relates to geophysics field, especially relates to one kind according to Oil/gas Well hydraulic fracturing micro-seismic monitoring knot
Fruit (i.e. three-dimensional point set) determines pressure break reservoir leading edge (english abbreviation SRF), calculates pressure break reservoir area (english abbreviation SRA) and pressure
Split the method for reservoir volume (english abbreviation SRV).
Background technology
Hydraulic fracturing is an important measures of well production increment, and the man-made fracture that hydraulic fracturing is formed can by micro-
Monitoring result is shaken to be described.Because high-pressure fluid injects reservoir during hydraulic fracturing, cause reservoir rock to rupture and discharge micro-
Weak seismic signal, these microseism signals can be received by monitoring array by underground or surface deployment, then pass through signal transacting
Focus coordinate can be determined with seismic source location.The result of hydraulic fracturing micro-seismic monitoring is the sequence of events in three dimensions
Point, the parameters such as the orientation, length, height in crack can be gone out by the fitting of distribution of case point for simple strip crack, still
For the chicken-wire cracking of complexity, the distribution of micro-seismic event point does not have clear and definite directionality, and distribution is more dispersed, Wu Fayong
The simple parameters such as orientation, length, height describe.
The content of the invention
The purpose of the present invention is in view of the shortcomings of the prior art, to propose one kind according to Oil/gas Well hydraulic fracturing micro-seismic monitoring
The three-dimensional coordinate data of result event point, particularly suitable complex fracture carry out crack intuitively and quantitative description needs based on micro-
Pressure break reservoir leading edge, area and the calculation method of physical volume of shake monitoring three-dimensional point set.
The present invention overall plan be：
Step 1, the three-dimensional coordinate data of micro-seismic monitoring result event point is projected into horizontal plane, it is total calculates crack area
The pressure break reservoir leading edge (SRF) and pressure break reservoir area (SRA) of body；
Step 2, micro-seismic monitoring result event point is divided into multilayer in vertical direction, calculates every layer of pressure break reservoir leading edge
And pressure break reservoir area (SRA) (SRF)；
Step 3, every layer of pressure break reservoir area (SRA) is multiplied by floor height, obtains every layer of pressure break reservoir volume (SRV),
All layers of pressure break reservoir volume (SRV) summation obtains overall pressure break reservoir volume (SRV).
Such scheme also includes：It is using the point set on horizontal plane, using computational geometry to calculate pressure break reservoir leading edge (SRF)
Alpha-shape methods in determine the Alpha shapes of plane event point, and the line of its boundary point means that the wheel of shape
Exterior feature, that is, obtain pressure break reservoir leading edge (SRF).
Obtained pressure break reservoir leading edge (SRF) is a triangle subdivision of plane event point, calculates each triangle
Area and summation just obtain the pressure break reservoir area (SRA) represented by pressure break reservoir leading edge (SRF) shape.
Such scheme further comprises：
In step 1, the micro-seismic monitoring result of use is a series of point set P of micro-seismic event three-dimensional focal point compositions
={ p1 p2 … pN, wherein N is the number of total case point, pi=(xi yi zi) it is i-th of case point and its three-dimensional seat
Mark, three-dimensional point set P is projected into horizontal plane and obtains corresponding plane point set Q={ q1 q2 … qN, wherein qi=(xi,yi) be
I-th of plane event point and its coordinate；
For plane point set Q, using the Alpha shapes of the Alpha-shape method construct plane point sets in computational geometry
Shape, by constantly adjusting the value of the adjustable parameter α in Alpha-shape methods, obtain a plane close to true shape
The triangle subdivision of point set, taking boundary point, line just obtains intuitively representing the pressure of fracturing fracture affected area shape successively
Split reservoir leading edge (SRF)；
Plane point set Q Alpha shapes are a triangle subdivisions of plane point set, are calculated according to following formula each
The area of triangle
Wherein,And a, b, c are respectively the length of side on Atria bar side, by plane point set Q Alpha shapes
All triangle area summations of shape, that is, obtain pressure break reservoir area (SRA)；
In step 2, three-dimensional point set P is layered in vertical direction, according to the scope of vertical coordinate, floor height, h points are M
Layer, every layer has NkIndividual case point, andThe three-dimensional point set of kth layer isProject to water
The plane point set of individual layer is obtained in planeUsing identical method in step 1, i.e., using Alpha-
Shape methods obtain individual layer plane point set QkAlpha shapes, taking boundary point, line must be before the pressure break reservoir of kth layer successively
Edge (SRF), then the pressure break reservoir area (SRA) of kth layer is calculated；
In step 3, just obtain every layer of crack multiplied by with floor height h using every layer of (SRA) that step 2 obtains and spread
The volume in region, all layers of volume summation just obtain total pressure break reservoir volume (SRV).
The SRF of all case point plane distributions of method construct in the present invention intuitively depicts fracturing fracture felt area
The flat shape in domain, the area SRA that SRF surrounds shape is calculated, can quantitatively describe the area that fracturing fracture involves；Vertical
After direction division multilayer, the SRF and SRA of every layer of case point can be calculated, in floor height and the adjustable parameter of Alpha-shape methods
When α values are preferable, overall SRV can be more accurately calculated, so as to which the volume in production capacity region after quantitative estimation pressure break is big
It is small.
Brief description of the drawings
Fig. 1 is the flow chart of a specific embodiment of the inventive method；
Fig. 2 is the three-dimensional spatial distribution figure of certain a bite well micro-seismic monitoring result event point；
Fig. 3 is the distribution for the plane point set that three-dimensional events point projects in horizontal plane；
Fig. 4 is the Alpha shapes of the plane point set of all micro-seismic event points；
Fig. 5 is the Alpha shapes of the plane point set of a certain layer；
Fig. 6 is that the SRA of a certain layer is multiplied by the volume schematic diagram that floor height obtains.
Embodiment
To enable objects, features and advantages of the present invention to become apparent, it is cited below particularly go out preferred embodiment, and coordinate institute
Accompanying drawings, it is described in detail below.
The first step, micro-seismic monitoring result are a series of point set P={ p of micro-seismic event three-dimensional focal point compositions1 p2
… pN, wherein N is the number of total case point, pi=(xi yi zi) it is i-th of case point and its three-dimensional coordinate.By three-dimensional point
Collection P projects to horizontal plane and obtains corresponding plane point set Q={ q1 q2 … qN, wherein qi=(xi,yi) it is i-th of plane thing
Part point and its coordinate.
For plane point set Q, using the Alpha shapes of the Alpha-shape method construct plane point sets in computational geometry
Shape.Alpha-shape methods are a kind of methods of the Delaunay triangle subdivisions based on plane point set, Delaunay triangles
The point set shape contour that subdivision method obtains can only be convex polygon, and Alpha-shape methods can obtain concave polygon
Shape, the shape of expression are more accurate.Alpha-shape methods have an adjustable parameter α, can by continuous adjusting parameter α value
To obtain the profile close to true shape.A triangle subdivision that Alpha-shape methods obtain and plane point set,
Take boundary point successively line must pressure break reservoir leading edge SRF, SRF intuitively illustrate the shaped wheel of fracturing fracture affected area
It is wide.
Plane point set Q Alpha shapes are a triangle subdivisions of plane point set, are calculated according to following formula each
The area of triangle
Wherein,And a, b, c are respectively the length of side on Atria bar side, by plane point set Q Alpha shapes
The area summation of all triangles of shape, that is, obtain pressure break reservoir area SRA, SRA and quantitatively illustrate fracturing fracture felt area
The size in domain.
Second step, three-dimensional point set P is layered in vertical direction, according to the scope of vertical coordinate, by certain floor height point
For multilayer, such as floor height h rice, it is divided into M layers, every layer has NkIndividual case point, andThe three-dimensional point set of kth layer isProject on horizontal plane and obtain the plane point set of individual layerWith step 1
In it is identical, individual layer plane point set Q is obtained using Alpha-shape methodskAlpha shapes, taking boundary point, line must successively
The pressure break reservoir leading edge SRF of kth layer, then the pressure break reservoir area SRA of kth layer is calculated
3rd step, using every layer of the SRA that step 2 obtains multiplied by the crack diffusion zone for floor height h, just obtaining every layer
Volume, all layers of volume summation just obtain total pressure break reservoir volume SRV.
The flow and effect of the present invention are described further with concrete application example below in conjunction with the accompanying drawings.
As shown in figure 1, Fig. 1 is SRF, SRA and SRV calculating side based on micro-seismic monitoring result three-dimensional point set of the present invention
The flow chart of method.
In a step 101, according to the three-dimensional coordinate data of micro-seismic monitoring result event point, the three-dimensional space of case point is drawn
Between distribution map, such as Fig. 2.
In a step 102, by micro-seismic monitoring result event spot projection to horizontal plane, the distribution of case point two dimensional surface is drawn
Figure, such as Fig. 3.
In step 103, for the plane point set of micro-seismic event point, using Alpha-shape method construct planar points
The Alpha shapes of collection.Fig. 4 is the Alpha shapes of constructed plane point set, and asterisk is the micro-seismic event point in plane, is put down
The Alpha shapes of millet cake collection are a subsets of the Delaunay triangle subdivisions of plane point set, the triangle of internal subdivision
While being represented by dashed line, the line of the marginal point of Alpha shapes shows the profile of shape, is exactly pressure break reservoir leading edge SRF.
At step 104, the Alpha shapes of the plane point set obtained according to step 103, the face of each triangle is calculated
Product, all triangle areas and be exactly pressure break reservoir area SRA.
In a step 101, the vertical range being distributed according to three-dimensional events point, case point is divided into vertical direction more
Layer, into step 201.
In step 201, be M layers by three-dimensional point set P points, it is assumed that floor height h, then calculate every layer of case point SRF,
SRA and SRV.
In step 202, using the Alpha shapes with step 103 identical method construct individual layer case point.Fig. 5 is it
In a certain layer plane event point and its Alpha- shapes, the line of border incident point forms the SRF of individual layer.
In step 203, the triangle subdivision of the individual layer case point Alpha shapes determined according to step 202, calculate all
The area of triangle and summation, obtain the SRA of individual layer.
In step 204, the individual layer SRA obtained according to step 203, multiplied by with the floor height of individual layer, that is, individual layer is obtained
SRV。
In step 205, all layers of individual layer SRV is added, just obtains total pressure break reservoir volume SRV.Total SRV amounts
Illustrate the volume that fracturing fracture involves reservoir with changing, be the important evidence for estimating production capacity size.
Claims (3)
1. pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set, it is characterised in that the party
Method includes：
Step 1, the three-dimensional coordinate data of micro-seismic monitoring result event point is projected into horizontal plane, calculates crack area totality
Pressure break reservoir leading edge (SRF) and pressure break reservoir area (SRA)；
Step 2, micro-seismic monitoring result event point is divided into multilayer in vertical direction, calculates every layer of pressure break reservoir leading edge
And pressure break reservoir area (SRA) (SRF)；
Step 3, every layer of pressure break reservoir area (SRA) is multiplied by floor height, obtains every layer of pressure break reservoir volume (SRV), owned
Pressure break reservoir volume (SRV) summation of layer obtains overall pressure break reservoir volume (SRV)；
It is using the point set on horizontal plane, using the Alpha-shape in computational geometry to calculate pressure break reservoir leading edge (SRF)
Method determines the Alpha shapes of plane event point, and the line of its boundary point means that the profile of shape, that is, obtains pressure break reservoir
Leading edge (SRF).
2. the pressure break reservoir leading edge according to claim 1 based on micro-seismic monitoring three-dimensional point set, area and volume calculate
Method, it is characterised in that：Obtained pressure break reservoir leading edge (SRF) is a triangle subdivision of plane event point, is calculated each
The area of triangle and summation just obtain the pressure break reservoir area (SRA) represented by pressure break reservoir leading edge (SRF) shape.
3. the pressure break reservoir leading edge based on micro-seismic monitoring three-dimensional point set, area and volume according to claim 2 calculate
Method, it is characterised in that：
In step 1, the micro-seismic monitoring result of use be a series of micro-seismic event three-dimensional focal points compositions point set P=
{p1 p2 … pN, wherein N is the number of total case point, pi=(xi yi zi) it is i-th case point and its three-dimensional
Coordinate, three-dimensional point set P is projected into horizontal plane and obtains corresponding plane point set Q={ q1 q2 … qN, wherein qi=
(xi,yi) it is i-th plane event point and its coordinate；
For plane point set Q, using the Alpha shapes of the Alpha-shape method construct plane point sets in computational geometry,
By constantly adjusting the value of the adjustable parameter α in Alpha-shape methods, a planar point close to true shape is obtained
The triangle subdivision of collection, taking boundary point, line just obtains intuitively representing the pressure break of fracturing fracture affected area shape successively
Reservoir leading edge (SRF)；
Plane point set Q Alpha shapes are a triangle subdivisions of plane point set, and each three are calculated according to following formula
Angular area
Wherein,And a, b, c are respectively the length of side on Atria bar side, by plane point set Q's
All triangle area summations of Alpha shapes, that is, obtain pressure break reservoir area (SRA)；
In step 2, three-dimensional point set P is layered in vertical direction, according to the scope of vertical coordinate, floor height, h points are M layers,
Every layer has Nk case point, andThe three-dimensional point set of kth layer isProjection
The plane point set of individual layer is obtained on to horizontal planeUsing identical method in step 1,
Individual layer plane point set Qk Alpha shapes are obtained using Alpha-shape methods, taking boundary point, line must kth successively
The pressure break reservoir leading edge (SRF) of layer, then the pressure break reservoir area (SRA) of kth layer is calculated；
In step 3, every layer of the pressure break reservoir area (SRA) obtained using step 2 just obtains every layer multiplied by with floor height h
Crack diffusion zone volume, all layers volume summation just obtains total pressure break reservoir volume (SRV).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410317648.1A CN105223611B (en) | 2014-07-04 | 2014-07-04 | Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410317648.1A CN105223611B (en) | 2014-07-04 | 2014-07-04 | Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105223611A CN105223611A (en) | 2016-01-06 |
| CN105223611B true CN105223611B (en) | 2018-02-02 |
Family
ID=54992666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410317648.1A Active CN105223611B (en) | 2014-07-04 | 2014-07-04 | Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105223611B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109507723A (en) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | The pressure break calculation method of physical volume and system of microseism fracturing fracture model |
| CN109001814A (en) * | 2018-08-07 | 2018-12-14 | 大庆时代宏业石油科技有限公司 | Underground micro-seismic monitoring method |
| CN111553504B (en) * | 2019-01-23 | 2022-06-03 | 中国石油化工股份有限公司 | Method and system for calculating volume of fracturing modified oil and gas reservoir |
| CN110288227A (en) * | 2019-06-24 | 2019-09-27 | 中国石油大学(北京) | A method of fracturing effect Dominated Factors are influenced for evaluating |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606124A (en) * | 2012-03-16 | 2012-07-25 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | Sectioned multi-cluster synchronous volume fracturing technology for double horizontal wells |
| CN103760602A (en) * | 2014-01-27 | 2014-04-30 | 中国石油天然气集团公司 | Method and system for forwards simulating three-dimensional hydrofracture and microseism event |
| CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8902710B2 (en) * | 2009-11-10 | 2014-12-02 | Microseismic, Inc. | Method for determining discrete fracture networks from passive seismic signals and its application to subsurface reservoir simulation |
-
2014
- 2014-07-04 CN CN201410317648.1A patent/CN105223611B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606124A (en) * | 2012-03-16 | 2012-07-25 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | Sectioned multi-cluster synchronous volume fracturing technology for double horizontal wells |
| CN103760602A (en) * | 2014-01-27 | 2014-04-30 | 中国石油天然气集团公司 | Method and system for forwards simulating three-dimensional hydrofracture and microseism event |
| CN103869363A (en) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | Micro-earthquake positioning method and device |
Non-Patent Citations (1)
| Title |
|---|
| "What Is Stimulated Reservior Volume";MJ.Mayerhofer el;《SPE Shale Gas Production Conference》;20100228;第89页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105223611A (en) | 2016-01-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105223611B (en) | Pressure break reservoir leading edge, area and calculation method of physical volume based on micro-seismic monitoring three-dimensional point set | |
| CN110147561B (en) | Method for predicting volume fracture network of tight oil and gas reservoir containing natural fracture | |
| CN103759677B (en) | Based on the method for triangle area ratio law amount rock mass discontinuity three-dimensional roughness | |
| CN105386756B (en) | A method of brittle formation porosity is calculated using dependent variable | |
| CN104483705B (en) | A kind of three-dimensional residual static corrections | |
| Jacquemyn et al. | Multi-scale three-dimensional distribution of fracture-and igneous intrusion-controlled hydrothermal dolomite from digital outcrop model, Latemar platform, Dolomites, northern Italy | |
| CN102879820B (en) | Based on the three-dimensional table layer model construction method of triangle gridding | |
| CN103400405A (en) | Multi-beam bathymetric chart construction method based on seabed digital depth model feature extraction | |
| Zawada et al. | A multiscale analysis of coral reef topographic complexity using lidar-derived bathymetry | |
| CN103437381B (en) | Clump of piles traction superposition settling amount evaluation process method | |
| CN103559368B (en) | Chemically mechanical polishing analogy method and clearance computational methods thereof | |
| CN104933208B (en) | A kind of different scale Fractured Reservoir numerical simulation integrated conduct method | |
| Coscarella et al. | Turbulence anisotropy and intermittency in open-channel flows on rough beds | |
| CN106033653A (en) | Ground-laser-scanning-based three-dimensional indoor map making method | |
| CN110068399A (en) | Reconstruction of temperature field algorithm based on radial basis function and regularization | |
| CN102778689B (en) | Wide curved line seismic data underground reflection line building method | |
| CN109884698B (en) | Quantitative evaluation method of seismic exploration observation system based on target layer | |
| CN105573963A (en) | Reconstruction method for horizontal nonuniform structure of ionized layer | |
| CN104361578A (en) | Hierarchical grid partition method under multi-scale precision control | |
| CN104020508A (en) | Direct ray tracking algorithm for geological radar chromatography detecting | |
| CN102736104A (en) | Near-surface surface layer model construction method | |
| JP5671184B1 (en) | Tomographic map creating apparatus, tomographic map creating method and tomographic map creating program | |
| CN104391272B (en) | The method and system of target positioning is carried out using direction finding data | |
| CN105223630B (en) | Omnibearing observation systematic parameter Demonstration Method based on geological model | |
| Langousis et al. | Theoretical framework to estimate spatial rainfall averages conditional on river discharges and point rainfall measurements from a single location: an application to western Greece |
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 |