CN107390264A - The characterizing method of the disconnected solution internal structure of carbonate rock - Google Patents
The characterizing method of the disconnected solution internal structure of carbonate rock Download PDFInfo
- Publication number
- CN107390264A CN107390264A CN201710595771.3A CN201710595771A CN107390264A CN 107390264 A CN107390264 A CN 107390264A CN 201710595771 A CN201710595771 A CN 201710595771A CN 107390264 A CN107390264 A CN 107390264A
- Authority
- CN
- China
- Prior art keywords
- tensor
- seismic data
- acoustic impedance
- attribute
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000011435 rock Substances 0.000 title claims abstract description 17
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 15
- 238000009499 grossing Methods 0.000 claims abstract description 18
- 239000002131 composite material Substances 0.000 claims description 14
- 238000005553 drilling Methods 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 11
- 238000009826 distribution Methods 0.000 abstract description 8
- 238000012512 characterization method Methods 0.000 abstract description 3
- 230000007797 corrosion Effects 0.000 description 4
- 238000005260 corrosion Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 241001074085 Scophthalmus aquosus Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003447 ipsilateral effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A kind of characterizing method for solution internal structure of breaking the invention discloses carbonate rock.Wherein method includes:Post-stack inversion is carried out to seismic data, obtains underground Acoustic Impedance Data;It is determined that representing the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development, through-carved work processing is carried out to the underground Acoustic Impedance Data according to the Acoustic Impedance Data threshold value;Calculate the tensor attribute of seismic data;Space smoothing processing is carried out to the tensor attribute, it is determined that representing the tensor attribute thresholds that disconnected solution crushes region, the result of tensor attribute space smoothing processing is handled according to the tensor attribute thresholds, obtains the profile of disconnected solution;Using the profile of the disconnected solution as border, the tensor attribute and the underground Acoustic Impedance Data are included in profile, the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.The technology that the present invention is stacked using figure carries out the comprehensive characterization of the disconnected solution of carbonate rock, can portray the border of disconnected solution oil reservoir and the spatial distribution of Reservoir Body.
Description
Technical field
The present invention relates to the production seismics method in oil reservoir development, and in particular to a kind of disconnected solution internal structure of carbonate rock
Characterizing method.
Background technology
Disconnected solution oil reservoir reservoir space is based on large-scale cave or corrosion hole, by karst water and fault development degree shadow
Ring, breaking, solution oil reservoir internal structure is extremely complex, and vertically and horizontally anisotropism is extremely strong for reservoir space, and vertically and horizontally seam hole structure is Compartmentalization
Substantially, same fracture belt causes oil reservoir filling to differ greatly because degree of crushing is different, to follow-up reservoir prediction and reservoir model-building
It is proposed larger challenge.
Current disconnected solution oil reservoir depicting method is broadly divided into following a few classes:
(1) reservoir prediction (patent No. in plane is carried out using amplitude change rate technology:CN201010534865.8).
The technology can be conveniently used for the selection of well location in plane, but can not characterize longitudinal architectural feature;
(2) energy body technique is utilized, it is special that distribution of the beading anomalous body in longitudinal direction can be portrayed by the calculating of amplitude envelope
Sign (Bao allusion quotation, the Compartmentalization description technique study of the disconnected solution oil reservoir of Zhang Huitao Carbonate Rocks In Tahe Oilfields, Xinjiang's oil and gas,
2017).The technology is influenceed by wavelet side lobe effect so that Reservoir Body scale is bigger than normal, it is more difficult to the longitudinal subdivision of Reservoir Body
Portrayed;
(3) using the successional seismic properties of earthquake are calculated, such as relevant, ant body, the fault development spy for portraying underground
Sign (Wu Bo etc., application of the ant tracer technique in fracture-pore reservoir FRACTURE PREDICTION, fault-blcok oil-gas field, 2014).But can not
The degree of crushing of fracture is portrayed;
(4) tensor attribute is utilized, calculating the discontinuity of the underground space, (Lee is first-class, and gradient-structure tensor analysis method is three
Application in dimension seismic data channel sand prediction, gas industry, 2011), fracture fracture area is carried out in conjunction with data on well
Portray, but portraying for Reservoir Body can not be characterized.Never ipsilateral has carried out the sign of carbonate rock Reservoir Body to above method,
But portray result codomain because of distinct methods to differ greatly, can not directly carry out the integrated forecasting of Reservoir Body.
The content of the invention
In view of the above problems, it is proposed that the present invention can portray disconnected solution profile and description Reservoir Body to provide one kind
The characterizing method of the disconnected solution internal structure of the carbonate rock of locus.
The characterizing method of solution internal structure according to an aspect of the invention, there is provided a kind of carbonate rock breaks, including:
Post-stack inversion is carried out to seismic data, obtains underground Acoustic Impedance Data;
It is determined that the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development is represented, according to the Acoustic Impedance Data threshold value to described
Underground Acoustic Impedance Data carries out through-carved work processing;
Calculate the tensor attribute of seismic data;
Space smoothing processing is carried out to the tensor attribute, it is determined that the tensor attribute thresholds that disconnected solution crushes region are represented,
The result of tensor attribute space smoothing processing is handled according to the tensor attribute thresholds, obtains the profile of disconnected solution;
Using the profile of the disconnected solution as border, the tensor attribute and the underground Acoustic Impedance Data are included taking turns
In exterior feature, the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.
Further, described to carry out post-stack inversion to seismic data, obtaining underground Acoustic Impedance Data is specially:Using based on
The inversion method of model carries out post-stack inversion to seismic data, obtains underground Acoustic Impedance Data.
Further, it is described that post-stack inversion is carried out to seismic data, obtain underground Acoustic Impedance Data and further comprise:
Step 1, using Geological Mode as guidance, well-log information is combined with seismic data, establishes initial geological model;
Step 2, statistical wavelet is extracted according to the seismic data, carrying out forward modeling to the initial geological model is closed
Into record;
Step 3, the composite traces is fitted with seismic data, judges whether fitting result meets preparatory condition,
If so, then terminate iteration;If it is not, then ask for the residual error between the composite traces and the seismic data;
Step 4, using the parameter of the residual error adjustment initial geological model, step 2 is continued executing with.
Further, described be fitted composite traces with seismic data further comprises:Will using least squares method
Composite traces is fitted with seismic data.
Further, the determination represents the Acoustic Impedance Data threshold value that fracture hole type Reservoir Body is developed and further comprised:According to
Well-log information determines to represent the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development.
Further, the tensor attribute for calculating seismic data includes:
Step 1:Structure tensor is established, calculates the every bit gradient vector of 3D seismic data;
Step 2:According to the result of calculation of the every bit gradient vector of 3D seismic data, the structure tensor is calculated;
Step 3:Each composition of the structure tensor is carried out using Gaussian window smooth;
Step 4:Smooth result is carried out to each composition of the structure tensor according to Gaussian window, calculates the structure
The characteristic value and characteristic vector of tensor.
Further, it is characterised in that described to determine that representing the tensor attribute thresholds that disconnected solution crushes region further wraps
Include:The tensor attribute for being spatially smoothed processing is demarcated using drilling time curve, it is determined that representing disconnected solution crushes region
Tensor attribute thresholds.
The present invention obtains underground Acoustic Impedance Data by carrying out post-stack inversion to seismic data;It is determined that represent the storage of fracture hole type
The Acoustic Impedance Data threshold value of collective's development, hollow out carving is carried out to the underground Acoustic Impedance Data according to the Acoustic Impedance Data threshold value
Quarter is handled;Calculate the tensor attribute of seismic data;Space smoothing processing is carried out to the tensor attribute, broken it is determined that representing disconnected solution
The tensor attribute thresholds in broken region, according to the tensor attribute thresholds to the result of tensor attribute space smoothing processing at
Reason, obtain the profile of disconnected solution;Using the profile of the disconnected solution as border, by the tensor attribute and the underground wave impedance
Data are shown in profile, and the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.The present invention utilizes geological data meter
Obtained wave impedance and tensor attribute, the comprehensive characterization for the disconnected solution of technology row carbonate rock being stacked using figure, can be carved
The border of disconnected solution oil reservoir and the spatial distribution of Reservoir Body are drawn, it is comprehensive to embody the profile and inner structural features for resolving solution,
Relative to the Predicting Technique of early stage, the technology can preferably depict the scope of disconnected solution and the otherness of internal Reservoir Body.
The invention can help reservoir geology personnel more clearly epistemically under oil pool distribution and border, be follow-up well site deployment and
Reservoir model-building provides data basis.
Described above is only the general introduction of technical solution of the present invention, in order to better understand the technological means of the present invention,
And can be practiced according to the content of specification, and in order to allow above and other objects of the present invention, feature and advantage can
Become apparent, below especially exemplified by the embodiment of the present invention.
Brief description of the drawings
By reading the detailed description of hereafter preferred embodiment, it is various other the advantages of and benefit it is common for this area
Technical staff will be clear understanding.Accompanying drawing is only used for showing the purpose of preferred embodiment, and is not considered as to the present invention
Limitation.And in whole accompanying drawing, identical part is denoted by the same reference numerals.In the accompanying drawings:
Fig. 1 shows the flow chart of the characterizing method according to a kind of disconnected solution internal structure of carbonate rock of the present invention;
Fig. 2 shows according to embodiments of the present invention one seismic data;
Fig. 3 shows according to embodiments of the present invention one wave impedance section figure;
Fig. 4 shows the profile of according to embodiments of the present invention one wave impedance after through-carved work;
Fig. 5 shows according to embodiments of the present invention one tensor attribute profile;
Fig. 6 shows according to embodiments of the present invention one result carried out to tensor attribute after space smoothing;
Fig. 7 shows the profile of the according to embodiments of the present invention one disconnected solution depicted according to tensor threshold value;
Fig. 8 shows according to embodiments of the present invention one disconnected solution internal structure phenogram.
Embodiment
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
Limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure
Completely it is communicated to those skilled in the art.
Fig. 1 shows the characterizing method flow chart for solution internal structure of breaking according to a kind of carbonate rock of the present invention, this method
Comprise the following steps:
Step S101:Post-stack inversion is carried out to seismic data, obtains underground Acoustic Impedance Data, the purpose is to carbonate rock
The Reservoir Body of disconnected solution is portrayed.
Optionally, post-stack inversion is carried out to seismic data using based on the inversion method of model, obtains underground wave impedance number
According to.
Further, post-stack inversion is carried out to seismic data using based on the inversion method of model, obtains underground wave impedance
Data include:
Step 1, using Geological Mode as guidance, well-log information is combined with seismic data, establishes initial geological model;
Step 2, statistical wavelet is extracted according to seismic data, carrying out forward modeling to initial geological model obtains composite traces;
Step 3, composite traces and seismic data are fitted, judge whether fitting result meets preparatory condition, if so,
Then terminate iteration;If it is not, then ask for the residual error between composite traces and seismic data;
Step 4, the parameter of initial geological model is adjusted using residual error, continues executing with step 2.
Optionally, composite traces and seismic data are fitted specially in step 3:It will be synthesized using least squares method
Record is fitted with seismic data.Judge whether composite traces and the result of seismic data fitting meet preparatory condition, for example,
Composite traces and seismic data then meet preparatory condition through over-fitting, during both difference minimums, now terminate iteration.
Step S102:It is determined that the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development is represented, according to Acoustic Impedance Data threshold value
Through-carved work processing is carried out to underground Acoustic Impedance Data.
Well logging is using geophysical properties such as the electrochemical properties of rock stratum, conductive characteristic, acoustic characteristic, radioactivity, is surveyed
The method for measuring geophysical parameterses, analyzed by the well-log information obtained to well logging, to obtain various petroleum geologies and skill
Art data.
The real drilling well well-log information for being bored in work area and meeting fracture hole type Reservoir Body is collected, is determined to represent fracture hole type according to well-log information
The Acoustic Impedance Data threshold value of Reservoir Body development area, using Graphics Software, according to threshold value by Reservoir Body development area underground wave resistance
Anti- data through-carved work comes out.
Step S103:Calculate the tensor attribute of seismic data.Tensor attribute judges seismic event using earthquake tensor field
Trend, its core be exactly with reference to judging that fluid potential energy identical method judges seismic event occurrence, it is different so as to reach fracture hole
General knowledge is other.
Optionally, the tensor attribute for calculating seismic data comprises the following steps:
Step 1:Structure tensor is established, calculates the every bit gradient vector of 3D seismic data;
Step 2:According to the result of calculation of the every bit gradient vector of 3D seismic data, structure tensor is calculated;
Step 3:Each composition of structure tensor is carried out using Gaussian window smooth;
Step 4:Smooth result is carried out to each composition of structure tensor according to Gaussian window, calculates the spy of structure tensor
Value indicative and characteristic vector.
Wherein, smooth purpose is carried out to each composition of structure tensor using Gaussian window to be to use smooth gradient structure
Tensor characterizes the textural characteristics in certain area and suppresses the mutation of the structure tensor caused by noise;In the present invention, meter
The characteristic value of the structure tensor drawn is tensor attribute, such as in actual application, with second characteristic value pair
Disconnected solution is portrayed.
Step S104:Space smoothing processing is carried out to tensor attribute, it is determined that representing the tensor attribute that disconnected solution crushes region
Threshold value, the result of tensor attribute space smoothing processing is handled according to tensor attribute thresholds, obtain the profile of disconnected solution, profit
Disconnected solution is portrayed with the tensor attribute after processing is spatially smoothed, portrays the form of result and actual disconnected solution more
It coincide.
Optionally, the tensor attribute for being spatially smoothed processing is demarcated using drilling time curve, it is determined that representing disconnected molten
Body crushes the tensor attribute thresholds in region.
Fracture horizontal broken scope, the positive sidetracking for crossing fracture are demarcated using the drilling time curve of sidetracked hole in work area and horizontal well
Horizontal well is mostly to be vented absorption well, typically no log, the changing rule of drilling time curve can be utilized to judge the core of fracture
Portion's width and crushed zone scope.The drilling time curve of basement rock normally behaves as stable changing rule, when entering broken scope, bores
When curve change, overall trend shows as linear change feature, when broken corrosion band shows as low brill, and along with emptying
Leakage phenomena.Space smoothing is carried out to tensor attribute, entered using when boring with having carried out the tensor attribute after space smoothing is handled
Rower is determined, it is determined that the tensor door threshold value in broken region, obtains the profile of disconnected solution.
Step S105:Using the profile for the solution that breaks as border, tensor attribute and underground Acoustic Impedance Data are included in profile
Interior, the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.
The present invention obtains underground Acoustic Impedance Data by carrying out post-stack inversion to seismic data;It is determined that represent the storage of fracture hole type
The Acoustic Impedance Data threshold value of collective's development, hollow out carving is carried out to the underground Acoustic Impedance Data according to the Acoustic Impedance Data threshold value
Quarter is handled;Calculate the tensor attribute of seismic data;Space smoothing processing is carried out to the tensor attribute, broken it is determined that representing disconnected solution
The tensor attribute thresholds in broken region, according to the tensor attribute thresholds to the result of tensor attribute space smoothing processing at
Reason, obtain the profile of disconnected solution;Using the profile of the disconnected solution as border, by the tensor attribute and the underground wave impedance
Data are shown in profile, and the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.The present invention utilizes geological data meter
Obtained wave impedance and tensor attribute, the comprehensive characterization for the disconnected solution of technology row carbonate rock being stacked using figure, can be carved
The border of disconnected solution oil reservoir and the spatial distribution of Reservoir Body are drawn, it is comprehensive to embody the profile and inner structural features for resolving solution,
Relative to the Predicting Technique of early stage, the technology can preferably depict the scope of disconnected solution and the otherness of internal Reservoir Body.
The invention can help reservoir geology personnel more clearly epistemically under oil pool distribution and border, be follow-up well site deployment and
Reservoir model-building provides data basis.
Fig. 2 shows according to embodiments of the present invention one seismic data.
Fig. 3 shows according to embodiments of the present invention one wave impedance section figure.
Can be seen that wave impedance inversion with reference to Fig. 2 and Fig. 3 can eliminate wavelet side lobe effect, and its longitudinal frame is higher than
Seismic data, the spatial shape and its Distribution Characteristics of fracture hole can be preferably portrayed, improve reservoir prediction precision.
Fig. 4 shows the profile of according to embodiments of the present invention one wave impedance after through-carved work.Utilize real drilling well
Well-log information is demarcated to Acoustic Impedance Data, it is determined that the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development is represented, according to threshold
Value carries out through-carved work to Acoustic Impedance Data.In the present embodiment, threshold value is chosen for 14500 after being demarcated according to actual conditions
~16100m/s*gcc.
Fig. 5 shows according to embodiments of the present invention one tensor attribute profile.The main body of otherness inside disconnected solution
On present degree of crushing, after breaking to form, later stage dissolution extent is strong, then can form larger Caves Reservoir Body, corrosion journey
Spend weak region and form fracture hole type, slit formation Reservoir Body, with reference to Fig. 2 and Fig. 5 can be seen that tensor attribute can not only show it is disconnected
The border of solution, additionally it is possible to portray disconnected solution internal fracture degree.
Fig. 6 shows according to embodiments of the present invention one result carried out to tensor attribute after space smoothing.
Fig. 7 shows the profile of the according to embodiments of the present invention one disconnected solution depicted according to tensor threshold value.
It is preferred that space smoothing processing is carried out to tensor attribute, using drilling time curve to being spatially smoothed the tensor category of processing
Property demarcated, it is determined that disconnected solution crushes region tensor attribute thresholds, using the threshold value to the tensor after space smoothing is handled
Attribute is handled, and obtains the profile of disconnected solution.Threshold value in the present embodiment is defined as 2~150.
Fig. 8 shows according to embodiments of the present invention one disconnected solution internal structure phenogram.
As can be seen from Figure 8, in disconnected solution border, solution internal fracture degree of breaking change and High-quality Reservoir distribution are equal
Clearly portrayed, and there is good Matching Relationship in the strong high region of corrosion degree of crushing with High-quality Reservoir development region, meets
Geological knowledge rule.From real drilling well reservoir development situation analysis, prediction result is coincide all right with real brill.
In the specification that this place provides, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention
Example can be put into practice in the case of these no details.In some instances, known method, structure is not been shown in detail
And technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the disclosure and help to understand one or more of each inventive aspect,
Above in the description to the exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention:I.e. required guarantor
The application claims of shield features more more than the feature being expressly recited in each claim.It is more precisely, such as following
Claims reflect as, inventive aspect is all features less than single embodiment disclosed above.Therefore,
Thus the claims for following embodiment are expressly incorporated in the embodiment, wherein each claim is in itself
Separate embodiments all as the present invention.
In the specification that this place provides, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention
Example can be put into practice in the case of these no details.In some instances, known method, structure is not been shown in detail
And technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the disclosure and help to understand one or more of each inventive aspect,
Above in the description to the exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention:I.e. required guarantor
The application claims of shield features more more than the feature being expressly recited in each claim.It is more precisely, such as following
Claims reflect as, inventive aspect is all features less than single embodiment disclosed above.Therefore,
Thus the claims for following embodiment are expressly incorporated in the embodiment, wherein each claim is in itself
Separate embodiments all as the present invention.
Those skilled in the art, which are appreciated that, to be carried out adaptively to the module in the equipment in embodiment
Change and they are arranged in one or more equipment different from the embodiment.Can be the module or list in embodiment
Member or component be combined into a module or unit or component, and can be divided into addition multiple submodule or subelement or
Sub-component.In addition at least some in such feature and/or process or unit exclude each other, it can use any
Combination is disclosed to all features disclosed in this specification (including adjoint claim, summary and accompanying drawing) and so to appoint
Where all processes or unit of method or equipment are combined.Unless expressly stated otherwise, this specification (including adjoint power
Profit requires, summary and accompanying drawing) disclosed in each feature can be by providing the alternative features of identical, equivalent or similar purpose come generation
Replace.
In addition, it will be appreciated by those of skill in the art that although some embodiments described herein include other embodiments
In included some features rather than further feature, but the combination of the feature of different embodiments means in of the invention
Within the scope of and form different embodiments.For example, in the following claims, embodiment claimed is appointed
One of meaning mode can use in any combination.
It should be noted that the present invention will be described rather than limits the invention for above-described embodiment, and ability
Field technique personnel can design alternative embodiment without departing from the scope of the appended claims.In the claims,
Any reference symbol between bracket should not be configured to limitations on claims.Word "comprising" does not exclude the presence of not
Element or step listed in the claims.Word "a" or "an" before element does not exclude the presence of multiple such
Element.The present invention can be by means of including the hardware of some different elements and being come by means of properly programmed computer real
It is existing.In if the unit claim of equipment for drying is listed, several in these devices can be by same hardware branch
To embody.The use of word first, second, and third does not indicate that any order.These words can be explained and run after fame
Claim.
Claims (7)
- The characterizing method of solution internal structure 1. a kind of carbonate rock breaks, it is characterised in that the described method comprises the following steps:Post-stack inversion is carried out to seismic data, obtains underground Acoustic Impedance Data;It is determined that the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development is represented, according to the Acoustic Impedance Data threshold value to the underground Acoustic Impedance Data carries out through-carved work processing;Calculate the tensor attribute of seismic data;Space smoothing processing is carried out to the tensor attribute, it is determined that the tensor attribute thresholds that disconnected solution crushes region are represented, according to The tensor attribute thresholds are handled the result of tensor attribute space smoothing processing, obtain the profile of disconnected solution;Using the profile of the disconnected solution as border, the tensor attribute and the underground Acoustic Impedance Data are included in profile Interior, the exterior contour of the disconnected solution of realization is portrayed and internal structure characterizes.
- 2. according to the method for claim 1, it is characterised in that it is described that post-stack inversion is carried out to seismic data, obtain underground Acoustic Impedance Data is specially:Post-stack inversion is carried out to seismic data using based on the inversion method of model, obtains underground wave impedance Data.
- 3. according to the method for claim 2, it is characterised in that it is described that post-stack inversion is carried out to seismic data, obtain underground Acoustic Impedance Data further comprises:Step 1, using Geological Mode as guidance, well-log information is combined with seismic data, establishes initial geological model;Step 2, statistical wavelet is extracted according to the seismic data, carrying out forward modeling to the initial geological model obtains synthesis note Record;Step 3, the composite traces is fitted with seismic data, judges whether fitting result meets preparatory condition, if so, Then terminate iteration;If it is not, then ask for the residual error between the composite traces and the seismic data;Step 4, using the parameter of the residual error adjustment initial geological model, step 2 is continued executing with.
- 4. according to the method for claim 3, it is characterised in that described that composite traces and seismic data are fitted into one Step includes:Composite traces and seismic data are fitted using least squares method.
- 5. according to the method for claim 1, it is characterised in that the determination represents the wave impedance of fracture hole type Reservoir Body development Data threshold further comprises:Determined to represent the Acoustic Impedance Data threshold value of fracture hole type Reservoir Body development according to well-log information.
- 6. according to the method for claim 1, it is characterised in that the tensor attribute for calculating seismic data includes:Step 1:Structure tensor is established, calculates the every bit gradient vector of 3D seismic data;Step 2:According to the result of calculation of the every bit gradient vector of 3D seismic data, the structure tensor is calculated;Step 3:Each composition of the structure tensor is carried out using Gaussian window smooth;Step 4:Smooth result is carried out to each composition of the structure tensor according to Gaussian window, calculates the structure tensor Characteristic value and characteristic vector.
- 7. according to the method for claim 1, it is characterised in that described to determine to represent the tensor attribute that disconnected solution crushes region Threshold value further comprises:The tensor attribute for being spatially smoothed processing is demarcated using drilling time curve, it is determined that representing disconnected molten Body crushes the tensor attribute thresholds in region.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710595771.3A CN107390264B (en) | 2017-07-20 | 2017-07-20 | The characterizing method of the disconnected solution internal structure of carbonate rock |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710595771.3A CN107390264B (en) | 2017-07-20 | 2017-07-20 | The characterizing method of the disconnected solution internal structure of carbonate rock |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107390264A true CN107390264A (en) | 2017-11-24 |
CN107390264B CN107390264B (en) | 2019-06-07 |
Family
ID=60336359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710595771.3A Active CN107390264B (en) | 2017-07-20 | 2017-07-20 | The characterizing method of the disconnected solution internal structure of carbonate rock |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107390264B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007345A (en) * | 2019-04-08 | 2019-07-12 | 中国石油化工股份有限公司 | A kind of disconnected solution oil-gas reservoir cave type reservoir finimeter calculation method and device |
CN110308487A (en) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | A kind of disconnected solution type oil reservoir quantitatively characterizing method |
CN110471107A (en) * | 2018-05-10 | 2019-11-19 | 中国石油天然气股份有限公司 | The depicting method and device of the disconnected solution of carbonate rock |
CN110488385A (en) * | 2019-09-20 | 2019-11-22 | 西南石油大学 | The production method of the disconnected solution model of three-dimensional |
CN110794476A (en) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | Inversion method based on phase control of broken solution phase |
CN111435173A (en) * | 2019-01-15 | 2020-07-21 | 中国石油天然气集团有限公司 | Rock reservoir post-stack seismic data denoising method and device |
CN111830558A (en) * | 2019-04-17 | 2020-10-27 | 中国石油化工股份有限公司 | Fracture zone engraving method |
CN112114357A (en) * | 2019-06-19 | 2020-12-22 | 中国石油天然气股份有限公司 | Method and device for predicting connectivity of fracture-cavity reservoir |
CN112649862A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and device for identifying broken solution based on stratum structure information separation |
CN112731527A (en) * | 2019-10-14 | 2021-04-30 | 中国石油化工股份有限公司 | Multi-attribute research-based method and device for enhancing characteristics of broken solution |
CN113514893A (en) * | 2020-04-10 | 2021-10-19 | 中国石油天然气股份有限公司 | Method and device for describing internal structure of broken zone of sliding fracture of carbonate rock |
RU2815465C1 (en) * | 2022-11-22 | 2024-03-18 | Саутвест Петролеум Юниверсити | System for semi-deterministic modeling of fracturing based on random matrix |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070280048A1 (en) * | 2006-06-06 | 2007-12-06 | Baker Hughes Incorporated | P-wave anisotropy evaluation by measuring acoustic impedance of the rock by beam-steering from within the borehole at different angles |
CN106896405A (en) * | 2017-02-28 | 2017-06-27 | 中国石油化工股份有限公司 | A kind of sand-conglomerate body spatial Forecasting Methodology and device based on chaos attribute |
-
2017
- 2017-07-20 CN CN201710595771.3A patent/CN107390264B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070280048A1 (en) * | 2006-06-06 | 2007-12-06 | Baker Hughes Incorporated | P-wave anisotropy evaluation by measuring acoustic impedance of the rock by beam-steering from within the borehole at different angles |
CN106896405A (en) * | 2017-02-28 | 2017-06-27 | 中国石油化工股份有限公司 | A kind of sand-conglomerate body spatial Forecasting Methodology and device based on chaos attribute |
Non-Patent Citations (1)
Title |
---|
梁生贤 等: "大地电磁阻抗张量模的拟地震解释法研究及应用", 《地球物理学进展》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110308487A (en) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | A kind of disconnected solution type oil reservoir quantitatively characterizing method |
CN110471107A (en) * | 2018-05-10 | 2019-11-19 | 中国石油天然气股份有限公司 | The depicting method and device of the disconnected solution of carbonate rock |
CN110794476B (en) * | 2018-08-01 | 2021-10-19 | 中国石油化工股份有限公司 | Inversion method based on phase control of broken solution phase |
CN110794476A (en) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | Inversion method based on phase control of broken solution phase |
CN111435173A (en) * | 2019-01-15 | 2020-07-21 | 中国石油天然气集团有限公司 | Rock reservoir post-stack seismic data denoising method and device |
CN111435173B (en) * | 2019-01-15 | 2023-04-25 | 中国石油天然气集团有限公司 | Method and device for denoising seismic data after lithologic oil reservoir stack |
CN110007345A (en) * | 2019-04-08 | 2019-07-12 | 中国石油化工股份有限公司 | A kind of disconnected solution oil-gas reservoir cave type reservoir finimeter calculation method and device |
CN111830558A (en) * | 2019-04-17 | 2020-10-27 | 中国石油化工股份有限公司 | Fracture zone engraving method |
CN111830558B (en) * | 2019-04-17 | 2022-07-22 | 中国石油化工股份有限公司 | Fracture zone engraving method |
CN112114357A (en) * | 2019-06-19 | 2020-12-22 | 中国石油天然气股份有限公司 | Method and device for predicting connectivity of fracture-cavity reservoir |
CN110488385A (en) * | 2019-09-20 | 2019-11-22 | 西南石油大学 | The production method of the disconnected solution model of three-dimensional |
CN112649862A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Method and device for identifying broken solution based on stratum structure information separation |
CN112731527A (en) * | 2019-10-14 | 2021-04-30 | 中国石油化工股份有限公司 | Multi-attribute research-based method and device for enhancing characteristics of broken solution |
CN113514893A (en) * | 2020-04-10 | 2021-10-19 | 中国石油天然气股份有限公司 | Method and device for describing internal structure of broken zone of sliding fracture of carbonate rock |
RU2815465C1 (en) * | 2022-11-22 | 2024-03-18 | Саутвест Петролеум Юниверсити | System for semi-deterministic modeling of fracturing based on random matrix |
Also Published As
Publication number | Publication date |
---|---|
CN107390264B (en) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107390264B (en) | The characterizing method of the disconnected solution internal structure of carbonate rock | |
CN107367757A (en) | The depicting method of the disconnected solution of carbonate rock | |
US10330825B2 (en) | Complex fracture network mapping | |
CN107515423B (en) | A kind of new 3 D geological model construction method of high-dip structure | |
CN104977617A (en) | Reservoir fracture identification method and imaging logging reservoir fracture identification method | |
CN104678436A (en) | Method for predicting thickness of compressional fault zone on coverage area | |
CN110658556A (en) | Seismic technology combination method for identifying and evaluating broken zone of sliding fracture of carbonate rock | |
CN106089191B (en) | Method for identifying pressure fracture zone structure well logging | |
CN111722277A (en) | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration | |
CN107798189A (en) | One kind is based on the geophysical precisely dynamic anti-burst method of Transparent space | |
CN110532507B (en) | Method for improving drilling rate of compact oil reservoir of horizontal well | |
Baoping et al. | A drilling technology guided by well-seismic information integration | |
Delorme et al. | A methodology to characterize fractured reservoirs constrained by statistical geological analysis and production: a real field case study | |
Ahmed et al. | Geological Model for Mauddud Reservoir Khabaz Oil Field, Kirkuk, Northern Iraq | |
CN104251135B (en) | Highly-deviated well space in-place method | |
CN110007345A (en) | A kind of disconnected solution oil-gas reservoir cave type reservoir finimeter calculation method and device | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
CN114740547A (en) | Back-flushing sliding fracture fold grading and partitioning analysis and three-dimensional characterization method | |
Al-Anazi et al. | Automatic fracture density update using smart well data and artificial neural networks | |
CN112394416A (en) | Heterogeneous fracture control reservoir prediction method and device | |
CN111963148A (en) | Method for determining pressure coefficient and drilling mud density of carbonate fracture-cave type oil and gas reservoir | |
US12031900B2 (en) | Corrosion analysis using magnetic flux leakage measurements | |
Khaksar et al. | Sanding evaluation and real-time sand management decisions in high-rate gas wells, Ichthys Field, Browse Basin Australia | |
US20240230513A9 (en) | Corrosion analysis using magnetic flux leakage measurements | |
Chinellato et al. | Logging-While-Drilling Solutions for Sustainable Underground Gas Storage Development |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |