CN105089662B - Method and system for correcting carbonate rock stratum logging comprehensive histogram - Google Patents

Method and system for correcting carbonate rock stratum logging comprehensive histogram Download PDF

Info

Publication number
CN105089662B
CN105089662B CN201510441364.8A CN201510441364A CN105089662B CN 105089662 B CN105089662 B CN 105089662B CN 201510441364 A CN201510441364 A CN 201510441364A CN 105089662 B CN105089662 B CN 105089662B
Authority
CN
China
Prior art keywords
lithology
section
log
well
carbonate formation
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
Application number
CN201510441364.8A
Other languages
Chinese (zh)
Other versions
CN105089662A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510441364.8A priority Critical patent/CN105089662B/en
Publication of CN105089662A publication Critical patent/CN105089662A/en
Application granted granted Critical
Publication of CN105089662B publication Critical patent/CN105089662B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention provides a method and a system for correcting a carbonate stratum logging comprehensive histogram, wherein the method comprises the following steps: acquiring logging information of a carbonate rock stratum to be subjected to logging comprehensive histogram correction and lithology classification of a core section carbonate rock stratum; generating a target discrimination function by using a Fisher discrimination method according to lithology classification and logging information of the carbonate rock stratum at the coring section; and identifying the lithology of the carbonate stratum of the non-coring section of the carbonate stratum to be subjected to the logging comprehensive histogram correction by using the target discrimination function, and correcting the logging comprehensive histogram of the carbonate stratum to be subjected to the logging comprehensive histogram correction according to the identification result. The method fully utilizes the drilling coring data and the logging data, establishes the relationship between the lithology and the logging response, realizes the lithology re-identification based on Fisher discriminant analysis, corrects the original logging histogram with lower precision, and improves the lithology profile precision.

Description

