CN109143399A - A method of identification carbonate rock sequence interface - Google Patents

A method of identification carbonate rock sequence interface Download PDF

Info

Publication number
CN109143399A
CN109143399A CN201710507954.5A CN201710507954A CN109143399A CN 109143399 A CN109143399 A CN 109143399A CN 201710507954 A CN201710507954 A CN 201710507954A CN 109143399 A CN109143399 A CN 109143399A
Authority
CN
China
Prior art keywords
log
porosity
sequence interface
interface
frequency sequence
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
Application number
CN201710507954.5A
Other languages
Chinese (zh)
Other versions
CN109143399B (en
Inventor
李宏涛
肖开华
胡向阳
史云清
李�浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710507954.5A priority Critical patent/CN109143399B/en
Publication of CN109143399A publication Critical patent/CN109143399A/en
Application granted granted Critical
Publication of CN109143399B publication Critical patent/CN109143399B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

This application provides a kind of methods for identifying carbonate rock sequence interface, this method comprises: the core sample using carbonate rock demarcates at least one Logging Curves;Lithology sensitivity log is determined from least one Logging Curves;It is combined according to the Resistivity log in lithology sensitivity log and determines tri-porosity logging curve combination;Identification low frequency sequence interface is combined according to lithology sensitivity log, tri-porosity logging curve combination and Resistivity log;High frequency sequence interface is identified according to low frequency sequence interface.According to this method, solve the problems, such as that coring is few and representative difference is difficult to accurately identify sequence interface, and this method strong operability, it is at low cost, it haves laid a good foundation for the sequence deposition understanding and overall merit of carbonate rock.

Description

A method of identification carbonate rock sequence interface
Technical field
The present invention relates to oil-gas explorations, development field, and identify more particularly, to a kind of carbonate rock sequence interface Method.
Background technique
The identification of sequence interface is basis and the key component of sequence stratigraphy study, the conventional method of sequence interface identification Mainly comprehensive descision is carried out using many data such as rock core, earthquake, well logging and extinct plants and animal.Well-log information is in Sequence Stratigraphic Units It has very important significance in the identification of interface, especially studies area in the covering of sedimentary basin, seismic data resolution is limited, takes Heart data is limited and discontinuous, without outcrop data, it is difficult to determine accurate SEQUENCE STRATIGRAPHIC Boundary.At this point, in conjunction with limited Rock core information carry out comprehensive analysis, with well-log information carry out SEQUENCE STRATIGRAPHIC DIVISION, comparison be particularly important.Based on well logging The method that data carries out Sequence Stratigraphic Units division can be divided into qualitative and two kinds of quantitative judge.Quantitative identification method frequency division sometimes Analysis technology, INPEFA technology, multi-scale Wavelet Analysis technology, Empirical mode decomposition, log fractals and well logging are more Scale data fusion method etc..But Quantitative research method geological Significance is indefinite, live practical application.Qualitative method has certainly Right current potential and natural potential image method, interval transit time and resistivity curve overlay method, sonic time difference, natural potential and apparent resistance Rate combined method and accumulative stratigraphic dip well logging method etc..However, above method is mainly used in clastic deposited sediments stratum, and exist It is certain to divide multi-solution.
Reservoir division comparison fineness determines layer description fine degree, but the fine Division and contrast of carbonate formation compared with Clastic rock is more complicated, because the size and circularity of carbonate rock particle cannot significantly reflect the process and energy of hydrodynamism Amount, diagenetic complexity often enhance reservoir heterogeneity, weaken connate deposit cyclicity.Studies have shown that utilizing well logging Identify that carbonate rock sequence interface is one of the effective way that solves the above problems.
Although forefathers have carried out certain identification application to sequence interface according to well-log information, generally directed to be low frequency The interface of third-band sequence and the above rank.Using Logging Curves to high frequency sequence interface identify and compare, document report compared with It is few, at present still in the trial stage, systematic research and application method are not yet summed up, in especially classical sequence stratigraphy The high frequency sequence (Pyatyi sequence or six grades of sequences) of subsequence rank.Peng Jun etc. is in " the wellblock TK625 hawk mountain group carbonate reservoir Fine correlation research " and trip Li Wei etc. " application of the high frequency sequence research in carbonate reservoir compares --- with Tahe For the reservoir of the oil field wellblock the TK625 top Ying Shanzu " in a text, selection is bent with resistivity, natural gamma, sound wave and density log Line has divided the wellblock the system in Tahe Oilfield TK625 top Ying Shanzu reservoir fourth level sequence, but description is excessively general, not to fourth level sequence The log at interface responds and variation characteristic carries out detailed analysis.
Therefore, it is necessary to a kind of improved methods for accurately identifying carbonate rock sequence interface.
Summary of the invention
For above-mentioned the problems of the prior art, present applicant proposes a kind of carbonate rock sequence interface identification methods, lead to It crosses and demarcates log using rock core, and then identify and determine stratigraphic sequence interface, solve carbonate rock complex area brill Well coring is few, it is difficult to the problem of accurately identifying high frequency sequence interface.
On the one hand, a kind of method for identifying carbonate rock sequence interface is provided, this method comprises: step 1, utilizing carbonic acid The core sample of rock salt demarcates at least one Logging Curves;Step 2 determines lithology from least one Logging Curves Sensitive log, wherein lithology sensitivity log includes Resistivity log combination;Step 3, according to resistivity logging Curve combination determines tri-porosity logging curve, wherein tri-porosity logging curve combination include acoustic porosity log, Density porosity log and neutron porosity log curve;Resistivity log combination includes that deep lateral resistivity is logged well Curve and shallow lateral Resistivity log;Step 4, according to lithology sensitivity log, tri-porosity logging curve combination and Resistivity log combination identification low frequency sequence interface, low frequency sequence interface are fourth level sequence or the boundary of its above rank sequence Face;Step 5 identifies high frequency sequence interface according to low frequency sequence interface, and wherein high frequency sequence interface is Pyatyi sequence interface or six Grade sequence interface.According to the method for this aspect, solve the problems, such as that coring is few and representative difference is difficult to accurately identify sequence interface, And this method operability is asked, it is at low cost, meet sequence stratigraphy and well logging learn principle, be carbonate rock sequence deposition understanding and Overall merit is had laid a good foundation.
According to a kind of possible implementation of this aspect, step 1 includes: step 11, the target zone using carbonate rock Core sample makes core sample thin slice;Step 12, observation core sample thin slice simultaneously demarcate at least one Logging Curves.Root According to the implementation, the principal lithologic and depositional environment of target zone carbonate rock can determine.
According to a kind of possible implementation of this aspect, step 3 includes: step 31, according to Resistivity log group It closes and determines compact limestone section;Step 32, the interval transit time skeleton value and density skeleton value for counting compact limestone section, and utilize sound wave Time difference skeleton value and density skeleton value calculate separately acoustic porosity and density porosity;Step 33, the sound wave hole by calculating Degree and the density porosity calculated and neutron porosity are placed in same log road, to form tri-porosity logging curve group It closes.According to a kind of possible implementation of this aspect, step 3 further include: step 34, using actual measurement core porosity to three holes The combination of porosity log is corrected.
According to a kind of possible implementation of this aspect, step 4 is specifically included: step 41, by tri-porosity logging song Line combination is overlapped at compact limestone section;Step 42, according to lithology sensitivity log, tri-porosity logging curve combination and electricity The combination of resistance rate log, bond area drilling stratum, lithologic character data identify low frequency sequence interface.According to this Possible implementation can be identified low using the compact limestone section as index bed by the mutation of feature on log Frequency sequence interface.
According to a kind of possible implementation of this aspect, step 5 is specifically included: step 51, observation target zone core sample Product;Lithology above and below step 52, the variation of lithological trend for summarizing high frequency sedimentary cycle and different high frequency sedimentary cycle contact surfaces Mutation differences identify high frequency sequence interface.
According to a kind of possible implementation of this aspect, step 5 further include: step 53, according to high frequency sequence interface, weight Multiple step 2 and step 3 are to identify full well section high frequency sequence interface.It, can be to the full well of individual well by the possible implementation Section high frequency sequence interface is identified that solving well section can not identify asking for sequence interface entirely by core sample in the prior art Topic.
According to a kind of possible implementation of this aspect, this method further include: step 7, to low frequency sequence interface and height Frequency sequence interface is verified.It, can be to identified low frequency sequence interface and high frequency layer by the possible implementation Sequence interface is identified, it is ensured that the accuracy of recognition result.
Above-mentioned technical characteristic may be combined in various suitable ways or be substituted by equivalent technical characteristic, as long as can reach To the purpose of the present invention.
Detailed description of the invention
The invention will be described in more detail below based on embodiments and refering to the accompanying drawings.Wherein:
Fig. 1 is the process block diagram according to the method at the identification carbonate rock sequence interface of the embodiment of the present invention.
Fig. 2 is the lithology sensitivity logging curve according to the embodiment of the present invention.
Fig. 3 is the logging curve demarcated according to the high frequency sequence interface of identification according to the embodiment of the present invention.
Fig. 4 is the logging curve demarcated according to the high frequency sequence interface of identification according to another embodiment of the present invention.
Fig. 5 is the high frequency sequence division result comparison diagram according to another embodiment of the present invention.
In the accompanying drawings, identical component uses identical appended drawing reference.Attached drawing is not according to actual ratio.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings.
It should be understood that in the present invention, low frequency sequence refers to first order sequence, second-order sequence, third-band sequence and fourth level sequence, high frequency Sequence refers to Pyatyi sequence and six grades of sequences.
Fig. 1 is the schematic flow chart of the method 100 at identification carbonate rock sequence interface provided by the invention.Such as Fig. 1 institute Show, this method 100 includes:
S110 demarcates at least one Logging Curves using the core sample of carbonate rock.
S120 determines lithology sensitivity log from least one Logging Curves, wherein the well logging of lithology sensitivity is bent Line includes Resistivity log combination, and Resistivity log combination is lateral including deep lateral resistivity log and shallowly Resistivity log.
S130 is combined according to Resistivity log and is determined tri-porosity logging curve combination, wherein tri-porosity logging Curve combination includes acoustic porosity log, density porosity log and neutron porosity log curve.
S140 is combined according to lithology sensitivity log, tri-porosity logging curve combination and Resistivity log and is known Other low frequency sequence interface, low frequency sequence interface are the interface of fourth level sequence and the above rank sequence.
S150 identifies high frequency sequence interface according to low frequency sequence interface, and wherein high frequency sequence interface is Pyatyi sequence interface Or six grades of sequence interfaces.
Specifically, the present invention by observing interval of interest rock core in detail, at least one Logging Curves of calibration, It is preferred that the sensitivity curve of reflection lithology, summarizes the log response characteristics at lithology and high frequency sequence interface, recalculates three Porosity logging curve, and permutation and combination are high to full well section under the constraint and deposition understanding guidance of low frequency sequence interface identification Frequency sequence interface is identified, Stratigraphic Division And Correlation fineness is improved, and solves the carbonate rock area of complex lithology, drilling well Coring is few, it is difficult to the problem of accurate full well section identifies high frequency sequence interface.
In S110, firstly, observation target zone core sample in detail, and sample, after grinding production, under the microscope in detail Thin observation, identification, determine target zone carbonate rock principal lithologic and substantially depositional environment.Such as the western carbonate gas reservoirs in river are ground Study carefully area's interval of interest (thickness about 130-140m), coring section rock core, thin slice lithologic character are observed in detail, determine lithology.With For W1 well and X1 well, W1 well and X1 well are relatively complete and continuous in the coring of interval of interest, therefore are seen by detailed rock core It examines with thin slice lens-belowed identifying the results show that target zone is mainly the cloud rock and limestone of opposite low energy environment, Yun Yan is mainly algae (sand) considers cloud rock to be worth doing, algal layer folds cloud rock, algae bonds cloud rock, powder crystal cloud rock, and limestone is based on algae calcarenite, micrite, in mesh Layer bottom develop a small amount of cream matter cloud rock or cloud matter cream rock etc.;Outside biological algae removal class, other biological is further shown with respect to agensis Show the medium-low energy, opposite highrank fuel salinity water body environment that primary sedimentary environment is opposite occlusion.
Then, qualification result is observed according to the lithology of rock core, chip sample, Logging Curves is combined, in detail At least one Logging Curves are demarcated, determine the response characteristic of the various logs of different lithology.In S120, from Found in proven at least one Logging Curves the log and lithology sensitivity log sensitive to lithology with The variation tendency of variation of lithological.It is illustrated in figure 2 the lithology sensitivity log of combination, it can be found that passing through lithology from Fig. 2 Qualification result demarcates log, and natural gamma (GR) value of most (cloud matter) cream rocks is lower than 30API, and the GR value of limestone class is slightly higher In 30API~40API, the GR value highest of cloud rock class, typically larger than 35API, generally between 35~50API, limestone and dolomite GR value there is certain overlapping, cream rock, limestone, Yun Yan GR value general morphologictrend to gradually increase.Interval transit time (AC), Density (DEN), neutron porosity (CNL) can distinguish to a certain extent cloud rock, limestone (such as limestone usually relatively low AC (< 48us/ft)、DEN(2.75g/m3) and CNL (< 4%)), but DEN, AC and CNL are expanding vulnerable to drilling well, crack and corrosion because Element is affected.And deep lateral resistivity (LLD) curve, corrosion expanding by drilling well etc. is influenced smaller response to lithology More sensitive, compact limestone highest, up to 50000 ohm or more, cream rock takes second place LLD, 10000 ohm or more reachable, pure white cloud Rock is minimum, and usually less than 2000 ohm, and transition lithology such as cloud matter limestone, grey matter cloud rock, cream matter cloud rock, generally between 50000 Between~2000 ohm.Target zone lithology is mainly limestone and Yun Yan, and log LLD is with dolomite content increase in rock It is obvious that resistance gradually decreases feature.Therefore, the Lithology Discrimination of conventional logging is based on LLD, supplemented by GR, AC, DEN and CNL.Cause This, chooses GR, hole diameter (CAL), LLD, shallow lateral seven resistivity (LLS), DEN, AC and CNL Logging Curves as lithology Sensitivity curve, and by CAL, GR curve as together, by DEN, AC and CNL curve and resistivity curve respectively as together with.
In S130, on the basis of lithologic analysis and log identification lithology, combined according to Resistivity log Compact limestone section is chosen, the interval transit time average value and density average of this section of respective depth point are counted, as compact limestone Sound wave and density skeleton value.Using the porosity model of sound wave and density, acoustic porosity and density hole are calculated separately Degree, and it is placed in same log road with neutron porosity, form tri-porosity logging curve combination.Utilize actual measurement rock core hole Porosity is corrected three porosity curve, is then overlapped three porosity curve at compact limestone section.As Fig. 2 W1 well exists At 6177.4m~6181.7m, resistance curve is high resistant, and up to 200000 ohm or more, core sample porosity value is between 0.26% ~1.47%, average 0.56%, it is typical compact limestone section.Density skeleton average value harmony wave skeleton average value is respectively 2.66g/cm3With 47 μ s/ft, by calculating density porosity and acoustic porosity, (process is the prior art, is not gone to live in the household of one's in-laws on getting married herein State), and it is placed in same log road with neutron porosity, porosity logging curve combination is formed, actual measurement rock core hole is utilized Porosity is corrected three porosity curve, then substantially overlaps three porosity curve at the interval, by the depth point Index bed of the compact limestone section as subsequent sequence dividing and comparison.
The base of combination is identified and rearranged with S130 porosity logging curve in S120 lithology sensitivity features of logging curve On plinth, in S140, bond area drilling stratum, lithologic character data, the lithology reflected using log Characteristics of Mutation Mutation has the low frequency sequences interfaces such as the karst plane of unconformity (or lithology conversion surface) of drilling well to research area's target zone, and most The system tracts interfaces such as big marine flooding surface carry out well logging recognition.As existing research shows that river four sections of totality of western thunderbolt ways can divide For a third-band sequence (Lv Yuzhen etc., 2013).Further study show that four sections of thunder may further be divided into three level Four layers Sequence, purpose of horizon is named as four sections of thunder upper sub- sections in the fourth level sequence at four sections of tops of thunder, thus by target zone, as shown in Fig. 2, W1 Well fourth level sequence top interface (SSB1) is Regional karst erosion surface, is to deposit than more typical karst palaeo-weathering crust under interface, mud Matter content increases, and the reference lamina easily identified is formd on log, shows as GR high level, LLD is reduced, on interface The low GR of thick-layer shape compact limestone, the electrical measurement response characteristic difference of high LLD it is obvious.Four sections of thunder upper sub- section bottom interface (i.e. level Four layers Sequence top interface (SSB2)) be lithology interface, i.e. the interface of dolomite and cream rock, interface upper and lower well log curve GR, LLD and AC, DEN, CNL are also in more apparent Characteristics of Mutation.Since four sections of Western Sichuan thunder upper sub- sections are whole inside tableland, lack tableland Edge and continental shelf slope, thus the agensis of low level domain, fourth level sequence are mainly high-order by secondary marine invasion system tract (STST) and secondary System tract (SHST) is constituted, and secondary system domain interface is the upper bound for the lithology conversion surface that stratiform dolomite is mutated to thick-layer shape limestone Face is deposited by the limestone that quick marine invasion forms maximum marine flooding surface (SMFS), it is phase that log, which shows as thick-layer limestone section GR, To high level, for maximum marine flooding surface argillaceous sediment it is relatively more caused by.On porosity logging curve, sequence interface or system with Interface leads to three porosity curve in coincidence or separates since upper and lower lithology or physical property change.
S150 summarizes the variation of lithological of the high frequency sedimentary cycle to shoal upwards in the low frequency sequence interface of S140 identification Lithology mutation differences above and below trend, and different high frequency sedimentary cycle contact surfaces, identify high frequency sequence interface.Finally utilize The high frequency sequence interface identified in S150 by rock core, further demarcates lithology sensitivity log and tri-porosity logging is bent Line summarizes the logging response character at high frequency sequence interface, and then completes the identification of full well section high frequency sequence interface.
Such as under the constraint of identified low frequency sequence interface, rock core is further observed in detail, and high frequency sequence interface is carried out Identification, local area Pyatyi high frequency sequence interface are generally limestone/Yun Yan lithology abrupt interface, well log curve respond Characteristics of Mutation compared with Obviously (such as Resistivity log and tri-porosity logging curve).Six grades of sequence interfaces remove above-mentioned lithology surface of discontinuity, also deposit In lithofacies conversion surface, such as the boundary at the interface or bird's-eye structure cloud rock and lamellar cloud rock of Dark grey crystallite cloud rock and blocky cloud rock Face, the feature with the subsequence boundary that sedimentary water body on interface is deepened suddenly.The above lithologic interface or lithofacies interface is detailed Thin calibration log, high frequency sequence interface show as sporting low-resistance by high resistant up and down, and porosity logging curve is often by weight The low-porosity of conjunction sports separated high porosity;In each high frequency sequence, log LLD in rock dolomite contain Amount increases resistance and gradually decreases, and porosity logging curve is totally in that (density porosity is vulnerable to well radius variations shadow for variation from low to high Ring, as a result for referring to), the high frequency sequence of the deposition that obviously shoals upwards is thus constituted, as shown in Figure 3.
Preferably, after identifying low frequency sequence interface and high frequency sequence interface according to the method described above, it is also necessary to low to this Frequency sequence interface and the accuracy at high frequency sequence interface are verified.Specifically, mainly according to the high frequency of other well core observations Sequence interface and lithology sensitivity log and porosity logging curve comparison verify whether reliable.For example, can use X1 well To the above method verified as a result, as shown in figure 4, the high frequency sequence interface identified using X1 well rock core, survey to X1 well Well curvilinear characteristic is demarcated.In Fig. 4, the results showed that, log LLD, the LLS at X1 well high frequency sequence interface and calculating Log LLD, the LLS at the high frequency sequence interface of the variation characteristic and W1 well (Fig. 2) of three porosity curve and three holes of calculating Porosity curvilinear motion feature is similar, and compares well, shows that verifying is reliable.
After above-mentioned steps are feasible to identified low frequency sequence interface and high frequency sequence interface checking, repetition S110~ The step of S150 and full well section high frequency sequence interface identify is proved that feasible method promotes and applies in research area for above-mentioned, right The mesh stratum of other wells carries out the identification of full well section high frequency sequence interface, and under the understanding guidance of the depositional environment of S110, carries out high Frequency matching of successions.If Fig. 5 is using the above method, after carrying out full well section high frequency sequence division to Y1 well, and to W1 in research area The result that the high frequency sequence of well, X1 well and Y1 well compares.Four sections of area's thunder of research upper sub- section formation thickness variations are small, it is seen that bright The aobvious facies marker that the tidal-flat facies such as stone construction, laminated structure, bird's-eye structure deposition is laminated.Therefore, research area is mainly arid-half The formation of tidal rhythmites under arid background, high frequency sequence is mainly influenced by sea level variability.Moreover, high frequency sequence is being studied Longitudinally upper high frequency sequence number is identical in area, has good comparative, and high frequency sequence thickness is roughly the same transverse direction in, with synchronization Property, disclose these high frequency sequences not instead of by facies tract lateral transfer formed from cycle, sedimentation setting (sea level) variation It is formed by different cycle.It is carried out as a result, using fourth level sequence and Pyatyi sequence interface as the correlation marker layer of six grades of high frequency sequences Constraint, the method that six grades of high frequency sequences compare one by one, the high frequency sequence screen work of the final company's of foundation well.
To sum up, this specification is by taking the four sections of carbonate gas reservoirs of western thunderbolt ways of river as an example, to utilization core sample and rock core The identification at carbonate rock sequence interface is illustrated in sensitive log, solves that coring is few and representative difference is difficult to accurately Identify sequence interface, or only by core sample can not full well section identify sequence interface the problem of.This method strong operability, at This is low, meets sequence stratigraphy and principle is learned in well logging, can be convenient quick and accurately identifies different sequence interfaces, is carbon Sequence deposition understanding, the overall merit of Carbonate Rocks are had laid a good foundation.
Although describing the present invention herein with reference to specific embodiment, it should be understood that, these realities Apply the example that example is only principles and applications.It should therefore be understood that can be carried out to exemplary embodiment Many modifications, and can be designed that other arrangements, without departing from spirit of the invention as defined in the appended claims And range.It should be understood that different appurtenances can be combined by being different from mode described in original claim Benefit requires and feature described herein.It will also be appreciated that the feature in conjunction with described in separate embodiments can be used In other described embodiments.

Claims (8)

1. a kind of method for identifying carbonate rock sequence interface characterized by comprising
Step 1 demarcates at least one Logging Curves using the core sample of the carbonate rock;
Step 2 determines lithology sensitivity log from least one Logging Curves, wherein the lithology is sensitive Log includes Resistivity log combination, and the Resistivity log combination includes deep lateral resistivity log With shallow lateral Resistivity log;
Step 3 combines determining tri-porosity logging curve combination according to the Resistivity log, wherein three hole Spending log combination includes acoustic porosity log, density porosity log and neutron porosity log curve;
It is step 4, bent according to the lithology sensitivity log, the tri-porosity logging curve combination and the resistivity logging Line combination identification low frequency sequence interface, the low frequency sequence interface are fourth level sequence interface or the interface of its above rank sequence;
Step 5 identifies high frequency sequence interface according to the low frequency sequence interface, wherein the high frequency sequence interface is Pyatyi layer Sequence interface or six grades of sequence interfaces.
2. the method according to claim 1, wherein step 1 includes:
Step 11 makes core sample thin slice using the target zone core sample of the carbonate rock;
Step 12, the observation core sample thin slice simultaneously demarcate at least one Logging Curves.
3. according to the method described in claim 2, it is characterized in that, step 3 includes:
Step 31 combines determining compact limestone section according to the Resistivity log;
The interval transit time skeleton value and density skeleton value of step 32, the statistics compact limestone section, and according to the interval transit time Skeleton value and the density skeleton value calculate separately acoustic porosity and density porosity;
Step 33, the density porosity and neutron porosity of the acoustic porosity of calculating and calculating are placed in it is same Log road, to form the tri-porosity logging curve combination.
4. according to the method described in claim 3, it is characterized in that, step 3 further include:
Step 34 is corrected the tri-porosity logging curve combination using actual measurement core porosity.
5. according to the method described in claim 4, it is characterized in that, step 4 specifically includes:
The tri-porosity logging curve combination after correction is overlapped by step 41 at the compact limestone section;
Step 42, according to the lithology sensitivity log, the tri-porosity logging curve combination and the resistivity logging Curve combination, bond area drilling stratum, lithologic character data identify the low frequency sequence interface.
6. according to the method described in claim 5, it is characterized in that, step 5 specifically includes:
Step 51, the observation target zone core sample;
Lithology above and below step 52, the variation of lithological trend for summarizing high frequency sedimentary cycle and different high frequency sedimentary cycle contact surfaces Mutation differences identify the high frequency sequence interface.
7. according to the method described in claim 6, it is characterized in that, step 5 further include:
Step 53, according to the high frequency sequence interface, repeat step 2 and step 3 to identify full well section high frequency sequence interface.
8. method according to any one of claim 1 to 7, which is characterized in that the method also includes:
Step 6 verifies the low frequency sequence interface and the high frequency sequence interface.
CN201710507954.5A 2017-06-28 2017-06-28 Method for identifying carbonate rock sequence interface Active CN109143399B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710507954.5A CN109143399B (en) 2017-06-28 2017-06-28 Method for identifying carbonate rock sequence interface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710507954.5A CN109143399B (en) 2017-06-28 2017-06-28 Method for identifying carbonate rock sequence interface

Publications (2)

Publication Number Publication Date
CN109143399A true CN109143399A (en) 2019-01-04
CN109143399B CN109143399B (en) 2021-11-09

Family

ID=64803411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710507954.5A Active CN109143399B (en) 2017-06-28 2017-06-28 Method for identifying carbonate rock sequence interface

Country Status (1)

Country Link
CN (1) CN109143399B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111751903A (en) * 2019-03-29 2020-10-09 中国石油化工股份有限公司 Thin slice identification method for carbonate gas reservoir dissolved pore cause type
CN112965137A (en) * 2021-02-05 2021-06-15 中国石油天然气股份有限公司 Method and device for identifying lower cut valley in marine carbonate rock stratum
CN113433589A (en) * 2021-08-03 2021-09-24 成都理工大学 Weathered crust reservoir bottom interface identification method based on mathematical statistics
CN113885096A (en) * 2020-07-01 2022-01-04 中国石油化工股份有限公司 High-frequency sequence partitioning and small layer comparison method and device, electronic equipment and medium
CN114415238A (en) * 2022-01-25 2022-04-29 大庆油田有限责任公司 Method, device, equipment and medium for explaining marine carbonate earthquake in same period and different phases

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060290350A1 (en) * 2005-06-27 2006-12-28 Hursan Gabor G Method and apparatus for reservoir fluid characterization in nuclear magnetic resonance logging
CN103901471A (en) * 2014-03-31 2014-07-02 成都理工大学 Method for recovering depositional sequence prototype section

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060290350A1 (en) * 2005-06-27 2006-12-28 Hursan Gabor G Method and apparatus for reservoir fluid characterization in nuclear magnetic resonance logging
CN103901471A (en) * 2014-03-31 2014-07-02 成都理工大学 Method for recovering depositional sequence prototype section

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
刘家瑾等: "《煤田测井资料数字处理》", 31 March 1991 *
史云清等: "川西坳陷中段西部雷四段气藏沉积层序及其对储层的控制", 《东北石油大学学报》 *
姚光庆等: "《油气储层地质学实训教程》", 31 January 2016, 武汉:中国地质大学出版社 *
王锋庆: "火山岩油气藏的储层描述与建模", 《地下水》 *
石玉江: "低渗透岩性油藏含油性与富集区测井评价研究-以鄂尔多斯盆地姬塬地区长8油层组为例", 《中国博士学位论文全文数据库 基础科学辑》 *
范中专: "邵家油田沙四段湖相碳酸盐岩储层测井评价", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
裘亦楠等: "《油气储层评价技术(修订版)》", 30 April 2001, 石油工业出版社 *
邓宏文等: "《高分辨率层序地层学-原理及应用》", 30 September 2002, 北京:地质出版社 *
陈明江等: "含沥青储层的测井识别及评价", 《测井技术》 *
黄隆基: "《放射性测井原理》", 30 April 1985, 北京:石油工业出版社 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111751903A (en) * 2019-03-29 2020-10-09 中国石油化工股份有限公司 Thin slice identification method for carbonate gas reservoir dissolved pore cause type
CN113885096A (en) * 2020-07-01 2022-01-04 中国石油化工股份有限公司 High-frequency sequence partitioning and small layer comparison method and device, electronic equipment and medium
CN112965137A (en) * 2021-02-05 2021-06-15 中国石油天然气股份有限公司 Method and device for identifying lower cut valley in marine carbonate rock stratum
CN112965137B (en) * 2021-02-05 2023-01-10 中国石油天然气股份有限公司 Method and device for identifying lower cut valley in marine carbonate rock stratum
CN113433589A (en) * 2021-08-03 2021-09-24 成都理工大学 Weathered crust reservoir bottom interface identification method based on mathematical statistics
CN114415238A (en) * 2022-01-25 2022-04-29 大庆油田有限责任公司 Method, device, equipment and medium for explaining marine carbonate earthquake in same period and different phases

Also Published As

Publication number Publication date
CN109143399B (en) 2021-11-09

Similar Documents

Publication Publication Date Title
Winker et al. Sedimentary dynamics of the salt-dominated continental slope, Gulf of Mexico: integration of observations from the seafloor, near-surface, and deep subsurface
Sayago et al. Characterization of a deeply buried paleokarst terrain in the Loppa High using core data and multiattribute seismic facies classification
CN110275210A (en) A kind of recognition methods of the sedimentary micro facies model of carbonate rock high frequency sequence screen work
Campion et al. Outcrop expression of confined channel complexes
CN109143399A (en) A method of identification carbonate rock sequence interface
Ferreira et al. Unsupervised seismic facies classification applied to a presalt carbonate reservoir, Santos Basin, offshore Brazil
Al Haddad et al. Reservoir characterization, modeling, and evaluation of Upper Jurassic Smackover microbial carbonate and associated facies in Little Cedar Creek field, southwest Alabama, eastern Gulf coastal plain of the United States
Infante-Paez et al. Seismic analysis of paleotopography and stratigraphic controls on total organic carbon: Rich sweet spot distribution in the Woodford Shale, Oklahoma, USA
Amerman et al. Submarine mass-transport deposits of the Permian Cutoff Formation, West Texas, USA: Internal architecture and controls on overlying reservoir sand deposition
Neuhaus et al. Quantitative seismic reservoir characterization of an Oligocene Miocene carbonate buildup: Malampaya field, Philippines
Kumar et al. Delineation of a buried volcanic system in Kora prospect off New Zealand using artificial neural networks and its implications
Xu et al. How much systems-tract scale, three-dimensional stratigraphic variability is present in sequence stratigraphy?: An answer from the middle Miocene Pearl River Mouth Basin
Mayall et al. Seismic interpretation workflows for deep-water systems: A practical guide for the subsurface
Fournillon et al. Seismic expression of carbonate build-up karstification: karst modelling strategies and insights from synthetic seismic
Ahmad et al. High-resolution seismic and resistivity profiling of a buried Quaternary subglacial valley: Northern Alberta, Canada
Bueno et al. Structural and stratigraphic feature delineation and facies distribution using seismic attributes and well log analysis applied to a Brazilian carbonate field
Duguid et al. Geological controls on Upper Permian Plattendolomit Formation reservoir prospectivity, Wissey Field, UK Southern North Sea
Xia et al. Sequence architecture and depositional evolution in Lufeng Depression, Pearl River Mouth Basin, South China Sea Part B: The shore deposits of Zhuhai Formation
Lindzey et al. Lithological and petrophysical controls on production of the Mississippian limestone, northeastern Woods County, Oklahoma
Niyazi et al. Hectometer-scale, shallow buried honeycomb-like structures on the continental shelf of the Otway Basin, southeastern Australia
CN113514884A (en) Compact sandstone reservoir prediction method
Le Heron Interpretation of Late Ordovician glaciogenic reservoirs from 3-D seismic data: an example from the Murzuq Basin, Libya
Duvbiama et al. 3D Static Modelling of an offshore field in the Niger-delta
Dash et al. Play Fairway Analysis and Depositional Model of Panna Formation using High Frequency Sequence Stratigraphic Framework and Pre-Stack Geostatistical Inversion: Screening tools for Geological Risk Constrained Exploration
Dull Integrating sequence stratigraphy and multiple three-dimensional geostatistical realizations in constructing a model of the second Eocene reservoir, Wafra field, Partitioned Neutral Zone, Kuwait and Saudi Arabia

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