CN106294900A - A kind of structure and methods for using them of the digital rock mass for drilling well - Google Patents
A kind of structure and methods for using them of the digital rock mass for drilling well Download PDFInfo
- Publication number
- CN106294900A CN106294900A CN201510276067.2A CN201510276067A CN106294900A CN 106294900 A CN106294900 A CN 106294900A CN 201510276067 A CN201510276067 A CN 201510276067A CN 106294900 A CN106294900 A CN 106294900A
- Authority
- CN
- China
- Prior art keywords
- grid
- layer
- summit
- attribute
- value
- 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
- 238000005553 drilling Methods 0.000 title claims abstract description 56
- 239000011435 rock Substances 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000010276 construction Methods 0.000 claims abstract description 13
- 241001269238 Data Species 0.000 claims abstract description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 10
- 238000005516 engineering process Methods 0.000 claims description 7
- 239000011148 porous material Substances 0.000 claims description 7
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims 1
- 235000011613 Pinus brutia Nutrition 0.000 claims 1
- 241000018646 Pinus brutia Species 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 6
- 241001074085 Scophthalmus aquosus Species 0.000 abstract description 3
- 230000008901 benefit Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000000205 computational method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000035800 maturation Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a kind of digital rock mass for drilling well and build and methods for using them, this construction method includes: set up the stratigraphic framework model of setting regions；Obtain the drilling engineering geological property value on wellbore trace；Calculate the variogram that every attribute is distributed in each layer of position；In each layer of position, being respectively directed to every attribute and travel through each grid, if there is this attribute in grid, then assigning it to each summit of this grid；If each summit of grid is all or part of sky, then utilizes variogram to be calculated property value, and be assigned to the summit of this null value；The all vertex datas obtained in each layer of position are stored and obtains numeral rock mass.During the Drilling Design or real brill of new well, use these geologic parameters, then can obtain each vertex parameter of this coordinate place grid by the coordinate of each point on well track, and then calculate every kind of engineering parameter of this coordinate points, finally can obtain all engineering parameters of full hole track.
Description
Technical field
This method belongs to petroleum drilling technology field, particularly relate to a kind of digital rock mass for drilling well structure and
Application process.
Background technology
In petroleum exploration and development process, the geological condition of underground must be studied by we, and just can have puts
Vow, and take reasonably to drill, mining type, exploration success ratio and the rate of oil and gas recovery could be improved.
For many years, for exploration and exploitation, each oil field consumption huge fund establishes the oil reservoir of seismic data cube and reservoir and retouches
State data volume etc., In Oil Field Exploration And Development serves pivotal role.But, in drilling engineering field,
Needing also exist for block as unit, the various geologic(al) factors that drilling well be must take into integrate, it is simple to optimize and bore
Well scheme, reduces construction risk and improves the purpose of drilling efficiency.
At present, some investigators proposes the computational methods of some drilling engineering geologic parameters, but does not proposes to have
Close the construction method of numeral rock mass (i.e. engineering geology data volume).
Summary of the invention
One of the technical problem to be solved is to need to provide the structure of a kind of digital rock mass for drilling well
Construction method.
In order to solve above-mentioned technical problem, embodiments herein provide firstly a kind of numeral for drilling well
The construction method of rock mass, including: set up the stratigraphic framework model of setting regions, described stratigraphic framework model every
The space of individual layer position is divided into Corner-point Grids form, and the summit cartesian coordinate of each grid represents；Obtain
Various drilling engineering geological property values on described setting regions wellbore trace；Calculate every attribute described
The variogram of distribution in each layer of position of layer Skeleton Model；In each layer of position of described stratigraphic framework model,
Being respectively directed to every attribute and travel through each grid, if there is this attribute in grid, then this attribute being composed simultaneously
To each summit of this grid；Travel through each grid, if each summit of grid is all or part of null value,
The variogram then utilizing the attribute corresponding to this null value to be distributed in this layer of position is calculated property value, and by it
It is assigned to the summit of this null value；The all vertex datas storage that will obtain in each layer of position of described stratigraphic framework model
Get up to obtain numeral rock mass.
Preferably, in setting up the step of stratigraphic framework model of setting regions, with every in described setting regions
The geological layering that a bite history borehole logging tool is explained is foundation, utilizes spatial interpolation technology to set up the space on each stratum
Geometry screen work.
Preferably, described drilling engineering geological property includes: comprcssive strength, tensile strength, cohesive strength, interior
Angle of friction, Poisson's ratio, elastic modelling quantity, bulk modulus, maximum horizontal principal stress, minimum horizontal principal stress, hang down
Normal stress, formation pore pressure, formation collapsed pressure and/or formation fracture pressure.
Preferably, utilize the variogram that attribute corresponding to null value is distributed in this layer of position, use Ke Lijin to insert
Value method is calculated property value, and this property value is assigned to the summit of null value.
Preferably, all tops that Octree mode will obtain in each layer of position of described stratigraphic framework model are utilized
Point data stores and obtains numeral rock mass.
Preferably, geostatistics method is utilized to calculate every attribute each layer at described stratigraphic framework model
The variogram of distribution in position.
According to a further aspect in the invention, additionally provide the application process of a kind of digital rock mass for drilling well,
Including: well track is divided into multiple tracing point, it is thus achieved that the degree of depth of each tracing point, hole angle and azimuth,
And it is calculated the cartesian coordinate value of each tracing point；For each tracing point, corresponding Descartes is utilized to sit
Scale value, searches for the grid residing for this coordinate from described numeral rock mass, and is calculated each brill of this tracing point
Well engineering geology property value.
Preferably, in the step of each drilling engineering geological property value being calculated this tracing point, if
This tracing point overlaps with arbitrary summit of certain grid, then the property value on this summit is directly assigned to this tracing point；
Otherwise, by calculating the arithmetic mean of instantaneous value of each summit respective attributes value of grid residing for this tracing point, this rail is obtained
Each drilling engineering geological property value of mark point.
Compared with prior art, the one or more embodiments in such scheme can have the advantage that or have
Benefit effect.
Embodiment of the present invention stratigraphic framework model based on setting regions, build this setting regions for drilling well
Digital rock mass, and the various engineering geology parameters along well track can be obtained from numeral rock mass, thus finely
Optimize wellbore construction scheme, improve drilling efficiency, reduce construction risk.
Other features and advantages of the present invention will illustrate in the following description, and, partly from explanation
Book becomes apparent, or understands by implementing technical scheme.The purpose of the present invention and its
He can come real by structure specifically noted in description, claims and accompanying drawing and/or flow process at advantage
Now and obtain.
Accompanying drawing explanation
Accompanying drawing is used for providing being further appreciated by of the technical scheme to the application or prior art, and constitutes
A part for bright book.Wherein, the accompanying drawing expressing the embodiment of the present application is used for explaining together with embodiments herein
The technical scheme of the application, but it is not intended that the restriction to technical scheme.
Fig. 1 is the schematic flow sheet of the construction method of the digital rock mass for drilling well of the embodiment of the present invention.
Fig. 2 is the schematic diagram of rock mass Skeleton Model.
Fig. 3 be resolution be the Corner-point Grids body schematic diagram of 3*3*3.
Fig. 4 is that tracing point A is through the schematic diagram of a grid on Corner-point Grids body.
Fig. 5 is the storage organization figure of Octree form.
Fig. 6 is the schematic flow sheet of the application process of the digital rock mass for drilling well of the embodiment of the present invention.
Fig. 7 is the schematic diagram of data pick-up.
Detailed description of the invention
Embodiments of the present invention are described in detail, whereby to the present invention how below with reference to drawings and Examples
Application technology means solve technical problem, and the process that realizes reaching relevant art effect can fully understand and evidence
To implement.Each feature in the embodiment of the present application and embodiment, can mutually tie under not colliding premise
Closing, the technical scheme formed is all within protection scope of the present invention.
It addition, the step shown in the flow chart of accompanying drawing can be at the computer of such as one group of computer executable instructions
System performs.And, although show logical order in flow charts, but in some cases, permissible
To be different from the step shown or described by order execution herein.
(embodiment one)
Fig. 1 is the schematic flow sheet of the construction method of the digital rock mass for drilling well of the embodiment of the present invention, below
Each step of the digital rock mass construction method for drilling well of the present embodiment is described with reference to Fig. 1.So-called " number
Word rock mass " mean the spatial data body of each attribute of rock mass " describe " (industry has " engineering rock mass " vocabulary,
Refer to physics rock mass).
Step S110, sets up stratigraphic framework model (the also referred to as rock mass Skeleton Model) V of setting regionsf。
It should be noted that stratigraphic framework model is by folding that coordinate data, individual-layer data and layer data are set up
Close FEM layer model.Will the identical layer group of each well couple together by isochronic correlation, cambium layer surface model, then profit
With layer data, tomography is combined with FEM layer model, sets up reservoir space screen work, and carry out three dimensional network block
Change.
Specifically when setting up stratigraphic framework model, with a bite history borehole logging tool every in setting regions and/or logging explanation
Geological layering be foundation, utilize spatial interpolation technology to set up the space geometry screen work on each stratum, i.e. obtain every
The horizontal spread of individual bed boundary.According to Fig. 2, it can be seen that the Skeleton Model constructed is by the rock mass in this region
It is divided into multilamellar.
In determining setting regions every a bite history borehole logging tool explain geological layering time, utilize seismic data cube,
History log analysis data and geological logging data interpretation obtain being layered the information such as lithology, fissure fault, utilize this
A little information together constitute stratigraphic framework model with strata division.
Step S120, obtains the various drilling engineering geological property values on this setting regions wellbore trace.
Specifically, utilize the existing computational methods of industry, binding tests data, respectively obtain each on wellbore trace
Plant drilling engineering geological property value (being referred to as x0).
More specifically, calculate first with log datas such as acoustic logging, density log and GRs, and tie
Close core experiment, obtain the formation rock property parameters (elastic and Intensity attribute parameter) along well track distribution,
This formation rock property parameters includes comprcssive strength, tensile strength, cohesive strength, internal friction angle, Poisson's ratio, bullet
Property modulus and bulk modulus etc..
Utilize the log data of each depth point of every mouthful of well (preferably each two point is spaced 1 meter) of this setting regions,
Calculating the rock properties parameter of this point respectively, finally give is that every attribute all can corresponding a series of take
Value, each depth point on corresponding well track.
The result that individual well calculates can represent such as following table:
Wherein, tracing point X-coordinate, Y coordinate and the Z coordinate of each depth point represent, by above-mentioned meter
Calculation method obtains the rock properties for each tracing point.
Secondly, for every mouthful of well, utilize industry maturation algorithm to calculate maximum horizontal principal stress, minimum water respectively
The mechanical attributes such as flat principal stress, vertical stress, formation pore pressure, formation collapsed pressure, formation fracture pressure.
More specifically, utilize density log data integral and calculating to obtain the overlying rock pressure along well track distribution
Power, then utilizes log data to combine triaxial stress test, and obtaining the maximum horizontal master along well track distribution should
Power, minimum level crustal stress and vertical stress.Utilize sound wave measuring well curve, engineering logging data or earthquake layer speed
Degrees of data calculates the formation pore pressure data along well track distribution.Log data is utilized to calculate along well
The formation collapsed pressure data of track distribution.Utilize log data to calculate and combine burst pressure tests and obtain along well
The formation fracture pressure data of eye track distribution.
By above-mentioned steps can obtain comprcssive strength, tensile strength, cohesive strength, internal friction angle, Poisson's ratio,
Elastic modelling quantity, bulk modulus, maximum horizontal principal stress, minimum horizontal principal stress, vertical stress, formation pore
These ten several drilling engineering geological property of pressure, formation collapsed pressure, formation fracture pressure.
However, it is necessary to explanation, above-mentioned drilling engineering geological property can be selected according to practical situation or needs
Calculate to selecting property, the combination of multiple drilling engineering geological property can be obtained, be not limited to that a kind of group above-mentioned
Close.
Step S130, calculates every attribute at stratigraphic framework model VfEach layer of position in distribution variogram.
Specifically, at stratigraphic framework model VfEach layer of position in, for each drilling engineering geological property, profit
Set up every kind of engineering geology attribute respectively by geostatistics method to divide in each layer of position of stratigraphic framework model
The variogram of cloth.
Variogram is a kind of tolerance of regionalized variable Spatial Variability, reflects spatial variability degree with distance
The feature of change.Variogram emphasizes the Data Shaping on three dimensions, thus can describe compartmentalization quantitatively and become
The spatial coherence of amount, the rock parameter that i.e. geologic rule is caused dependency spatially and randomness.
Typically prepare (data decimation, quality examination, data statistics) through data during concrete calculating, experiment becomes
The calculating of difference function (utilizing the variogram of the serial variance that well data calculates), theoretical variogram
Optimal fitting and structure fit (select suitable variogram model that serial experiment variogram is fitted,
Obtain the continuous variogram curve of optimum；Variogram model includes spherical model, exponential model, Gauss
Model, power function model, cavity effect model etc., can select according to actual count feature), variogram ginseng
The four big steps such as the optimality inspection of number.
Step S140, by stratigraphic framework model VfThe space of each layer of position be divided into Corner-point Grids form (this
Grid described in is 3 D stereo grid, lower same).
In the present embodiment, the floor projection of each grid soffit XY is square, and the direction of Z is walked with stratum
To vertically.Fig. 3 be resolution be the Corner-point Grids body schematic diagram of 3*3*3, it can be seen that pin
To stratigraphic framework model VfEach layer of position, divide the space into the Corner-point Grids body that resolution is 3*3*3, often
Each summit (totally 8 summits) of individual grid represents with cartesian coordinate respectively.
It should be noted that the single grid that becomes, multiple grids composition grid body, 3*3*3 net such as above
Lattice body comprises 27 grids.It addition, the grid body resolution of above-mentioned division is only an example, it is also possible to
The resolution of grid body is reasonably set according to computer attribute, thickness and required precision.Easy to understand, if
In the space of a particular size, the biggest then data of resolution are the most intensive, and precision is the highest, then can requirement therewith
Calculating and the display performance of computer are the highest.
Generate space screen work network model time, need to arrange I, J direction meshes number, arrange net boundary,
Two-dimensional grid I, the Trendline etc. in J direction are set.At I, J direction can define the size of grid cell, use
The quality of grid is improved in Trendline and direction, and screen work grid has been become fault block by tomography with boundary segmentation, and each breaks
Block has a given grid cell number, thus it is possible to vary this number is with local cypher or takes out dilute grid.?
Eventually on the basis of FAULT MODEL, according to combine fault information two-dimensional grid be distributed, set up end face, middle
With space lattice face, 3, bottom surface, just constituting a screen work grid model, formation beds can be inserted into afterwards
Wherein.
Step S150, at stratigraphic framework model VfEach layer of position in, be respectively directed to each drilling engineering above-mentioned
Geological property travels through each grid, if there is attribute x0 in grid, then this attribute x0 is assigned to this net simultaneously
Each summit (totally 8 summits) of lattice.
In other words, successively at stratigraphic framework model VfEach layer of position in, travel through each grid, if in grid
There is tracing point (depth point of well), then all drilling engineering geological property values of this tracing point are assigned to this net
Each summit on 8 summits of lattice.
For example, as shown in Figure 4, A point is a tracing point on well track, and it falls within this grid
In, then all properties value of this point is assigned to simultaneously each summit on 8 summits of grid.In the present embodiment,
Due to the comprcssive strength of A point, tensile strength, cohesive strength, internal friction angle, Poisson's ratio, elastic modelling quantity, volume
Modulus, overburden pressure, maximum horizontal principal stress, minimum horizontal principal stress, vertical stress, formation pore
These 14 drilling engineering geological property of pressure, formation collapsed pressure, formation fracture pressure are all determined, then
Then these 14 drilling engineering geological property values are assigned on each summit of A point place grid respectively.
Step S160, travels through each grid, if the 8 of grid summits are all or part of null value and (i.e. do not have
It is endowed x0), then the variogram utilizing the attribute corresponding to this null value to be distributed in this layer of position is calculated
Property value, and assign it to the summit of this null value.
Specifically, utilize the variogram that attribute corresponding to null value is distributed in this layer of position, use Kriging regression
Method is calculated property value, and this property value is assigned to the summit of null value.
If certain stratum is divided for 100 grids, there are 10 mouthfuls of offset wells (i.e. 10 tracks) this area,
The grid wherein passed by these 10 tracks assumes there are 45, also remains 55 grids and is not worn by any track
Cross, then these abortive haul lattice are accomplished by using said method to be endowed property value.
In addition, it is necessary to explanation, Kriging regression algorithm is industry maturation algorithm, and its utilization has built up
Grid screen work and corresponding variogram, can calculate the attribute data of each abortive haul lattice at random.Finally make institute
Some grid vertexes have property value.
Step S170, by stratigraphic framework model VfEach layer of position in all vertex datas of obtaining store
Obtain complete digital rock mass VE。
Store more specifically, all vertex datas are carried out tissue in Octree mode in the present embodiment.
If as it is shown in figure 5, a rock mass is divided into N number of grid, then the leaf of bottommost layer (the 1st layer)
Node is the summit (8, such as node 410～417 is 8 summits of certain grid) of each grid, meter
Calculate the attribute average property value as this grid (such as node 41) on 8 summits.Then by the most adjacent 8
Individual grid (such as node 40～47) is combined into a father node (such as node 4), calculates this 8 grid property
The meansigma methods of value is as the property value of father node, until complete for the combination of all original meshes, forming the 2nd layer.
Combine the most layer by layer, until all collecting a root node position, it is assumed that finally have M shell.With this
Plant Octree form can varigrained data be stored in data base.
The advantage of octotree data structure is: data extract carry out three-dimensional rendering time, can setting data flexibly
The number of plies that reads thus obtain varigrained display result, be conducive to adapting to different application demand and different
Computing power；Or it is optimized on display algorithm, such as first starts to read and render from root node, by
Layer renders out, allows user can see render process from coarse to fine rather than waits as long for.
Certainly, in addition to using Octree form to carry out data storage, it is also possible to use raster data structure etc.
Grid (is not combined, but all grid positions and property value is carried out coequally by other data structures
Record), therefore the present invention does not limit concrete memory module.
The method of the present embodiment, based on stratigraphic framework model, utilizes part to survey logging explanation data and experimental data,
Building the drilling engineering geological property body in region, i.e. numeral rock mass, these attributes include that comprcssive strength, tension are strong
Degree, cohesive strength, internal friction angle, Poisson's ratio, elastic modelling quantity, bulk modulus, maximum horizontal principal stress, minimum
Flatly stress, vertical stress, formation pore pressure, fracture pressure, caving pressures etc., these attributes again may be used
Realize extracting, for drilling plan optimization and risk profile along any new well track.
Further, use utilize that above-mentioned steps obtains for the digital rock mass of drilling well time, mainly at new well
When Drilling Design or construction, extract from numeral rock mass along well designed path or drilled wellbore trajectories be distributed each
Plant engineering geology parameters, for prediction or the scheme optimization of drilling well risk.The process of parameter extraction and numeral rock mass
Building process contrary, see Fig. 6, specifically include following steps.
Step S610, is divided into multiple tracing point, it is thus achieved that the degree of depth of each tracing point, hole angle by well track
And azimuth, and it is calculated the cartesian coordinate value of each tracing point.
Step S620, for each tracing point, utilizes corresponding cartesian coordinate value, from numeral rock mass VEIn
Search for the grid residing for this coordinate, be calculated each drilling engineering geological property value of this tracing point respectively.
By calculating arithmetic mean of instantaneous value (8 vertex value of each summit respective attributes value of grid residing for this tracing point
Arithmetic mean of instantaneous value), obtain each drilling engineering geological property value of this tracing point.
More specifically, in this step, if this tracing point overlaps with arbitrary summit of certain grid, then should
The property value on summit is directly assigned to this tracing point；Otherwise, by calculating each summit phase of grid residing for this tracing point
Answer the arithmetic mean of instantaneous value of property value, obtain each drilling engineering geological property value of this tracing point.
Owing to using Octree mode that all vertex datas carry out tissue storage, therefore, this number of acquisition before
Word rock mass have the advantage that is exactly, not high enough in computing power or that rock mass is described required precision is the highest
In the case of, during data pick-up, (above-mentioned rock mass has M shell altogether, and retrieval is to the 1st can suitably to reduce retrieval level
Grid corresponding to Ceng is the thinnest, obtains data precision the highest, and data volume is maximum) to realize data weight reducing.
For example, as shown in Figure 7, it is assumed that certain numeral rock mass has 64 grids, the 100m of well track
Have passed through wherein 4 grids (234/323/322/411) of this rock mass, then these 4 grids to 103 meters of sections
4 depth points that attribute data is assigned on well track respectively.
The present invention devises a kind of digital rock mass for drilling well and builds and methods for using them, surveys first with traditional
Visit exploitation Geological Modeling, utilize geological data and survey logging explanation data construct to go out tectonic geology model, so
Rear utilization drilling well, as constraints, constructs stratigraphic framework model；Each layer of stratigraphic framework model is drawn
It is divided into several 8 mesh of vertices, utilizes historical data homing method to build the change of every kind of engineering geology parameters simultaneously
Difference function；Existing method is utilized to calculate all engineering geology parameters according to the log data boring history warp afterwards
(referred to as initial parameter), utilizes corresponding variogram to extrapolate each grid body vertex position on this basis
Engineering geology parameters value, finally make the summit of all grid body all impart every kind of engineering geology parameters
Value, forms a data volume, i.e. numeral rock mass.During Drilling Design or real brill, these geology to be used
Parameter, then can obtain each vertex parameter of this coordinate place grid body by the coordinate of each point on well track,
And then calculate every kind of engineering parameter of this coordinate points, finally can obtain all engineering parameters of full hole track.
(example)
YB block is the important survey area produced in natural gas increasing the storage, but encounters a series of in this area's drilling well
A drilling well difficult problem, slimhole drilling low including bottom Continental Facies Stratigraphy conventional drilling rate of penetration raises speed difficulty greatly, again
Miscellaneous situation is multiple.For an above-mentioned difficult problem, technical staff needs flower great effort to carry out the excellent of wellbore construction scheme
Change design, it is therefore an objective to improve rate of penetration, reduce the generation of complex situations.And to carry out scheme targetedly
Optimizing, it is necessary to understand the engineering geological condition of full well section, main method before this is the situation with reference to offset well.Profit
Construct the digital rock mass of YB by the method for the invention after, often carry out the conceptual design of the new well of a bite, design
Personnel i.e. extract the various engineering geology parameters along borehole track (as lithology, stratum three are pressed from numeral rock mass
Power etc.), owing to it has merged this all offset wells of well periphery and earthquake related data, therefore these engineering geology
Parameter compared with for simply with reference to several mouthfuls of offset wells, more comprehensively and accurately, scheme optimization effect is obviously improved.
Although the embodiment that disclosed herein is as above, but described content is only to facilitate understand the present invention
And the embodiment used, it is not limited to the present invention.Technology people in any the technical field of the invention
Member, on the premise of without departing from the spirit and scope that disclosed herein, can be in the formal and details implemented
On make any amendment and change, but the scope of patent protection of the present invention, still must be with appending claims institute
Define in the range of standard.
One of ordinary skill in the art will appreciate that all or part of step realizing in above-described embodiment method is
Can instruct relevant hardware by program to complete, described program can be stored in an embodied on computer readable
In storage medium, this program upon execution, including each step of said method, described storage medium, such as:
ROM/RAM, magnetic disc, CD etc..
Claims (8)
1. it is used for a construction method for the digital rock mass of drilling well, including:
Setting up the stratigraphic framework model of setting regions, the space of each layer of position of described stratigraphic framework model is divided
For Corner-point Grids form, the summit cartesian coordinate of each grid represents；
Obtain the various drilling engineering geological property values on described setting regions wellbore trace；
Calculate the variogram that every attribute is distributed in each layer of position of described stratigraphic framework model；
In each layer of position of described stratigraphic framework model, it is respectively directed to every attribute and travels through each grid, if
There is this attribute in grid, then this attribute is assigned to each summit of this grid simultaneously；
Travel through each grid, if each summit of grid is all or part of null value, then utilize this null value institute right
The variogram that the attribute answered is distributed in this layer of position is calculated property value, and assigns it to the top of this null value
Point；
The all vertex datas obtained in each layer of position of described stratigraphic framework model are stored and obtains numeral
Rock mass.
Method the most according to claim 1, it is characterised in that setting up the Stratigraphic framework of setting regions
In the step of model,
The geological layering that every a bite history borehole logging tool is explained in described setting regions, as foundation, utilizes space interpolation
Technology sets up the space geometry screen work on each stratum.
Method the most according to claim 1 and 2, it is characterised in that
Described drilling engineering geological property includes: comprcssive strength, tensile strength, cohesive strength, internal friction angle, pool
Pine ratio, elastic modelling quantity, bulk modulus, maximum horizontal principal stress, minimum horizontal principal stress, vertical stress,
Layer pore pressure, formation collapsed pressure and/or formation fracture pressure.
The most according to the method in any one of claims 1 to 3, it is characterised in that
Utilize the variogram that attribute corresponding to null value is distributed in this layer of position, use Kriging regression method to calculate
To property value, and this property value is assigned to the summit of null value.
Method the most according to any one of claim 1 to 4, it is characterised in that
Octree mode is utilized to be deposited by all vertex datas obtained in each layer of position of described stratigraphic framework model
Storage gets up to obtain numeral rock mass.
Method the most according to any one of claim 1 to 5, it is characterised in that
Utilize geostatistics method to calculate every attribute to be distributed in each layer of position of described stratigraphic framework model
Variogram.
7. it is used for an application process for the digital rock mass of drilling well, including:
Well track is divided into multiple tracing point, it is thus achieved that the degree of depth of each tracing point, hole angle and azimuth,
And it is calculated the cartesian coordinate value of each tracing point；
For each tracing point, utilize corresponding cartesian coordinate value, from described numeral rock mass, search for this coordinate
Residing grid, and it is calculated each drilling engineering geological property value of this tracing point.
Method the most according to claim 7, it is characterised in that be calculated each of this tracing point
In the step of drilling engineering geological property value,
If this tracing point overlaps with arbitrary summit of certain grid, then the property value on this summit is directly assigned to this
Tracing point；Otherwise, by calculating the arithmetic mean of instantaneous value of each summit respective attributes value of grid residing for this tracing point,
Obtain each drilling engineering geological property value of this tracing point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510276067.2A CN106294900A (en) | 2015-05-26 | 2015-05-26 | A kind of structure and methods for using them of the digital rock mass for drilling well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510276067.2A CN106294900A (en) | 2015-05-26 | 2015-05-26 | A kind of structure and methods for using them of the digital rock mass for drilling well |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106294900A true CN106294900A (en) | 2017-01-04 |
Family
ID=57634634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510276067.2A Pending CN106294900A (en) | 2015-05-26 | 2015-05-26 | A kind of structure and methods for using them of the digital rock mass for drilling well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106294900A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109902360A (en) * | 2018-08-16 | 2019-06-18 | 清能艾科（深圳）能源技术有限公司 | Method, device and machine equipment for optimizing engineering parameters in drilling site operations |
| CN113139263A (en) * | 2020-01-17 | 2021-07-20 | 中国石油化工股份有限公司 | Method and system for improving drilling simulation calculation speed and accuracy |
| CN116152460A (en) * | 2023-04-14 | 2023-05-23 | 瞳见科技有限公司 | Method, device, terminal and medium for generating rock stratum model based on UE4 |
| CN119227382A (en) * | 2024-09-25 | 2024-12-31 | 西南石油大学 | A method for optimizing wellbore structure based on three-dimensional formation pressure |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726754A (en) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | Image simulation method based on continuous multiple-point geostatistics method, soft data and hard data |
| CN102313790A (en) * | 2011-07-19 | 2012-01-11 | 北京师范大学 | Submarine geologic body carbon dioxide sequestration potential assessment method |
| CN104182571A (en) * | 2014-08-12 | 2014-12-03 | 电子科技大学 | Kriging interpolation method based on Delaunay and GPU |
| CN104252549A (en) * | 2013-06-28 | 2014-12-31 | 中国石油天然气股份有限公司 | Analytical well spacing method based on kriging interpolation |
-
2015
- 2015-05-26 CN CN201510276067.2A patent/CN106294900A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726754A (en) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | Image simulation method based on continuous multiple-point geostatistics method, soft data and hard data |
| CN102313790A (en) * | 2011-07-19 | 2012-01-11 | 北京师范大学 | Submarine geologic body carbon dioxide sequestration potential assessment method |
| CN104252549A (en) * | 2013-06-28 | 2014-12-31 | 中国石油天然气股份有限公司 | Analytical well spacing method based on kriging interpolation |
| CN104182571A (en) * | 2014-08-12 | 2014-12-03 | 电子科技大学 | Kriging interpolation method based on Delaunay and GPU |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109902360A (en) * | 2018-08-16 | 2019-06-18 | 清能艾科（深圳）能源技术有限公司 | Method, device and machine equipment for optimizing engineering parameters in drilling site operations |
| CN109902360B (en) * | 2018-08-16 | 2022-12-09 | 清能艾科（深圳）能源技术有限公司 | Method, device and mechanical equipment for optimizing engineering parameters in drilling field operation |
| CN113139263A (en) * | 2020-01-17 | 2021-07-20 | 中国石油化工股份有限公司 | Method and system for improving drilling simulation calculation speed and accuracy |
| CN116152460A (en) * | 2023-04-14 | 2023-05-23 | 瞳见科技有限公司 | Method, device, terminal and medium for generating rock stratum model based on UE4 |
| CN116152460B (en) * | 2023-04-14 | 2024-03-29 | 瞳见科技有限公司 | Method, device, terminal and medium for generating rock stratum model based on UE4 |
| CN119227382A (en) * | 2024-09-25 | 2024-12-31 | 西南石油大学 | A method for optimizing wellbore structure based on three-dimensional formation pressure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817535B (en) | Method and device for determining thin interlayer | |
| CN104992468B (en) | Three-dimensional geological modeling method for fracture-cave carbonate reservoir | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| CN103514630B (en) | Fault structure three-dimensional modeling method | |
| AU2005287335B2 (en) | Method for creating facies probability cubes based upon geologic interpretation | |
| CN105651676B (en) | A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern | |
| CN102222365A (en) | Method for reconstructing curved surface of complex space | |
| Harris | The role of geology in reservoir simulation studies | |
| CN103632397B (en) | Multi-block multi-unit three-dimensional geological survey method | |
| CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
| Ou et al. | 3D visualization of hydraulic fractures using micro-seismic monitoring: Methodology and application | |
| CN115469361A (en) | Three-dimensional geological modeling method for clastic rock stratum | |
| CN106294900A (en) | A kind of structure and methods for using them of the digital rock mass for drilling well | |
| CN114740547A (en) | Back-flushing sliding fracture fold grading and partitioning analysis and three-dimensional characterization method | |
| Abeed et al. | Modeling of petrophysical properties and reserve estimation of Mishrif formation-Garraf oil field | |
| Al Rassas et al. | Application of 3d reservoir geological model on es1 formation, block nv32, shenvsi oilfield, China | |
| CN114820959A (en) | Method for digitally modeling multiple geological elements of rock burst mine | |
| Gonzalez et al. | 3D modeling of multistage hydraulic fractures and two-way-coupling geomechanics/fluid-flow simulation of a horizontal well in the nikanassin tight gas formation, Western Canada sedimentary basin | |
| Al-Rubaye et al. | A high resolution 3D geomodel for giant carbonate reservoir-A field case study from an Iraqi oil field | |
| CN111859687A (en) | Mixed geological modeling method and system for depicting geological structure of uranium-bearing sand layer | |
| CN106154326A (en) | The method and device that a kind of vertical buckling fold fracture spacing is evaluated | |
| CN114690245B (en) | A forward modeling method for geological model of seismic anomalies without drilling | |
| Gaddipati et al. | An Integrated Reservoir Modeling Case Study to Simulate Multi-Stage Hydraulically Fractured Horizontal Wells, based on Seismic, Petrophysical and Geological data for Pinedale Tight Gas Fluvial Reservoir | |
| Cox et al. | Integrated modeling for optimum management of a giant gas condensate reservoir, Jurassic eolian Nugget Sandstone, Anschutz Ranch East field, Utah overthrust (USA) | |
| Zellou et al. | Integrated fractured reservoir characterization using neural networks and fuzzy logic: Three case studies |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170104 |