The bearing calibration of carbonate formation well logging composite columnar section and system
Technical field
The present invention relates to oil-gas exploration technical field, specifically, being on a kind of carbonate formation well logging synthetic column The bearing calibration of figure and system.
Background technology
Stratum well logging composite columnar section is most basic data in petroleum geology research.Individual well sedimentary facies research relies on accurate Lithological profile, well logging is more accurate, single well facies division it is more accurate.At present, the acquisition of well logging depends on landwaste Well logging.And the well logging work of situ of drilling well is limited by time, experience, environment, the original well logging of acquisition often low precision, It is even more so for the less carbonate formation of lithology difference.At present, carbonate formation well logging lithology is very It is single, the content and variation of lithological of carbonate rock particle, grey shale, crystal grain can not be reflected, constrain sedimentary facies and SEQUENCE STRATIGRAPHIC Learn research.Although landwaste thin slice can be corrected to well logging, time-consuming and expense is expensive, and depth playback has error.It is right In the well that Lao Jing, mineral rights change, even more it is difficult by landwaste and carries out well logging correction.
The content of the invention
The invention provides a kind of bearing calibration of carbonate formation well logging composite columnar section, methods described includes following Step:
S100, obtain the well-log information of the carbonate formation of pending well logging composite columnar section correction and a section carbon of coring The lithology breakdown on Carbonate Rocks stratum;
S200, according to the lithology breakdown of section carbonate formation of coring and the well-log information, sentenced using Fisher Other method generates target-recognition function;
S300, using the target-recognition function, the carbonate rock corrected to the pending well logging composite columnar section The lithology of the carbonate formation of the section of not coring of layer is identified, and is integrated according to the recognition result correction pending well logging The well logging composite columnar section of the carbonate formation of block diagram correction.
Correspondingly, the present invention also provides a kind of correction system of carbonate formation well logging composite columnar section, including obtains Module, generation module and correction module;
The acquisition module, the well-log information of the carbonate formation for obtaining pending well logging composite columnar section correction And the lithology breakdown for section carbonate formation of coring;
The generation module, lithology breakdown and the well-log information for section carbonate formation of being cored according to, Target-recognition function is generated using Fisher method of discrimination;
The correction module, for utilizing the target-recognition function, the pending well logging composite columnar section is corrected The lithology of carbonate formation of section of not coring of carbonate formation be identified, and treated according to correcting recognition result Carry out the well logging composite columnar section of the carbonate formation of well logging composite columnar section correction.
The effect of the present invention:The survey for the carbonate formation that the present invention is corrected by obtaining pending well logging composite columnar section The lithology breakdown of well data and section carbonate formation of coring, then according to core section carbonate formation lithology breakdown and Well-log information, target-recognition function is generated using Fisher method of discrimination, finally using target-recognition function, to pending well logging The lithology of the carbonate formation of the section of not coring of the carbonate formation of composite columnar section correction is identified, and according to identification The well logging composite columnar section of the carbonate formation of the pending well logging composite columnar section correction of calibration of the output results.It takes full advantage of brill Well coring data and well-log information, the relation between lithology and log response is established, realized based on Fisher discriminant analyses Lithology identify that the well logging block diagram relatively low to original precision is corrected, and improves lithological profile precision again.
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 There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those skilled in the art, without having to pay creative labor, can be with root Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the flow chart of the bearing calibration of carbonate formation well logging composite columnar section provided in an embodiment of the present invention;
Fig. 2 is the structure chart of the correction system of carbonate formation well logging composite columnar section provided in an embodiment of the present invention;
Fig. 3 is the concrete structure schematic diagram of the correction system of the carbonate formation well logging composite columnar section shown in Fig. 2;
Fig. 4 is in the embodiment of the present invention and 26, the 4 shaft bottom heart core log composite columnar section;
Fig. 5 is in the embodiment of the present invention and 4 well Cambrian system Qiulitag group carbonate formations pending corrections The well logging composite columnar section of carbonate formation;
Fig. 6 is in the embodiment of the present invention and 4 well Cambrian system Qiulitag group carbonate formation logs corrections Result schematic diagram;
Fig. 7 is the knowledge that in the embodiment of the present invention and 4 well Cambrian system Qiulitag group carbonate rocks identify sample set Other effect diagram;
Fig. 8 is comprehensive for after the correction in the embodiment of the present invention and 4 well Cambrian system Qiulitag group carbonate formation well loggings Close block diagram.
Embodiment
Coordinate schema and presently preferred embodiments of the present invention below, the present invention is expanded on further to reach predetermined goal of the invention institute The technological means taken.
As shown in figure 1, the embodiments of the invention provide a kind of bearing calibration of carbonate formation well logging composite columnar section, Comprise the following steps:
S100, obtain the well-log information of the carbonate formation of pending well logging composite columnar section correction and a section carbon of coring The lithology breakdown on Carbonate Rocks stratum;
S200, according to the lithology breakdown and well-log information of section carbonate formation of coring, given birth to using Fisher method of discrimination Into target-recognition function;
S300, using target-recognition function, the carbonate formation corrected to pending well logging composite columnar section does not take The lithology of the carbonate formation of heart section is identified, and corrects what pending well logging composite columnar section corrected according to recognition result The well logging composite columnar section of carbonate formation.
Wherein, step S100 comprises the following steps:
S110, obtain the log of the carbonate formation of pending well logging composite columnar section correction;
S120, the lithology for section carbonate formation of coring is classified using Deng Hamu sorting techniques.
Step S120 mainly accurately identifies to drilling and coring delivery section lithology.By core observation and description, rock is carried out Property it is preliminary name, be sampled for the rock core section of different lithology, sampling interval is advisable with 0.5m.It is frequent for variation of lithological Well section, should suitably reduce sampling interval.Find the interface of variation of lithological accurately by core observation, draw the lithology column of rock core Figure.The rock sample of collection is prepared into petrographic thin section, accurately named for lithology under an optical microscope.Name mode and use Deng The classification schemes that Ha Mu (Dunham, 1962) is proposed, are divided into CRYSTALLINE ROCKS, particle rock, shale by carbonate rock by the difference of structure Particle rock, particle matter mud stone, mud stone.When carbonate formation lithology is based on CRYSTALLINE ROCKS, crystal grain can be according to its granularity division For gravel crystalline substance, sand crystal, powder crystal, micrite.
Step S200 comprises the following steps:
S210, the quality of the log of acquisition is judged and corrected;
S220, Lithology Discrimination sample is established according to the log after the lithology breakdown for section carbonate formation of coring and correction This collection;
S230, using the lithology classification for section carbonate formation of coring as packet variable, with the survey of the log after correction Well value is independent variable, and discriminant function is generated using Fisher method of discrimination;
S240, using Lithology Discrimination sample set, the Lithology Discrimination accuracy of discriminant function is examined, is chosen from discriminant function Meet the target-recognition function of setting Lithology Discrimination accuracy requirement.
Further, step S210 comprises the following steps:
S211, when the hole diameter value of log exceedes baseline value 50%, and other logs also go out in corresponding well section When being now mutated, judge that there occurs expanding for current well section;
S212, when the expanding well section thickness of generation is less than 20m, utilize the value pair of the log of adjacent non-expanding well section It is corrected;
S213, when expanding well section thickness occurring being more than 20m, reject the lithology classification of the expanding well section of generation with it is corresponding Log value.
For the quality of log except unexpected by instrumental effects, wall quality is also its important influence factor.If The borehole wall occurs significantly to collapse in drilling process, and in well logging, logger can not be adjacent to the borehole wall well, cause to survey The data value obtained can not reflect real data value.In actual applications, if calliper log curve shows that expanding surpass occurs for well section Bit diameter 50% (the hole diameter value of log exceedes baseline value 50%) is crossed, and form occur accordingly in other logs During mutation, show that the log quality at this is low, it is necessary to correct.When the thickness of continuous extension diameter section is less than 20m, it is believed that Formation variation is little in the thickness range, and it is corrected using the log value of adjacent non-expanding well section.Expand when continuous Footpath section is more than 20m, then is difficult to ensure that the homogeneity on stratum, can not accurately be corrected.Expanding well section, which occurs, in these to enter The correction of row Lithology Discrimination and well logging block diagram, coring data can not be brought into Lithology Discrimination sample set.
Further, step S220 comprises the following steps:
S221, according to the lithology classification for section carbonate formation of coring, select the well logging after correspondingly 5 to 7 corrections bent Line, build Lithology Discrimination sample set;
S222, judge to be less than 10 with the presence or absence of the depth point of the corresponding log of lithology in Lithology Discrimination sample set Individual sample point, if in the presence of Rejection of samples point.
Specifically, using lithology as the log value of class label li (i=1,2 ..., k) and its corresponding log (X1, X2 ..., Xm) establish lithology sample set (li, X1, X2 ..., Xm).Every mouth well establishes single Lithology Discrimination sample set.If rock The cumulative thickness of certain a kind of lithology is too thin in the heart or to cause certain a kind of lithology to be difficult to because logging data quality is low corresponding enough Log value (being less than 10 log values), then such lithology can not participate in the calculating of Lithology Discrimination.
Log should try one's best the selection curve sensitive to petrophysical property, for example, density log, acoustic travel time logging, The serial well logging of gamma ray log, resistivity, neutron porosity log etc..In order to farthest utilize well-log information, Ying Xuan Select 5 to 7 logs and form Lithology Discrimination sample set with corresponding lithology classification.
In Lithology Discrimination sample set, the well logging of resistivity series is taken the logarithm value, and the unit of all kinds of logs is made in Great Britain Or the metric system.
Further, step S240 comprises the following steps:
S241, Lithology Discrimination sample set is updated to discriminant function, is grouped variate-value by calculating, obtains Lithology Discrimination sample The lithology classification of this collection;
S242, according to the lithology classification of obtained Lithology Discrimination sample set and the lithology breakdown for section carbonate formation of coring As a result, the Lithology Discrimination accuracy of current discriminant function is calculated;
S243, when the Lithology Discrimination accuracy of discriminant function is more than or equal to 90%, chooses recognition function and sentence for target Other function;
S244, when the Lithology Discrimination accuracy of discriminant function is less than 90%, reject the minimum rock of Lithology Discrimination accuracy Property classification, perform step S230.
The general principle of Fisher discriminant analyses is to find most suitable axis of projection and dimension so that group difference is maximum, Group difference is minimum.Realize that dimension about subtracts by Fisher discriminant analyses, obtain n (m<N) individual typical discriminant function, makes it both protect The information of original m parameter is stayed, and can realizes effective identification of unfiled sample.This method integrated application contains much information, Differentiate that the degree of accuracy is high in the case where data volume is abundant.
As a kind of embodiment, first lithology is represented with digital code, as grainstone is 1, packstone 2, Micrite is 3, and Lithology Discrimination sample set data then are imported into SPSS systems.Under SPSS system discriminant analysis modules, Lithology row are defined as to be grouped variable, the maxima and minima of definition packet variable is respectively the maximum of lithology digital code With minimum value, log value row are defined as independent variable.The differentiation letter of Lithology Discrimination sample can be obtained by being calculated using discriminant analysis Number.For example, when the classification accuracy rate of discriminant function is less than 90%, it is good to show that discriminant function can not be carried out to all lithology Classification, a kind of minimum lithology of accuracy is rejected.Then Lithology Discrimination is re-started using SPSS systems, obtains new sentencing Other function, until lithology breakdown accuracy is more than 90%.
After obtaining the discriminant function of classification accuracy rate more than 90%, the log value on stratum to be corrected is inputted into SPSS systems, Wherein lithology row vacancy, the identification of lithology, the lithology classification predicted, according to the lithology of prediction are carried out using same method The well logging composite columnar section for the carbonate formation that classification corrects to pending well logging composite columnar section is corrected.
The bearing calibration of carbonate formation well logging composite columnar section provided in an embodiment of the present invention, can make full use of brill Well coring data and well-log information, the related pass established using the method for Fisher discriminant analyses between lithology and log response System, the Lithology Discrimination of carbonate rock is realized, so as to make correction to well logging composite columnar section, improve lithological profile precision.
Correspondingly, referring to Fig. 2, based on same inventive concept, the present invention also provides a kind of carbonate formation well logging synthesis The correction system of block diagram, including acquisition module 100, generation module 200 and correction module 300;
Acquisition module 100 be used to obtain the well-log information of the carbonate formation of pending well logging composite columnar section correction with And the lithology breakdown for section carbonate formation of coring;
Generation module 200 is used for lithology breakdown and well-log information according to section carbonate formation of coring, utilizes Fisher Method of discrimination generates target-recognition function;
Correction module 300 is used to utilize target-recognition function, to the carbonate rock of pending well logging composite columnar section correction The lithology of the carbonate formation of the section of not coring on stratum is identified, and corrects pending well logging according to recognition result and integrate post The well logging composite columnar section of the carbonate formation of shape figure correction.
Wherein, acquisition module 100 includes acquiring unit 110 and taxon 120, acquiring unit 110 be used to obtain treat into The log of the carbonate formation of row well logging composite columnar section correction, taxon 120, for using Deng Hamu classification sides Method is classified to the lithology for section carbonate formation of coring.
Referring to Fig. 3, in the correction system of carbonate formation well logging composite columnar section provided by the invention, generation module 200 establish unit 220, generation unit 230 including judging unit 210, sample set and choose unit 240.
Judging unit 210 is used to the quality of the log of acquisition is judged and corrected;
Sample set establish unit 220 be used for according to core section carbonate formation lithology breakdown and correction after well logging song Line establishes Lithology Discrimination sample set;
Generation unit 230 is used for using the lithology classification for section carbonate formation of coring as packet variable, with the survey after correction The log value of well curve is independent variable, and discriminant function is generated using Fisher method of discrimination;
Choose unit 240 to be used to utilize Lithology Discrimination sample set, the Lithology Discrimination accuracy of discriminant function is examined, from differentiation The target-recognition function for meeting the accuracy requirement of setting Lithology Discrimination is chosen in function.
Further, judging unit 210 includes the correction syndrome of subelement 212 and second of judgment sub-unit 211, first Unit 213.
Judgment sub-unit 211 is used to exceed baseline value 50% in the hole diameter value of log, and other logs also exist When corresponding well section is mutated, judge that there occurs expanding for current well section;
First correction subelement 212 is used to, when the expanding well section thickness of generation is less than 20m, utilize adjacent non-expanding well section The value of log it is corrected;
Second correction subelement 213 is used to, when the expanding well section thickness of generation is more than 20m, reject and expanding well section occurs Lithology classification and corresponding log value.
Further, sample set, which establishes unit 220, includes the structure correction subelement 222 of subelement 221 and the 3rd.Structure Unit 221 is used for the lithology classification according to section carbonate formation of coring, and selects the log after correspondingly 5 to 7 corrections, Build Lithology Discrimination sample set;3rd correction subelement 222 is used to judge in Lithology Discrimination sample set with the presence or absence of lithology and its The depth point of corresponding log is less than the sample point of 10, if in the presence of Rejection of samples point.
Further, choosing unit 240 includes the firstth computation subunit 241, the second computation subunit 242, chooses son Unit 243 and rejecting subelement 244.
First computation subunit 241 is used to Lithology Discrimination sample set being updated to discriminant function, and variable is grouped by calculating Value, obtain the lithology classification of Lithology Discrimination sample set;
Second computation subunit 242 is used for lithology classification and a section carbonate of coring according to obtained Lithology Discrimination sample set The lithology breakdown result on rock stratum, calculate the Lithology Discrimination accuracy of current discriminant function;
Subelement 243 is chosen to be used to, when the Lithology Discrimination accuracy of discriminant function is more than or equal to 90%, choose identification Function is target-recognition function;
Subelement 244 is rejected to be used to, when the Lithology Discrimination accuracy of discriminant function is less than 90%, reject Lithology Discrimination just The minimum lithology classification of true rate.
In the correction system of carbonate formation well logging composite columnar section provided by the invention, log is can be anti- The log of petrophysical property is reflected, it is bent density log curve, acoustic travel time logging curve, gamma ray log can be included Line, resistivity series log and neutron porosity log curve etc..
The present invention takes full advantage of drilling and coring delivery data and well-log information, establishes the pass between lithology and log response System, realizes the lithology based on Fisher discriminant analyses and identifies again, the well logging block diagram relatively low to original precision has carried out school Just, lithological profile precision is improved.
In order to which the bearing calibration of carbonate formation well logging composite columnar section provided the invention described above and system are carried out Apparent explanation, illustrated with reference to a specific embodiment, however, it should be noted that the embodiment is only It is in order to which the present invention is better described, does not form and the present invention is improperly limited.
Tarim Basin Cambrian system extensive development Marine Carbonate Rocks deposition, and generally by white clouds lithification.Cold force at present The Songliao basin of system is very low, and drilling hole number is few.In terms of the well logging composite columnar section of stratum, the lithology base of Carbonate Rocks of The Cambrian This is all named as dolomite.As Bachu Area and 4 well Cambrian system Qiulitag groups develop more than 730 rice carbonate sediment, Lithology is very single on the well logging of stratum, and in addition to 50 meters of thick gabbros are developed, remaining is carbonate rock, including 35 The thick micrite of rice.Single lithological profile is all unfavorable for sedimentary facies, SEQUENCE STRATIGRAPHIC and reservoir study.The well does not have Have and drilling and coring delivery is carried out to limestone section, only dolomite section is cored.From coring, dolomitic lithology classification is enriched. The well has carried out the well logging of a variety of logging programs simultaneously, therefore, can utilize the lithology of rock core and well-log information to dolomite section Identification, so as to being corrected to stratum well logging composite columnar section.
Cored altogether 4 times with 4 well Cambrian system Qiulitag groups, add up 17.87 meters of drilling depth, it is accumulative to core 16.99 meters.With rock The mode of heart observation carries out detailed observation to rock core and described, and the lithology of rock core is named and describes the change sequence of lithology Row.Rock core is sampled every 20cm, prepares core wafer.Observation by light microscope thin slice is utilized in laboratory, to lithology Progress is final to be named, and draws core log composite columnar section (as shown in Figure 4).Observed by rock core, thin slice, dolomite is Crystal dolomite, dolomicrite, powder crystal dolomite, fine grain dolomite, middle brilliant dolomite can be subdivided into.
Analyze the logging data quality for the well section to be corrected.It is to judge the borehole wall according to the change of calliper log curve It is no to collapse.When hole diameter value is more than baseline value more than 50%, while other logs are also in corresponding depth value occurrence value Mutation, then it is assumed that the logging data quality of the well section is low, be unsatisfactory for Lithology Discrimination prediction requirement.For continuous expanding well section Thickness is less than 20m, it is believed that lithology, without obvious change occurs, utilizes adjacent non-expanding well section in this depth bounds Log value is corrected to it.When continuous expanding well section thickness is more than 20m, should by the lithology classification of the well section with it is corresponding Log value is rejected, and is not brought into Lithology Discrimination sample set, while the identification and correction of lithology are no longer carried out to the well section. It is 9.6in, wherein 4475.875m to 4476.625m as knowable to the bit size with 4 well 4470m-4477m by calliper log The borehole wall occurs significantly to collapse, expanding more than 50%, while also incidence is significantly mutated other logs, is thus believed that 4475.875m to 4476.625m logging data quality is low, it is necessary to be corrected (as shown in Figure 5).In operation, it is necessary to right The well-log information of the well section of whole pending well logging block diagram correction is all evaluated, and (such as Fig. 6 is corrected or rejected by evaluation result It is shown).
More than 20 kind well loggings are carried out altogether with 4 wells, have selected to carry out lithology to 5 sensitive logs of petrophysical property Identification and correction, including density log (DEN), acoustic travel time logging (AC), gamma ray log (GR), deep lateral resistivity survey Well (RD), neutron porosity log (CNL).Lithology Discrimination sample set is established, is shown in Table one, wherein being logged well for deep lateral resistivity Value is taken the logarithm.Notice that the middle dolomitic thickness of crystalline substance is only 1 meter in rock core, every kind of log has only corresponded to 8 log values.By In the available log data amount of this kind of lithology very little, therefore fine grain dolomite and the dolomitic lithology of middle crystalline substance can be merged and are referred to as Detail crystalline substance dolomite.In Lithology Discrimination sample, every log of dolomicrite has 46 sample points, and powder crystal dolomite is then There are 42 sample points, detail crystalline substance dolomite shares 32 sample points.Lithology is represented with digital code, dolomicrite 1, powder Brilliant dolomite is 2, detail crystalline substance dolomite is 3.
Table one
Mud Powder It is brilliant It is deep Degree GR LnR D DE N CN L AC Powder Carefully It is brilliant It is deep Degree GR LnR D DE N CN L AC Carefully In It is brilliant It is deep Degree GR LnR D DE N CN L AC
1 435 4.8 16. 011 5.3 028 2. 46 5 0. 03 5 49. 36 2 492 4 11. 534 4.1 693 2. 68 2 0. 22 8 45. 317 3 481 5.9 25. 326 5.9 65 2. 78 4 0. 04 8 44.3 36
1 435 4.9 16. 441 5.6 007 2. 38 1 0. 03 9 48. 217 2 492 4.1 12. 75 4.2 123 2. 72 5 0. 22 7 45. 126 3 481 6 22. 835 5.5 878 2. 79 2 0. 04 8 44.3 76
1 435 5 17. 12 5.8 503 2. 32 3 0. 04 3 46. 93 2 492 4.3 13. 739 4.2 467 2. 75 9 0. 22 7 45. 241 3 481 6.1 20. 355 5.3 336 2. 80 1 0. 04 8 44.3 81
1 435 5.1 18. 259 6.0 334 2. 30 2 0. 04 4 45. 843 2 492 4.4 14. 2 4.2 644 2. 78 2 0. 23 45. 539 3 481 6.3 18. 306 5.1 759 2. 80 7 0. 04 8 44.3 51
1 435 5.3 19. 625 6.1 811 2. 30 8 0. 04 6 45. 258 2 492 4.5 13. 99 4.2 64 2. 79 2 0. 23 3 45. 869 3 481 6.4 16. 941 5.1 074 2. 81 1 0. 04 7 44.2 85
1 435 5.4 20. 679 6.3 152 2. 33 2 0. 05 45. 127 2 492 4.6 13. 188 4.2 531 2. 79 7 0. 23 5 46. 081 3 481 6.5 16. 356 5.1 077 2. 81 2 0. 04 7 44.2
1 435 5.5 20. 789 6.4 285 2. 35 5 0. 05 6 45. 293 2 492 4.8 12. 198 4.2 48 2. 80 7 0. 23 7 46. 141 3 481 3.8 24. 867 6.4 735 2. 82 3 0. 02 8 44.0 04
1 435 5.6 19. 339 6.4 914 2. 35 9 0. 06 3 45. 515 2 492 4.9 11. 275 4.2 56 2. 82 1 0. 23 7 46. 085 3 481 3.9 25. 143 6.6 164 2. 82 0. 03 1 44.1 21
1 435 5.8 16. 827 6.4 983 2. 34 1 0. 06 9 45. 662 2 492 5 10. 597 4.2 684 2. 83 5 0. 23 8 45. 99 3 481 4 24. 684 6.6 664 2. 81 9 0. 03 4 44.2 44
1 435 5.9 13. 956 6.4 521 2. 31 0. 07 3 45. 662 2 492 5.1 10. 345 4.2 505 2. 84 9 0. 24 45. 939 3 481 4.1 23. 704 6.6 019 2. 81 9 0. 03 7 44.3 5
1 435 6 11. 485 6.3 581 2. 27 8 0. 07 6 45. 517 2 492 5.3 10. 452 4.1 758 2. 86 0. 24 3 45. 956 3 481 4.3 22. 585 6.4 887 2. 81 9 0. 03 9 44.4 2
1 435 6.1 9.9 97 6.2 254 2. 26 6 0. 07 6 45. 493 2 492 5.4 10. 774 4.0 405 2. 86 5 0. 24 6 46. 016 3 481 4.4 21. 524 6.3 679 2. 81 6 0. 04 44.4 44
1 435 6.3 9.3 16 6.0 797 2. 27 9 0. 07 6 45. 91 2 492 5.5 11. 124 3.8 683 2. 86 2 0. 24 2 46. 06 3 481 4.5 20. 604 6.2 597 2. 81 4 0. 04 44.4 13
1 435 6.4 9.2 66 5.9 443 2. 31 6 0. 07 6 46. 409 2 492 5.6 11. 297 3.7 063 2. 84 8 0. 23 6 46. 046 3 481 4.6 19. 889 6.1 744 2. 81 0. 04 44.3 22
1 435 6.5 9.7 08 5.8 39 2. 36 5 0. 07 6 46. 514 2 492 5.8 11. 242 3.5 794 2. 82 8 0. 23 46. 006 3 481 4.8 19. 461 6.1 135 2. 80 6 0. 04 44.1 93
1 435 6.6 10. 622 5.7 811 2. 40 2 0. 07 5 46. 259 2 481 6.6 16. 604 5.1 348 2. 81 1 0. 04 8 44. 148 3 447 7.5 14. 712 2.9 606 2. 74 4 0. 04 1 48.6 04
1 435 6.8 11. 809 5.7 617 2. 41 1 0. 07 2 46. 165 2 481 6.8 17. 298 5.1 206 2. 81 0. 04 9 44. 182 3 447 7.6 18. 331 3.1 367 2. 77 8 0. 04 1 47.3 08
1 435 6.9 13. 041 5.7 769 2. 39 4 0. 06 6 46. 252 2 481 6.9 18. 045 5.0 22 2. 81 0. 05 44. 347 3 447 7.8 22. 803 3.3 484 2. 79 2 0. 04 1 46.9 84
1 435 7 14. 004 5.8 266 2. 36 6 0. 06 46. 347 2 481 7 18. 483 4.8 411 2. 81 1 0. 05 1 44. 657 3 447 7.9 27. 851 3.5 583 2. 79 9 0. 04 3 47.2 4
1 435 7.1 14. 263 5.9 272 2. 35 3 0. 05 4 46. 352 2 481 7.1 18. 467 4.6 352 2. 81 5 0. 05 1 45. 081 3 447 8 32. 879 3.7 184 2. 80 3 0. 04 4 47.7 45
1 435 7.3 13. 809 6.0 812 2. 36 3 0. 05 46. 185 2 481 7.3 18. 432 4.4 599 2. 81 9 0. 05 1 45. 498 3 447 8.1 36. 574 3.7 878 2. 80 3 0. 04 5 48.0 6
1 435 7.4 12. 787 6.2 889 2. 39 3 0. 05 45. 995 2 481 7.4 18. 849 4.3 432 2. 82 1 0. 05 45. 807 3 447 8.3 37. 962 3.8 083 2. 80 1 0. 04 4 48.0 17
1 435 7.5 11. 582 6.5 3 2. 43 0. 05 5 46. 021 2 481 7.5 20. 001 4.2 913 2. 82 0. 04 8 45. 938 3 447 8.4 36. 572 3.8 316 2. 8 0. 04 1 47.7 14
1 435 7.6 10. 756 6.7 423 2. 45 9 0. 06 4 46. 187 2 481 7.6 21. 629 4.2 929 2. 81 6 0. 04 5 45. 854 3 447 8.5 32. 828 3.8 916 2. 80 1 0. 03 8 47.4 01
1 435 7.8 10. 57 6.8 308 2. 46 9 0. 07 3 46. 233 2 481 7.8 22. 826 4.2 984 2. 81 2 0. 04 2 45. 663 3 447 8.6 28. 137 4.0 056 2. 80 3 0. 03 6 47.3 28
1 435 7.9 11. 014 6.7 686 2. 45 3 0. 08 46. 211 2 481 7.9 22. 906 4.2 654 2. 81 0. 04 45. 478 3 447 8.8 24. 172 4.1 445 2. 80 5 0. 03 3 47.3 77
1 435 8 11. 805 6.6 069 2. 41 4 0. 08 46. 254 2 481 8 21. 744 4.1 762 2. 81 0. 03 9 45. 383 3 447 8.9 21. 57 4.2 762 2. 80 4 0. 03 1 47.3 2
1 435 8.1 12. 554 6.4 449 2. 36 1 0. 07 6 46. 356 2 481 8.1 19. 972 4.0 578 2. 80 9 0. 03 9 45. 41 3 447 9 20. 186 4.3 732 2. 80 2 0. 03 2 47.0 19
1 435 8.3 13. 051 6.3 455 2. 30 8 0. 07 46. 475 2 481 8.3 18. 691 3.9 609 2. 80 9 0. 03 9 45. 498 3 447 9.1 19. 553 4.4 232 2. 8 0. 03 6 46.5 59
1 435 8.4 13. 352 6.3 38 2. 26 0. 06 5 46. 624 2 481 8.4 18. 43 3.9 141 2. 80 8 0. 03 9 45. 592 3 447 9.3 19. 234 4.4 419 2. 79 8 0. 04 2 46.1 47
1 435 8.5 13. 645 6.4 293 2. 22 3 0. 06 2 46. 752 2 481 8.5 19. 11 3.9 241 2. 80 9 0. 03 9 45. 651 3 447 9.4 19. 264 4.4 238 2. 80 1 0. 05 1 45.9 09
1 435 8.6 13. 976 6.6 088 2. 20 3 0. 06 47. 011 2 481 8.6 20. 148 3.9 823 2. 80 9 0. 03 8 45. 646 3 447 9.5 19. 738 4.3 314 2. 80 5 0. 06 46.0 5
Lithology Discrimination sample set is loaded into data statistics system SPSS, under discriminant analysis module, lithology is arranged and defined To be grouped variable, the maxima and minima of definition packet variable is respectively the maxima and minima of lithology digital code, will Log value row are defined as independent variable, are shown in Table two.
Table two
By Fisher discriminant analyses, following two discriminant functions can be obtained:
F1=0.97GR-0.329RD+15.467DEN+5.587CNL-0.102AC-36.258
F2=0.171GR+1.107RD+2.697DEN-2.560CNL+0.552AC-41.082
Lithology breakdown, classification accuracy rate 92.9% are carried out to Lithology Discrimination sample set using above-mentioned discriminant function, it is believed that The discriminant function being calculated is effective, and the lithology of other well sections can be identified (as shown in Figure 7).Wherein micrite is white Yun Yan classification accuracy rate 100%, the dolomitic classification accuracy rate of powder crystal are 89.7%, the brilliant dolomitic classification accuracy rate of detail For 85.2%, three are shown in Table.The log value for treating Lithology Discrimination stratomere is inputted into SPSS systems, it is right using identical discriminant function Lithology is predicted.Finally, carbonate formation well logging composite columnar section to be corrected is carried out using the lithology result of prediction Correction, the well logging composite columnar section after being corrected, as shown in Figure 8.
Table three
The present invention is by obtaining the well-log information of the carbonate formation of pending well logging composite columnar section correction and taking The lithology breakdown of heart section carbonate formation, then according to the lithology breakdown and well-log information of section carbonate formation of coring, profit Target-recognition function is generated with Fisher method of discrimination, finally using target-recognition function, to pending well logging composite columnar section The lithology of the carbonate formation of the section of not coring of the carbonate formation of correction is identified, and is treated according to recognition result correction Carry out the well logging composite columnar section of the carbonate formation of well logging composite columnar section correction.It takes full advantage of drilling and coring delivery data With well-log information, the relation between lithology and log response is established, the lithology based on Fisher discriminant analyses is realized and knows again Not, the well logging block diagram relatively low to original precision is corrected, and improves lithological profile precision.
Particular embodiments described above, the purpose of the present invention, technical scheme and beneficial effect are carried out further in detail Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., should be included in this Within the protection domain of invention.

Claims (10)

1. a kind of bearing calibration of carbonate formation well logging composite columnar section, it is characterised in that comprise the following steps:
S100, obtain the well-log information of the carbonate formation of pending well logging composite columnar section correction and a section carbonate of coring The lithology breakdown on rock stratum;
S200, according to the lithology breakdown of section carbonate formation of coring and the well-log information, utilize Fisher differentiation sides Method generates target-recognition function;
S300, using the target-recognition function, to the carbonate formation of the pending well logging composite columnar section correction The lithology of the carbonate formation for section of not coring is identified, and corrects the pending well logging synthetic column according to recognition result Scheme the well logging composite columnar section of the carbonate formation of correction;
Wherein, the step S100 comprises the following steps:
S110, obtain the log of the carbonate formation of the pending well logging composite columnar section correction;
S120, the lithology for coring section carbonate formation is classified using Deng Hamu sorting techniques;
Wherein, the step S200 comprises the following steps:
S210, the quality of the log of acquisition is judged and corrected;
S220, lithology is established according to the log after the lithology breakdown of section carbonate formation of coring and correction and known Other sample set;
S230, using the lithology classification of section carbonate formation of coring as packet variable, with the log after correction Log value be independent variable, utilize Fisher method of discrimination generation discriminant function;
S240, using the Lithology Discrimination sample set, the Lithology Discrimination accuracy of the discriminant function is examined, from the differentiation letter The target-recognition function for meeting the accuracy requirement of setting Lithology Discrimination is chosen in number.
2. the bearing calibration of carbonate formation well logging composite columnar section according to claim 1, it is characterised in that described Step S210 comprises the following steps:
S211, when the hole diameter value of the log exceedes baseline value 50%, and other logs also go out in corresponding well section When being now mutated, judge that there occurs expanding for current well section;
S212, when the expanding well section thickness of generation is less than 20m, utilize the value pair of the log of adjacent non-expanding well section It is corrected;
S213, when the expanding well section thickness of generation is more than 20m, reject the lithology classification that expanding well section occurs and corresponding survey The value of well curve.
3. the bearing calibration of carbonate formation well logging composite columnar section according to claim 1, it is characterised in that step S220 comprises the following steps:
S221, according to the lithology classification of section carbonate formation of coring, select the survey after correspondingly 5 to 7 corrections Well curve, build the Lithology Discrimination sample set;
S222, judge few with the presence or absence of the depth point of the corresponding log of lithology in the Lithology Discrimination sample set In the sample point of 10, if in the presence of rejecting the sample point.
4. the bearing calibration of carbonate formation well logging composite columnar section according to claim 1, it is characterised in that described Step S240 comprises the following steps:
S241, the Lithology Discrimination sample set is updated to the discriminant function, is grouped variate-value by calculating, obtains the rock Property identification sample set lithology classification;
S242, according to the lithology classification of the obtained Lithology Discrimination sample set and the lithology of section carbonate formation of coring Classification results, calculate the Lithology Discrimination accuracy of presently described discriminant function;
S243, when the Lithology Discrimination accuracy of the discriminant function is more than or equal to 90%, it is mesh to choose the recognition function Mark discriminant function;
S244, when the Lithology Discrimination accuracy of the discriminant function is less than 90%, reject the minimum rock of Lithology Discrimination accuracy Property classification, perform step S230.
5. the bearing calibration of the carbonate formation well logging composite columnar section according to any one of Claims 1-4, its feature It is, the log is that can reflect the log of petrophysical property, including density log curve, interval transit time are surveyed Well curve, Natural Gamma-ray Logging Curves, resistivity series log and neutron porosity log curve.
6. a kind of correction system of carbonate formation well logging composite columnar section, it is characterised in that including acquisition module, generation mould Block and correction module;
The acquisition module, for the carbonate formation that obtains the correction of pending well logging composite columnar section well-log information and Core the lithology breakdown of section carbonate formation;
The generation module, lithology breakdown and the well-log information for section carbonate formation of being cored according to, is utilized Fisher method of discrimination generates target-recognition function;
The correction module, for utilizing the target-recognition function, to the carbon of the pending well logging composite columnar section correction The lithology of the carbonate formation of the section of not coring on Carbonate Rocks stratum is identified, and described pending according to recognition result correction The well logging composite columnar section of the carbonate formation of well logging composite columnar section correction;
Wherein, the acquisition module includes acquiring unit and taxon;
The acquiring unit, it is bent for obtaining the well logging of carbonate formation of the pending well logging composite columnar section correction Line;
The taxon, for being divided using Deng Hamu sorting techniques the lithology for coring section carbonate formation Class;
Wherein, the generation module establishes unit, generation unit including judging unit, sample set and chooses unit;
The judging unit, the quality for the log to acquisition are judged and corrected;
The sample set establishes unit, described in after the lithology breakdown for section carbonate formation of being cored according to and correction Log establishes Lithology Discrimination sample set;
The generation unit, for using the lithology classification of section carbonate formation of coring as packet variable, after correction The log value of the log is independent variable, and discriminant function is generated using Fisher method of discrimination;
The selection unit, for utilizing the Lithology Discrimination sample set, the Lithology Discrimination accuracy of the discriminant function is examined, The target-recognition function for meeting the accuracy requirement of setting Lithology Discrimination is chosen from the discriminant function.
7. the correction system of carbonate formation well logging composite columnar section according to claim 6, it is characterised in that described Judging unit includes judgment sub-unit, the first correction subelement and the second correction subelement;
The judgment sub-unit, for exceeding baseline value 50%, and other logs in the hole diameter value of the log When corresponding well section is mutated, judge that there occurs expanding for current well section;
The first correction subelement, for when the expanding well section thickness of generation is less than 20m, utilizing adjacent non-expanding well section The value of the log is corrected to it;
The second correction subelement, for when the expanding well section thickness of generation is more than 20m, rejecting and expanding well section occurring The value of lithology classification and corresponding log.
8. the correction system of carbonate formation well logging composite columnar section according to claim 6, it is characterised in that described Sample set, which establishes unit, includes structure subelement and the 3rd correction subelement;
The structure subelement, for the lithology classification for section carbonate formation of being cored according to, selection correspondingly 5 to 7 The log after correction, build the Lithology Discrimination sample set;
It is described 3rd correction subelement, for judge in the Lithology Discrimination sample set with the presence or absence of lithology it is corresponding described in The depth point of log is less than the sample point of 10, if in the presence of rejecting the sample point.
9. the correction system of carbonate formation well logging composite columnar section according to claim 6, it is characterised in that described Choosing unit includes the firstth computation subunit, the second computation subunit, chooses subelement and rejects subelement;
First computation subunit, for the Lithology Discrimination sample set to be updated into the discriminant function, divided by calculating Group variate-value, obtain the lithology classification of the Lithology Discrimination sample set;
Second computation subunit, for the lithology classification according to the obtained Lithology Discrimination sample set and the section of coring The lithology breakdown result of carbonate formation, calculate the Lithology Discrimination accuracy of presently described discriminant function;
The selection subelement, for when the Lithology Discrimination accuracy of the discriminant function is more than or equal to 90%, choosing institute It is target-recognition function to state recognition function;
The rejecting subelement, for when the Lithology Discrimination accuracy of the discriminant function is less than 90%, rejecting Lithology Discrimination The minimum lithology classification of accuracy.
10. the correction system of the carbonate formation well logging composite columnar section according to any one of claim 6 to 9, it is special Sign is, the log is the log that can reflect petrophysical property, including density log curve, interval transit time Log, Natural Gamma-ray Logging Curves, resistivity series log and neutron porosity log curve.
CN201510441364.8A 2015-07-24 2015-07-24 Method and system for correcting carbonate rock stratum logging comprehensive histogram Active CN105089662B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510441364.8A CN105089662B (en) 2015-07-24 2015-07-24 Method and system for correcting carbonate rock stratum logging comprehensive histogram

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510441364.8A CN105089662B (en) 2015-07-24 2015-07-24 Method and system for correcting carbonate rock stratum logging comprehensive histogram

Publications (2)

Publication Number Publication Date
CN105089662A CN105089662A (en) 2015-11-25
CN105089662B true CN105089662B (en) 2018-04-06

Family

ID=54570893

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510441364.8A Active CN105089662B (en) 2015-07-24 2015-07-24 Method and system for correcting carbonate rock stratum logging comprehensive histogram

Country Status (1)

Country Link
CN (1) CN105089662B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108072915A (en) * 2016-11-16 2018-05-25 中国石油天然气股份有限公司 Method for identifying carbonate rock particle beach phase
CN109388816A (en) * 2017-08-07 2019-02-26 中国石油化工股份有限公司 A kind of hierarchical identification method of complex lithology
CN109989746A (en) * 2017-12-29 2019-07-09 中国石油天然气股份有限公司 Method and device for evaluating carbonate reservoir
CN111335870B (en) * 2018-12-18 2023-04-25 中国石油天然气股份有限公司 Method and device for determining oil and gas potential
CN114086948B (en) * 2020-12-08 2023-11-28 中国石油天然气集团有限公司 Carbonate rock classification and identification method and device, computer equipment and storage medium
CN116953797B (en) * 2023-08-14 2024-07-02 成都理工大学 Lake-phase thin-layer sandstone reservoir prediction method based on typical high-frequency gyratory isochrone contrast

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2313668C1 (en) * 2006-03-10 2007-12-27 Открытое акционерное общество "Объединенная энергетическая группа "Петросервис" Method for horizontal well bore drilling on the base of geological investigations
CN102071928A (en) * 2009-11-25 2011-05-25 中国石油天然气股份有限公司 Three-dimensional space volcanic lithology identification method
CN102436013A (en) * 2011-07-19 2012-05-02 北京师范大学 Method for partitioning glutenite sedimentary period based on Fischer diagram
CN102606150A (en) * 2012-03-30 2012-07-25 中国海洋石油总公司 Method and system for identifying fragmental rock lithological characters
CN103485758A (en) * 2013-08-12 2014-01-01 中国石油天然气股份有限公司 Complex lithology interpretation method and device utilizing multiple data fitting inversion
CN103615236A (en) * 2013-11-26 2014-03-05 中国石油集团川庆钻探工程有限公司 Method for monitoring formation pressure in real time through remote logging information
CN103711482A (en) * 2014-01-06 2014-04-09 山东科技大学 Method for quickly identifying clastic rock reservoir based on two well log overlapped picture
CN103775075A (en) * 2014-01-26 2014-05-07 中国海洋石油总公司 Method for identifying lithology of whole well section
CN104463686A (en) * 2014-10-29 2015-03-25 中国石油集团川庆钻探工程有限公司 Method for identifying shale gas reservoir while drilling by using discriminant analysis method
CN104632206A (en) * 2015-02-05 2015-05-20 中国石油天然气股份有限公司 Method and device for determining petroleum logging standard layer

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2313668C1 (en) * 2006-03-10 2007-12-27 Открытое акционерное общество "Объединенная энергетическая группа "Петросервис" Method for horizontal well bore drilling on the base of geological investigations
CN102071928A (en) * 2009-11-25 2011-05-25 中国石油天然气股份有限公司 Three-dimensional space volcanic lithology identification method
CN102436013A (en) * 2011-07-19 2012-05-02 北京师范大学 Method for partitioning glutenite sedimentary period based on Fischer diagram
CN102606150A (en) * 2012-03-30 2012-07-25 中国海洋石油总公司 Method and system for identifying fragmental rock lithological characters
CN103485758A (en) * 2013-08-12 2014-01-01 中国石油天然气股份有限公司 Complex lithology interpretation method and device utilizing multiple data fitting inversion
CN103615236A (en) * 2013-11-26 2014-03-05 中国石油集团川庆钻探工程有限公司 Method for monitoring formation pressure in real time through remote logging information
CN103711482A (en) * 2014-01-06 2014-04-09 山东科技大学 Method for quickly identifying clastic rock reservoir based on two well log overlapped picture
CN103775075A (en) * 2014-01-26 2014-05-07 中国海洋石油总公司 Method for identifying lithology of whole well section
CN104463686A (en) * 2014-10-29 2015-03-25 中国石油集团川庆钻探工程有限公司 Method for identifying shale gas reservoir while drilling by using discriminant analysis method
CN104632206A (en) * 2015-02-05 2015-05-20 中国石油天然气股份有限公司 Method and device for determining petroleum logging standard layer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
常规、成像测井等资料砂砾岩岩性识别研究;张洪江;《中国优秀硕士学位论文全文数据库 基础科学辑》;20120115;摘要,第36页 *

