CN103886129B - By the discrete method and apparatus to reservoir grid model of log data - Google Patents
By the discrete method and apparatus to reservoir grid model of log data Download PDFInfo
- Publication number
- CN103886129B CN103886129B CN201410057529.7A CN201410057529A CN103886129B CN 103886129 B CN103886129 B CN 103886129B CN 201410057529 A CN201410057529 A CN 201410057529A CN 103886129 B CN103886129 B CN 103886129B
- Authority
- CN
- China
- Prior art keywords
- grid cell
- point
- grid
- well track
- well
- 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
Abstract
The present invention provides a kind of by the discrete method and apparatus to reservoir grid model of log data, is related to geological property modeling technique field.The reservoir grid model is Corner-point Grids, and the method comprises the following steps：The grid cell that S1, acquisition well track are passed through；S2, the intersection point of the well track with the grid cell is calculated, and the well track is divided into multiple blocks using the intersection point；S3, it is weighted averagely to falling into the log property value on the well track of each block, tries to achieve attribute average, and the attribute average is assigned to corresponding grid cell.According to log data value, the grid assignment passed through by well track with the time as few as possible, so as to reduce the stand-by period of modeling personnel, improves operating efficiency to the present invention.
Description
Technical field
The present invention relates to geological property modeling technique field, and in particular to a kind of that log data is discrete to reservoir grid mould
The method and apparatus of type.
Background technology
Parameter of the geologic body inside containing characteristics such as various reflection formation lithologies, petroleum resources distributions, such as porosity, infiltration
Rate, containing saturation degree, deposit equal, the modeling of geologic body physical property is exactly utilization space interpolation method, such as Kriging, sequence Gauss,
Solve geologic body internal physical, chemical attribute parameter model.Achievement is modeled using geologic body physical property, storage can be effectively found
Layer, and Accurate Prediction is carried out to petroleum resources distribution situation, so as to instruct oil field development.In view of the stratum origin cause of formation, in geologic body
Some region of physical property, with upper strata or the lower floor that the relation of same layer physical property is far longer than the region.For more accurate expression
Geologic body physical property, in actual modeling process, on the basis of to geologic body subdivision substratum, can further segment, after subdivision
Each layer, using the physical property infomation on this layer obtained by borehole log data, carries out physical property and builds using spatial interpolation algorithm
Mould.By the discrete reservoir grid cell assignment for being to utilize log data to reservoir grid model, being well is passed through of log data.Assigned
The grid cell of value, when physical property modeling is carried out, as known quantity, and then determines other not measured property values for crossing region.
It can be seen that, it is the basis for carrying out geologic body physical property modeling.
At present, Corner-point Grids are the grid models of the description reservoir being widely used.As structured grid model, it
Have the advantages that to be easily achieved.The flexibility of height allows it to characterize reservoir geology feature well, such as heterogeneous
Property and tomography.
However, because Corner-point Grids do not use cartesian coordinate system, it is extremely complex cumbersome in discrete log data.Again plus
On in order to meet the demand to modeling region fine description, present mesh generation is all very fine and closely woven.So, it is a kind of rapidly
The discrete method to reservoir model of log data is by exigence.
The content of the invention
（One）The technical problem of solution
In view of the shortcomings of the prior art, the invention provides a kind of by the discrete method to reservoir grid model of log data
And device, according to log data value, the grid assignment passed through by well track with the time as few as possible, so as to reduce modeler
The stand-by period of member, improve operating efficiency.
（Two）Technical scheme
To realize object above, the present invention is achieved by the following technical programs：
It is a kind of by the discrete method to reservoir grid model of log data, the reservoir grid model is Corner-point Grids, bag
Include the following steps：
The grid cell that S1, acquisition well track are passed through；
S2, the intersection point of the well track with the grid cell is calculated, and divided the well track using the intersection point
It is multiple blocks；
S3, it is weighted averagely to falling into the log property value on the well track of each block, tries to achieve attribute
Average, and the attribute average is assigned to corresponding grid cell.
Preferably, step S1 includes：
S11, the interval that tested point selection in well track is set according to the web thickness information of reservoir grid model；
S12, according to it is described interval sequentially take out tested point；
S13, find each tested point where grid cell；
S14, the grid cell is arranged and merged.
Preferably, step S11 includes：
S111, according to well head coordinate, and reservoir grid model boundary information and grid length and width configuration information, estimate
Grid cell below well head；
S112, the average thickness for calculating well head lower section grid cell；
In the sampling interval of S113, surge well track data, calculate the interval of tested point selection；
Wherein, step S13 includes：Correspondence grid cell is found out using reference key, centered on grid cell, using facing
Nearly search method chooses grid cell, and inclusive test is carried out to tested point and the grid cell chosen；
Wherein, step S14 includes：By the index information of grid cell, using binary search, find out what is be not measured
Grid cell；For the situation that multiple tested points fall into same grid cell, only record first falls into treating for the grid cell
Measuring point information, the tested point information of the grid cell is fallen into last.
Preferably, the well track is calculated in step S2 includes step with the intersection point of the grid cell：
According to the face for having intersection point in grid cell index information determining unit with well track；The tested point at selection face two ends,
Calculate the intersection point of the well track and grid cell face between two tested points.
Preferably, quick approximatioss is used during the intersection point for calculating well track and grid cell face between two tested points, is realized
Mode is：
2 adjacent well track data points composition straight line is chosen, is calculated with the intersection point on grid cell surface, calculate the intersection point
With the distance of data point；
The equal data point of selected distance on well track, and on the basis of this data point, chooses adjacent thereto next
Individual well track data point, using the 2 adjacent well track data points for retrieving, carry out next time find intersection computing；
When intersection point is located between the 2 adjacent well track data points chosen, terminate computing.
Preferably, when well track being divided into multiple blocks in step S2, by two adjacent wells and grid cell face
Region between intersection point, as a block, a grid cell of each block correspondence reservoir grid model.
Preferably, when calculating falls into the log property value average of each block in step S3, unique numerical value is obtained,
And using this numerical value as the numerical value of block correspondence grid cell.
Present invention also offers a kind of by the discrete device to reservoir grid model of log data, including：
Grid cell searching modul, for finding out the grid cell that well track is passed through,
Log data division module, for calculating intersection point of the well track with grid cell face, and to log data dividing regions
Block；
Attribute computing module, for calculating weighted average to all log property values for falling into each block, and
It is corresponding grid cell assignment；
Interactive module, conducts interviews to reservoir grid model and log data.
Preferably, the grid cell searching modul includes：
Tested point interval generator, for producing tested point to choose interval；
Point sampling collector to be measured, the sampling collection for choosing tested point according to appointed interval；
Grid cell searcher, for the test equipment and net of the grid cell residing for accurate, quick judgement tested point
Lattice unit searches；
Determinant in grid cell, for judging whether tested point is located in grid cell；
Wherein, the log data division module includes：
Intersection point solver, for the quick intersection point for calculating well track and grid cell face；
Log data block device, for dividing block to log data；
Wherein, the attribute computing module includes：
Attribute solver, the log attribute weight average value of same block is fallen into for calculating；
Attribute assignor, for block attributes value to be set into correspondence grid cell property value；
Wherein, the interactive module includes：
Reservoir model accesses submodule, for the three-dimensional index value according to reservoir grid model, obtains reservoir grid model
In specifiable lattice unit information, and by specified physical property property value set to correspondence grid cell in；
Log data accesses submodule, and the access for providing well type of gesture data and log categorical data connects
Mouthful.
Preferably, the log data accesses submodule according to sampled point in corresponding well type of gesture data and log
Index in categorical data, obtains the customizing messages of sampled point；Wherein described well track data is included：Fathom information, big
D coordinates value and vertical depth under ground coordinate system；The borehole log data includes：Fathom information, porosity value,
Natural gamma value and sedimentary facies phase value.
（Three）Beneficial effect
The invention provides a kind of by the discrete method and apparatus to reservoir grid model of log data, by obtaining well rail
The grid cell that mark is passed through, calculates intersection point of the well track with the grid cell, and utilize the intersection point by the well rail
Mark is divided into multiple blocks；And then to fall into the log property value on the well track of each block be weighted it is flat
, attribute average is tried to achieve, and the attribute average is assigned to corresponding grid cell.
The present invention can efficiently calculate the category of the grid cell passed through by well in reservoir grid model using log data
Property value, reduces the stand-by period of modeling personnel, there is provided operating efficiency.
The present invention can accurately calculate the property value of corresponding units in reservoir grid model using log data, overcome
It is difficult to calculate the shortcoming of locus under Corner-point Grids system, improves the accuracy of modeling result.Considerably increase oil storage
The confidence level of layer fine description system, for the work such as oil well harvesting, oil field secondary development provide reliable support, so as to significantly save
About oilfield exploitation cost.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is by the schematic flow sheet of the discrete method to reservoir grid model of log data in the embodiment of the present invention；
Fig. 2 is the schematic diagram of bilinear surface；
Fig. 3 is by the structural representation of the discrete device to reservoir grid model of log data in the embodiment of the present invention；
Fig. 4 is the implementing procedure figure of grid cell searching modul；
Fig. 5 is the schematic diagram for calculating Corner-point Grids cell-average thickness；
Fig. 6 is the schematic diagram for calculating ray and bilinear surface intersection method；
Fig. 7 is to judge whether tested point is located at the schematic diagram for specifying grid cell method；
Fig. 8 is the implementing procedure figure of log data division module；
Fig. 9 is the schematic diagram that well track and grid cell face intersection point process are solved with quick approximatioss；
Figure 10 is the implementing procedure figure of attribute computing module；
Figure 11 is that the method that grid cell property value is determined when a certain block is fallen into without borehole log data point is illustrated
Figure.
Specific embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Embodiment 1：
As shown in figure 1, the embodiment of the invention provides a kind of by the discrete method to reservoir grid model of log data, institute
Reservoir grid model is stated for Corner-point Grids, is comprised the following steps：
The grid cell that S1, acquisition well track are passed through；
S2, the intersection point of the well track with the grid cell is calculated, and divided the well track using the intersection point
It is multiple blocks；
S3, it is weighted averagely to falling into the log property value on the well track of each block, tries to achieve attribute
Average, and the attribute average is assigned to corresponding grid cell.
The embodiment of the present invention can efficiently calculate the grid list passed through by well in reservoir grid model using log data
The property value of unit, reduces the stand-by period of modeling personnel, there is provided operating efficiency.Accurately storage can be calculated using log data
The property value of corresponding units, overcomes the shortcoming for being difficult to calculate locus under Corner-point Grids system in layer grid model, improves
The accuracy of modeling result.Considerably increase the confidence level of petroleum reservoir fine description system, be oil well harvesting, oil field it is secondary
The work such as exploitation provide reliable support, so as to significantly save oilfield exploitation cost.
The embodiment of the present invention is described in detail below：
Preferably, step S1 includes：
S11, the interval that tested point selection in well track is set according to the web thickness information of reservoir grid model；
S12, according to it is described interval sequentially take out tested point；
S13, find each tested point where grid cell；
S14, the grid cell is arranged and merged.
Preferably, step S11 includes：
S111, according to well head coordinate, and reservoir grid model boundary information and grid length and width configuration information, estimate
Grid cell below well head；
S112, the average thickness for calculating well head lower section grid cell；
In the sampling interval of S113, surge well track data, calculate the interval of tested point selection.
Preferably, step S13 includes：Correspondence grid cell is found out using reference key, centered on grid cell, is used
Close on search method and choose grid cell, inclusive test is carried out to tested point and the grid cell chosen.
Preferably, step S14 includes：
By the index information of grid cell, using binary search, the grid cell not being measured is found out；For multiple
Tested point falls into the situation of same grid cell, and only record first falls into the tested point information of the grid cell, with last
The individual tested point information for falling into the grid cell.
Preferably, the well track is calculated in step S2 includes step with the intersection point of the grid cell：
According to the face for having intersection point in grid cell index information determining unit with well track；The tested point at selection face two ends,
Calculate the intersection point of the well track and grid cell face between two tested points.
Preferably, quick approximatioss is used during the intersection point for calculating well track and grid cell face between two tested points, is realized
Mode is：
2 adjacent well track data points composition straight line is chosen, is calculated with the intersection point on grid cell surface, calculate the intersection point
With the distance of data point, the equal data point of selected distance on well track, on the basis of this data point, selection is adjacent thereto
Next well track data point.Using the 2 adjacent well track data points for retrieving, carry out next time find intersection computing.Directly
It is located between the 2 adjacent well track data points chosen to intersection point, algorithm terminates.
Preferably, when well track being divided into multiple blocks in step S2, by two adjacent wells and grid cell face
Region between intersection point, as a block, a grid cell of each block correspondence reservoir grid model.
Preferably, when calculating falls into the log property value average of each block in step S3, unique numerical value is obtained,
And using this numerical value as the numerical value of block correspondence grid cell.
Embodiment 2：
As shown in figure 3, present invention also offers a kind of by the discrete device to reservoir grid model of log data, including：
Grid cell searching modul, for finding out the grid cell that well track is passed through,
Log data division module, for calculating intersection point of the well track with grid cell face, and to log data dividing regions
Block；
Attribute computing module, for calculating weighted average to all log property values for falling into each block, and
It is corresponding grid cell assignment；
Interactive module, conducts interviews to reservoir grid model and log data.
Preferably, the grid cell searching modul includes：
Tested point interval generator, for producing tested point to choose interval；
Point sampling collector to be measured, the sampling collection for choosing tested point according to appointed interval；
Grid cell searcher, for the test equipment and net of the grid cell residing for accurate, quick judgement tested point
Lattice unit searches；
Determinant in grid cell, for judging whether tested point is located in grid cell.
Preferably, the log data division module includes：
Intersection point solver, for the quick intersection point for calculating well track and grid cell face；
Log data block device, for dividing block to log data.
Preferably, the attribute computing module includes：
Attribute solver, the log attribute weight average value of same block is fallen into for calculating；
Attribute assignor, for block attributes value to be set into correspondence grid cell property value.
Preferably, the interactive module includes：
Reservoir model accesses submodule, for the three-dimensional index value according to reservoir grid model, obtains reservoir grid model
In specifiable lattice unit information, and by specified physical property property value set to correspondence grid cell in；
Log data accesses submodule, and the access for providing well type of gesture data and log categorical data connects
Mouthful.
Preferably, the log data accesses submodule according to sampled point in corresponding well type of gesture data and log
Index in categorical data, obtains the customizing messages of sampled point；Wherein described well track data is included：Fathom information, big
D coordinates value and vertical depth under ground coordinate system；The borehole log data includes：Fathom information, porosity value,
Natural gamma value and sedimentary facies phase value.
The specific work process to each module is described below：
The implementing procedure figure of grid cell searching modul, as shown in figure 4, comprising the steps of：
1）The mouth coordinate information of discrete well logging is treated in acquisition, the net specified during using the boundary information of reservoir model with modeling
Lattice obtain reference index information in x directions, the increment in y directions.Reference index information is meant that the storage obtained according to the index
Layer model grid cell is located at or the approximate underface positioned at well head；
2）Using corresponding grid cell in reference key acquisition reservoir, grid cell average thickness in a z-direction is calculated
Degree zSpan；By taking Fig. 5 as an example, computing formula selection：
zi0、zi1It is respectively the upper angle point of grid cell, the depth value of lower angle point in a grid posts
The sampling interval t of well track data is obtained, and is spaced using the selection that below equation calculates tested point：
3）With s as step-length, tested point is chosen to the leapfrog of big value since well track data fathoms small value；
4）Its grid cell being located at is judged each tested point；
Judge a little whether in grid cell when, use improved ray intersection statistical method.Ray asks friendship statistics side
Method principle is：A ray is drawn from tested point, the ray is calculated with polyhedral intersection point number.If odd number, then the point
In polyhedron；If even number, then the point is external positioned at multiaspect.As shown in Fig. 2 being the schematic diagram of bilinear surface.This hair
The bright improvement to ray intersection statistical method based on the fact grid cell face is bilinear surface, when intersection is carried out
The computing formula of friendship is sought using ray and bilinear surface, it is public rather than asking for the ray and plane used on ordinary meaning
Formula.Ray and bilinear surface ask the friendship derivation of equation as follows：
As shown in fig. 6, the parametric form method for expressing of ray is：
Bilinear surface is by four angle pointsWhat interpolation was obtained, in mapping point space,
Parameter u, v is used to represent that its linear equation is：
By two formula simultaneous, three t of unknown parameter, the equation of u, v are obtained：
4 angle points are substituted into respectivelyAnd pointX, y, z coordinate value, formed three equations
Group, obtains the value of t, and then obtain intersecting point coordinate by quadratic equation with one unknown is solved after the unit that disappears.
Determination methods of the point in grid cell as shown in fig. 7,It is tested point, fromDraw one vertically upward
Ray, calculates the intersection point number of ray and grid cell face.
Specific implementation is：
（1）As shown in Fig. 7 (b), respectively willProjected on x0y faces with the one of face of grid, it is rightCarry out a little
Judgement in polygon：
IfIn the quadrangle that grid surface border determines, then continue to judge ray whether there is intersection point with the grid surface, if
There is intersection point, then intersection point number adds 1；
IfOn quadrangle side and not in 4 corner locations, then continue to judge ray whether there is friendship with the grid surface
Point, if there is intersection point, intersection point number adds 1/2；
IfPositioned at 4 corner locations of quadrangle, then continue to judge ray whether there is intersection point with the grid surface, if in the presence of
Intersection point, then intersection point number add 1/3；
IfOutside quadrangle, then skip.
（2）Statistics intersection point sum, if odd number is then in point in grid cell；If even number is then put outside grid cell.
When grid cell is chosen, using closing on search strategy.One grid cell is obtained by reference to index.And with this
Centered on grid cell, the grid cell chosen in cubical area nearby judges successively.If tested point is located in one of net
In lattice unit, terminate the search of grid cell.If tested point is not in arbitrary mess, the nearest grid list of selected distance tested point
Unit is reference.The above method is repeated, the grid cell that tested point falls into is found.If the grid of reference unit for using twice is same
Individual, tested point is without falling into any one grid cell.Judgement to the tested point terminates.
Calculating tested point mentioned above is implemented as to the method for grid cell distance, calculates 8 angle points of grid cell
The weighted average of D coordinates value, the position of grid cell is represented with this, obtains distance of this with tested point.
Whenever a tested point falls into grid cell, the positional information of the grid cell index value and tested point is recorded,
Index information in well track.
5）The grid cell obtained in the way of to choosing tested point by interval is arranged.
（1）Continuity check and modification are carried out to all grid cells passed through by well track found out, step is implemented
It is rapid as follows：
Step1:By record above order take two neighboring grid cell successively, calculate they index between Manhattan away from
From going to Step2；
Step2:If manhatton distance is more than 1, Step3 is gone to, otherwise, return to Step1；
Step3:The middle data of pair two tested points being associated with two grid cells judges, judges it at which
In grid cell, insert result into record, the grid cell sequence number for making Step1 judge retracts one, returns to Step1.
（2）Redundancy inspection and modification are carried out to all grid cells passed through by well track found out.
Step1:Take two neighboring grid cell successively by new record order, calculate they index between Manhattan away from
From going to Step2；
Step2:If manhatton distance is 0, Step3. is gone to otherwise, return to Step1；
Step3:What two grid cells were represented is a grid cell, is merged.And close latter grid cell
The tested point information of connection is associated with previous grid cell.One grid cell possesses two tested point information, and edge is represented respectively
Well track data fathom it is small be worth to the big value lookup that fathoms, first tested point for falling into the grid cell and last
The individual tested point for falling into the unit.The current grid cell sequence number for judging of order retracts one, returns to step1.
Make two indexed coordinates of grid cell respectively (i1,j1,k1),(i2,j2,k2), then Manhattan mentioned above away from
It is from computational methods：
Manhattan Distance=|i1–i2|+|j1–j2|+|k1–k2|。
The implementing procedure figure of log data division module, as shown in figure 8, comprising the steps of：
Intersection point of the well track with grid cell is calculated using log data division module, and according to above intersection point information to surveying
Well data divide block.
1）Calculate intersection point of the well track with grid cell.
The grid cell face for having intersection point with well track is determined according to two index values from grid cell.Take previous net
The tested point information of the unit is finally fallen into lattice unit, first tested point for falling into the unit in latter grid cell is taken
Information.In the interval chosen, the intersection point of well track and grid cell face is calculated using quick approximatioss.
Quick approximatioss using well track in local curvature's very little, tend to the characteristic of straight line, since tested point one end just
After step estimates intersection point, the well track sampled point line segment that there can not possibly be intersection point with grid cell face is skipped, direct selection most has can
The sampled point line segment that can there is intersection point judged, and using result of calculation as feedback, and then more accurate instructs sampling dotted line
The selection of section.
The implementation example description as shown in Figure 9 of quick approximatioss.pspeIt is the calculating border of well track, p1It is well rail
With p on marksAdjacent point.Method is initial with psp1As line segment to be calculated, line segment p is calculatedsp1With the intersection point p ' in grid cell face,
If p ' is in line segment psp1On, then find intersection point.If p ' is in line segment psp1Outward, p ' and p is calculated1Distance, found on well track away from
From p1Equal point p2, and p2Latter point p3.To line segment p2p3Residing straight line calculates new intersection point with grid cell face.
The above method is repeated, until finding intersection point.
2）Well track is divided into block, the region between the adjacent intersection point of each two, as an area using above-mentioned intersection point
Block, represents the well data that the block included and entirely falls in corresponding grid cell.
The workflow diagram of attribute computing module, as shown in Figure 10, comprises the steps of：
1）Weighted average to falling into the borehole log data point treatment aiming at the pathogenesis value of same block, obtains unique attribute
Value；
2）The property value of grid cell corresponding with each block is set using attribute assignor.
Particularly, situation as shown in figure 11, when the grid cell passed through by well track is fallen into without borehole log data,
Their midpoint is tried to achieve first by two intersection point informations on grid cell face, and tries to achieve it and fathomed；Again from profit
With the data point that the midpoint two ends are located in depth information searching borehole log data；Property value using two data points enters line
Property interpolation, tries to achieve the property value at the midpoint, and in this, as the property value of the grid cell.
To sum up, the use of the embodiment of the present invention quick approach method by the method for interval sampling and when asking friendship is significantly
Degree is reduced on calculation times of the result without influence, and what grid cell searcher was based on close on search strategy utilizes well
Spatial locality principle, also greatly reduces amount of calculation, improves calculating speed；Friendship is asked by bilinear surface and straight line
Whether algorithm allows judging a little obtain accurate analytic solutions in grid cell and when solving the point of intersection so that meter
Calculate result accuracy high.
It should be noted that herein term " including ", "comprising" or its any other variant be intended to non-row
His property is included, so that process, method, article or equipment including a series of key elements not only include those key elements, and
And also include other key elements being not expressly set out, or also include for this process, method, article or equipment institute are intrinsic
Key element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that including institute
Also there is other identical element in process, method, article or the equipment of stating key element.
The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although with reference to the foregoing embodiments
The present invention has been described in detail, it will be understood by those within the art that：It still can be to foregoing each implementation
Technical scheme described in example is modified, or carries out equivalent to which part technical characteristic；And these modification or
Replace, do not make the spirit and scope of the essence disengaging various embodiments of the present invention technical scheme of appropriate technical solution.
Claims (8)
1. it is a kind of by the discrete method to reservoir grid model of log data, it is characterised in that the reservoir grid model is angle
Dot grid, comprises the following steps：
The grid cell that S1, acquisition well track are passed through；
S2, the intersection point of the well track with the grid cell is calculated, and be divided into the well track using the intersection point many
Individual block；
S3, it is weighted averagely to falling into the log property value on the well track of each block, tries to achieve attribute average,
And the attribute average is assigned to corresponding grid cell；
Wherein, the grid cell for passing through well track that obtains includes：
S11, the interval that tested point selection in well track is set according to the web thickness information of reservoir grid model；
S12, according to it is described interval sequentially take out tested point；
S13, find each tested point where grid cell；
S14, the grid cell is arranged and merged.
2. the method for claim 1, it is characterised in that wherein, step S11 includes：
S111, according to well head coordinate, and reservoir grid model boundary information and grid length and width configuration information, estimate and be located at
Grid cell below well head；
S112, the average thickness for calculating well head lower section grid cell；
In the sampling interval of S113, surge well track data, calculate the interval of tested point selection；
Wherein, step S13 includes：Correspondence grid cell is found out using reference key, centered on grid cell, is searched using closing on
Rope method chooses grid cell, and inclusive test is carried out to tested point and the grid cell chosen；
Wherein, step S14 includes：By the index information of grid cell, using binary search, the grid not being measured is found out
Unit；For the situation that multiple tested points fall into same grid cell, first tested point for falling into the grid cell is only recorded
Information, the tested point information of the grid cell is fallen into last.
3. the method for claim 1, it is characterised in that characterized in that, calculating the same institute in the well track in step S2
The intersection point for stating grid cell includes step：
According to the face for having intersection point in grid cell index information determining unit with well track；The tested point at selection face two ends, calculates
The intersection point of well track and grid cell face between two tested points.
4. method as claimed in claim 3, it is characterised in that calculate the friendship of the well track and grid cell face between two tested points
Quick approximatioss is used during point, is achieved in that：
2 adjacent well track data points composition straight line is chosen, is calculated with the intersection point on grid cell surface, calculate the intersection point and number
The distance at strong point；
The equal data point of selected distance on well track, and on the basis of this data point, choose next well adjacent thereto
Track data point, using the 2 adjacent well track data points for retrieving, carry out next time find intersection computing；
When intersection point is located between the 2 adjacent well track data points chosen, terminate computing.
5. the method for claim 1, it is characterised in that when well track being divided into multiple blocks in step S2, by two
Region between individual adjacent well and the intersection point in grid cell face, as a block, each block correspondence reservoir grid model
A grid cell.
6. the method for claim 1, it is characterised in that the log attribute for falling into each block is calculated in step S3
During value average, unique numerical value is obtained, and using this numerical value as the numerical value of block correspondence grid cell.
7. a kind of by the discrete device to reservoir grid model of log data, it is characterised in that including：
Grid cell searching modul, for finding out the grid cell that well track is passed through,
Log data division module, for calculating intersection point of the well track with grid cell face, and divides block to log data；
Attribute computing module, for calculating weighted average to all log property values for falling into each block, and is right
The grid cell assignment answered；
Interactive module, conducts interviews to reservoir grid model and log data；
Wherein, the grid cell searching modul includes：
Tested point interval generator, for producing tested point to choose interval；
Point sampling collector to be measured, the sampling collection for choosing tested point according to appointed interval；
Grid cell searcher, for the test equipment and grid list of the grid cell residing for accurate, quick judgement tested point
Meta Search Engine；
Determinant in grid cell, for judging whether tested point is located in grid cell；
Wherein, the log data division module includes：
Intersection point solver, for the quick intersection point for calculating well track and grid cell face；
Log data block device, for dividing block to log data；
Wherein, the attribute computing module includes：
Attribute solver, the log attribute weight average value of same block is fallen into for calculating；
Attribute assignor, for block attributes value to be set into correspondence grid cell property value；
Wherein, the interactive module includes：
Reservoir model accesses submodule, for the three-dimensional index value according to reservoir grid model, in acquisition reservoir grid model
The information of specifiable lattice unit, and specified physical property property value is set in correspondence grid cell；
Log data accesses submodule, the access interface for providing well type of gesture data and log categorical data.
8. device as claimed in claim 7, it is characterised in that the log data accesses submodule according to sampled point in correspondence
Index in well type of gesture data and log categorical data, obtains the customizing messages of sampled point；Wherein described well track
Packet contains：The D coordinates value and vertical depth fathomed under information, earth coordinates；The borehole log data bag
Include：Fathom information, porosity value, natural gamma value and sedimentary facies phase value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410057529.7A CN103886129B (en) | 2014-02-19 | 2014-02-19 | By the discrete method and apparatus to reservoir grid model of log data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410057529.7A CN103886129B (en) | 2014-02-19 | 2014-02-19 | By the discrete method and apparatus to reservoir grid model of log data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103886129A CN103886129A (en) | 2014-06-25 |
| CN103886129B true CN103886129B (en) | 2017-06-13 |
Family
ID=50955020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410057529.7A Active CN103886129B (en) | 2014-02-19 | 2014-02-19 | By the discrete method and apparatus to reservoir grid model of log data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103886129B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109118557B (en) * | 2018-08-08 | 2023-06-30 | 中国石油天然气集团有限公司 | Drawing method and device for two-dimensional profile of oil well |
| CN111259465B (en) * | 2018-11-30 | 2022-09-23 | 中国石油化工股份有限公司 | Drilling process simulation method and system |
| CN111562630B (en) * | 2019-02-13 | 2023-06-20 | 中国石油化工集团有限公司 | Reservoir parameter logging evaluation method based on grid division |
| CN110334317B (en) * | 2019-07-11 | 2020-07-17 | 东北石油大学 | Automatic extraction method for picture type logging curve data |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1910620A (en) * | 2004-01-13 | 2007-02-07 | 贝克休斯公司 | 3-d visualized data set for all types of reservoir data |
| CN103400020A (en) * | 2013-04-16 | 2013-11-20 | 西南石油大学 | Numerical reservoir simulation method for measuring flowing situation of a plurality of intersected discrete cracks |
-
2014
- 2014-02-19 CN CN201410057529.7A patent/CN103886129B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1910620A (en) * | 2004-01-13 | 2007-02-07 | 贝克休斯公司 | 3-d visualized data set for all types of reservoir data |
| CN103400020A (en) * | 2013-04-16 | 2013-11-20 | 西南石油大学 | Numerical reservoir simulation method for measuring flowing situation of a plurality of intersected discrete cracks |
Non-Patent Citations (3)
| Title |
|---|
| Pettrel建模教程;百度文库;《URL：http://wenku.baidu.com/view/243bc94b336c1eb91a375d54.html》;20130117;第13-14页 * |
| Unstructured 3D gridding and upscaling for coarse modeling of geometrically complex reservoirs;Mathieu Prevost等;《2005 EAGE/Geological Society of London》;20051031;第339-345页 * |
| 水平井地质设计系统的研究与实现;李贵敏;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑 》;20120131(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103886129A (en) | 2014-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103454686B (en) | Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction | |
| CN104635262B (en) | A kind of positive reversed fault isopleth automatic generation method based on extended rectangular grid | |
| CN102222365A (en) | Method for reconstructing curved surface of complex space | |
| US10221659B2 (en) | Automated well placement for reservoir evaluation | |
| CN103886129B (en) | By the discrete method and apparatus to reservoir grid model of log data | |
| CN105934689A (en) | Geologic feature splitting | |
| CN108397130B (en) | Boring method and device | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN102222366B (en) | Method for fitting complex space curved surfaces | |
| CN104597494B (en) | Seismic strata body analysis method and device | |
| CN103969683B (en) | Method for picking position faces in batched mode based on constraint in three-dimensional seismic interpretation | |
| CN105093277A (en) | Shallow-medium-deep strata velocity fusion method in seismic modeling | |
| CN105445789A (en) | Three-dimensional Fresnel volume travel-time tomographic method based on multiple reflected refraction wave constraint | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| CN112147684A (en) | Method and device for characterizing activity intensity of co-deposition fault | |
| CN108845350A (en) | The method and device of inverting two-dimension speed model | |
| CN106338764B (en) | The super wire stripping recognition methods of bio-inspired computing stratigraphic trap oil-gas reservoir | |
| CN105301638B (en) | A kind of method and apparatus in extraction base of weathering face | |
| CN104570070B (en) | A kind of method and apparatus for setting up two-dimentional near surface geological model | |
| CN111505713A (en) | Pre-stack seismic inversion method based on multi-point geological statistics | |
| CN106777117A (en) | A kind of automatic identifying method of horizontal stratum tectonic landform | |
| CN103543478A (en) | Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method | |
| CN1292263C (en) | Ray traction in earthquake prospection | |
| CN108508481A (en) | A kind of method, apparatus and system of longitudinal wave converted wave seismic data time match | |
| CN104111476B (en) | Build the method and device of formation velocity field |
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 |