Also Published As

Publication number Publication date
CN105089662A (en) 2015-11-25

Similar Documents

Publication Publication Date Title
CN105089662B (en) Method and system for correcting carbonate rock stratum logging comprehensive histogram
CN109061765B (en) Trap evaluation method for heterogeneous thin sandstone interbed reservoir
CN104612675B (en) A kind of carbonate formation Lithology while drilling method for quickly identifying
Kadkhodaie-Ilkhchi et al. A review of reservoir rock typing methods in carbonate reservoirs: Relation between geological, seismic, and reservoir rock types
CN110318745B (en) Particle size lithology logging evaluation method under deposition microphase constraint
Lia et al. Uncertainties in reservoir production forecasts
CN103670388A (en) Method for evaluating organic carbon content of shale
Chehrazi et al. Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies
CN103775075A (en) Method for identifying lithology of whole well section
CN107807410A (en) Reservoir geophysical sensitive parameter method for optimizing and system
Askari et al. A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir
CN117365458A (en) Rock logging identification method for mixed-accumulation shale reservoir
KR101893800B1 (en) Method of sedimentary environment interpretation through electrofacies construction
Wheatley et al. Sedimentology, diagenesis, and reservoir characterization of the Permian White Rim Sandstone, southern Utah: Implications for carbon capture and sequestration potential
Eftekhari et al. Electrical facies of the Asmari Formation in the Mansouri oilfield, an application of multi-resolution graph-based and artificial neural network clustering methods
Asgari et al. A fully integrated approach for the development of rock type characterization, in a middle east giant carbonate reservoir
CN110399649A (en) A kind of metamorphic rock quantization recognition methods based on diagenesis indicator element
CN114059999B (en) Gravity flow deposition cause logging identification method
CN116930023A (en) Fine interpretation method and device for dense sandstone phase-control classified porosity logging
Taheri et al. Quantifying fracture density within the Asmari reservoir: an integrated analysis of borehole images, cores, and mud loss data to assess fracture-induced effects on oil production in the Southwestern Iranian Region
CN112147713B (en) Shale total organic carbon content segmented prediction method
Masoud et al. Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya
Amanipoor Providing a subsurface reservoir quality maps in oil fields by geostatistical methods
CN113514879A (en) &#39;ear&#39; layer identification method
Li et al. Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China

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