CN106869915A - A kind of Horizontal Well spacer interlayers Forecasting Methodology and device - Google Patents

A kind of Horizontal Well spacer interlayers Forecasting Methodology and device Download PDF

Info

Publication number
CN106869915A
CN106869915A CN201710149370.5A CN201710149370A CN106869915A CN 106869915 A CN106869915 A CN 106869915A CN 201710149370 A CN201710149370 A CN 201710149370A CN 106869915 A CN106869915 A CN 106869915A
Authority
CN
China
Prior art keywords
well
interlayer
horizontal
lead
type
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
CN201710149370.5A
Other languages
Chinese (zh)
Other versions
CN106869915B (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 Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas 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 China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201710149370.5A priority Critical patent/CN106869915B/en
Publication of CN106869915A publication Critical patent/CN106869915A/en
Application granted granted Critical
Publication of CN106869915B publication Critical patent/CN106869915B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a kind of Horizontal Well spacer interlayers Forecasting Methodology and device, it is related to clastic rock oil gas development technology field.Method includes:Determine to lead a well every the type of interlayer, distributing position and occurrence according to core data, it is determined that leading level of a well every interlayer;Explained using Logging Curves and lead type of a well every interlayer, and determine to lead interface depth of a well every interlayer using dipmeter logging curve;It is determined that leading tendency and inclination angle of a well every interlayer;Determine scope and type of the horizontal well every interlayer, determine horizontal well every the position of interlayer relative level well well;A well will be led to be changed to the section of horizontal well path place every the tendency of interlayer and inclination angle;Determine the production feature of the water injection well and producing well in purpose survey region, the production feature of the water injection well in purpose survey region and producing well and varying level well are matched every the corresponding production feature of Distribution Pattern of Interlayer, determines the net horizontal section of water injection well in purpose survey region and producing well every interbed distribution.

Description

A kind of Horizontal Well spacer interlayers Forecasting Methodology and device
Technical field
The present invention relates to clastic rock oil gas development technology field, more particularly to a kind of Horizontal Well spacer interlayers Forecasting Methodology And device.
Background technology
The main clastic sedimentary rock such as including various sandstone, conglomerate, siltstone of clastic reservoir, is the most important storage of China Collection channel type, is also one of main type of reservoir of world's oil gas field.Currently for example, the oil such as grand celebration, triumph, Liaohe River of China Field is mainly clastic reservoir;In the world boolean's glycerine field of famous Kuwait, Russia Sa Moteluoer oil fields It is clastic reservoir.In oil and gas development collection, horizontal well application is relatively broad;Horizontal well is to be applied to clastic reservoir Advanced technology technique in the middle of recovering the oil, horizontal well refers to hole angle and reaches or close to 90 °, and well bore creeps into one along horizontal direction The well of measured length, well is long through the length of oil reservoir, the contact surface area of well and oil reservoir is considerably increased, so that the list of oil well Well production is high, and the speed of production of oil well is fast, reduces the production time.
During clastic rock oil and gas development, the anisotropism of clastic reservoir is enhanced every the distribution of interlayer, be shadow Ring profit and preserving intrastratal flow, control the key factor of remaining oil distribution.At present, it is main both at home and abroad to use individual well well logging recognition With reference between Geological Mode, and well phased Geologic modeling method characterize every interlayer the regularity of distribution, the former has certain subjectivity Property, quantitative relationship effect is poor;The latter then more relies on algorithm, the characteristics of with realize more.Simultaneously as horizontal well is surveyed Instrument is widely different to be disposed vertically in the middle of horizontal positioned, with straight well during well, and both logging environment differences cause to be difficult to Every the well log interpretation of interlayer, the research to Horizontal Well spacer interlayers Forecasting Methodology then more lacks horizontal well.It can be seen that, it is current to lack It is few accurately to explain every interlayer and the method for portraying its spatial distribution.
The content of the invention
Embodiments of the invention provide a kind of Horizontal Well spacer interlayers Forecasting Methodology and device, to solve currently to lack essence Really explain the problem every interlayer and the method for portraying its spatial distribution.
To reach above-mentioned purpose, the present invention is adopted the following technical scheme that:
A kind of Horizontal Well spacer interlayers Forecasting Methodology, including:
The core data in purpose survey region is obtained, and is determined to lead class of a well every interlayer according to the core data Type, distributing position and occurrence, and determine to lead level of a well every interlayer every the distributing position of interlayer and occurrence according to a well is led;
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, is explained using Logging Curves and is led eye Well every interlayer type, and using dipmeter logging curve determine lead interface depth of a well every interlayer;
According to type and interface depth of a well every interlayer is led, the corresponding window of a well long, step-length and exploration are led in constantly adjustment Angle so that determine that lead a well leads eye every the interface depth of interlayer with according to core data determination using dipmeter logging curve Well is corresponding every interlayer occurrence, it is determined that leading tendency and inclination angle of a well every interlayer;
Net horizontal section log response abnormal data is obtained, and horizontal well is determined according to net horizontal section log response abnormal data Every the scope and type of interlayer, according to horizontal well every interlayer depth resistivity data, the phase of shallow resistivity data and resistivity reservoir Determine horizontal well every the position of interlayer relative level well well to variable quantity;
A well will be led on the section where the tendency of interlayer and inclination angle are changed to horizontal well path, to complete to lead a well every folder Layer matching every interlayer with horizontal well;
Predefine the creation data of water injection well and producing well and horizontal well every the scope of interlayer, type and horizontal well every The relation of the position of interlayer relative level well well, determines varying level well every the corresponding production feature of Distribution Pattern of Interlayer;
The creation data of the water injection well and producing well in purpose survey region is obtained, the water filling in purpose survey region is determined The production feature of well and producing well, according to lead a well every interlayer and horizontal well every interlayer matching result, by purpose survey region The production feature of interior water injection well and producing well and varying level well are matched every the corresponding production feature of Distribution Pattern of Interlayer, Determine the net horizontal section of water injection well in purpose survey region and producing well every interbed distribution.
Specifically, a well of leading includes shale, calcareous and calcium shale every the type of interlayer;Lead level of a well every interlayer Including in interlayer, layer, between sand body and in single sand body.
Specifically, lower limit standard strategy is responded every interlayer according to the different type for pre-setting, using Logging Curves Type of a well every interlayer is led in explanation, and determines to lead interface depth of a well every interlayer using dipmeter logging curve, including:
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, using leading natural gal of a well every interlayer Horse logging response character value, density log response characteristic value and formation resistivity logging response characteristic value are explained and lead a well every folder The type of layer;Wherein, the different type for pre-setting includes different type every interlayer correspondence every interlayer response lower limit standard strategy Gamma ray log response characteristic value span, the span and formation resistivity of density log response characteristic value The span of logging response character value;
Flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading a well every the interface of interlayer Depth.
Specifically, according to leading type and interface depth of a well every interlayer, constantly adjustment leads that the corresponding window of a well is long, step-length And search angle so that determine that lead a well determines every the interface depth of interlayer with according to the core data using dipmeter logging curve Lead that a well is corresponding every interlayer occurrence, it is determined that lead tendency and inclination angle of a well every interlayer, including:
Stratigraphic dip polar plot is determined according to dipmeter logging curve;
The corresponding window of a well long, step-length and search angle are led in constantly adjustment so that determine to lead eye using stratigraphic dip polar plot Well every interlayer interface depth with determined according to the core data to lead a well corresponding every interlayer occurrence;
It is determined that leading azimuth frequency diagram and Schmidt polar plot of a well every interlayer, main peak side in the azimuth frequency diagram is counted To as leading tendency and inclination angle of a well every interlayer.
Specifically, a well will be led on the section where the tendency of interlayer and inclination angle are changed to horizontal well path, to complete to lead Eye well matching every interlayer with horizontal well every interlayer, including:
According to formula:A well will be led every the tendency of interlayer and inclined Angle is changed to section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well;Wherein, α is Apparent dip;θ is true dip angle;ω is the angle between line of dip and true line of dip.
Specifically, predefine the creation data of water injection well and producing well and horizontal well every the scope of interlayer, type and Horizontal well determines varying level well every the corresponding production of Distribution Pattern of Interlayer every the relation of the position of interlayer relative level well well Feature, including:
According to water injection well and the creation data of producing well, by perforation correspondence and effective response for water flood analysis, production curve feature Research, it is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous dynamic every the corresponding individual well production of interlayer spread State response characteristic.
A kind of Horizontal Well spacer interlayers prediction meanss, including:
Core analysis unit, for obtaining the core data in purpose survey region, and determines according to the core data A well is led every the type of interlayer, distributing position and occurrence, and determines to lead eye every the distributing position of interlayer and occurrence according to a well is led Level of the well every interlayer;
A well is led every sandwich type and interface depth determining unit, for being rung every interlayer according to the different type for pre-setting Lower limit standard strategy is answered, is explained using Logging Curves and is led type of a well every interlayer, and determined using dipmeter logging curve Lead interface depth of a well every interlayer;
Tendency and inclination angle determining unit of a well every interlayer are led, type and interface depth of a well every interlayer are led for basis The corresponding window of a well long, step-length and search angle are led in degree, constantly adjustment so that determine to lead a well every interlayer using dipmeter logging curve Interface depth a well is corresponding every interlayer occurrence with leading of being determined according to the core data, it is determined that leading a well inclining every interlayer To and inclination angle;
Horizontal well every interlayer information determination unit, for obtaining net horizontal section log response abnormal data, and according to level Well section log response abnormal data determines scope and type of the horizontal well every interlayer, according to horizontal well every interlayer depth resistivity number Relative variation according to, shallow resistivity data and resistivity reservoir determines horizontal well every the position of interlayer relative level well well;
Tendency and inclination angle converting unit, change to horizontal well path place for will lead a well every the tendency of interlayer and inclination angle On section, with complete lead a well every interlayer the matching every interlayer with horizontal well;
Production characteristics determining unit, for predefining the creation data and horizontal well of water injection well and producing well every interlayer Scope, type and horizontal well determine varying level well every interbed distribution every the relation of the position of interlayer relative level well well The corresponding production feature of pattern;
Net horizontal section is every interbed distribution predicting unit, the life for obtaining water injection well and producing well in purpose survey region Data are produced, the production feature of the water injection well and producing well in purpose survey region is determined, according to leading a well every interlayer and horizontal well Every the matching result of interlayer, by the production feature of the water injection well in purpose survey region and producing well and varying level well every interlayer The corresponding production feature of distribution pattern is matched, determine the net horizontal section of water injection well in purpose survey region and producing well every Interbed distribution.
Specifically, the well of leading in the core analysis unit includes shale, calcareous and calcium shale every the type of interlayer;Lead Eye well includes interlayer every the level of interlayer, in layer, between sand body and in single sand body.
Additionally, a well is led every sandwich type and interface depth determining unit, specifically for:
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, using leading natural gal of a well every interlayer Horse logging response character value, density log response characteristic value and formation resistivity logging response characteristic value are explained and lead a well every folder The type of layer;Wherein, the different type for pre-setting includes different type every interlayer correspondence every interlayer response lower limit standard strategy Gamma ray log response characteristic value span, the span and formation resistivity of density log response characteristic value The span of logging response character value;
Flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading a well every the interface of interlayer Depth.
Additionally, tendency and inclination angle determining unit of a well every interlayer are led, specifically for:
Stratigraphic dip polar plot is determined according to dipmeter logging curve;
The corresponding window of a well long, step-length and search angle are led in constantly adjustment so that determine to lead eye using stratigraphic dip polar plot Well every interlayer interface depth with determined according to the core data to lead a well corresponding every interlayer occurrence;
It is determined that leading azimuth frequency diagram and Schmidt polar plot of a well every interlayer, main peak side in the azimuth frequency diagram is counted To as leading tendency and inclination angle of a well every interlayer.
In addition, tendency and inclination angle converting unit, specifically for:
According to formula:A well will be led every the tendency of interlayer and inclined Angle is changed to section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well;Wherein, α is Apparent dip;θ is true dip angle;ω is the angle between line of dip and true line of dip.
In addition, production characteristics determining unit, specifically for:
According to water injection well and the creation data of producing well, by perforation correspondence and effective response for water flood analysis, production curve feature Research, it is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous dynamic every the corresponding individual well production of interlayer spread State response characteristic.
A kind of Horizontal Well spacer interlayers Forecasting Methodology provided in an embodiment of the present invention and device, obtain purpose research first Core data in region, and determine to lead a well every the type of interlayer, distributing position and occurrence, and root according to the core data Determine to lead level of a well every interlayer every the distributing position of interlayer and occurrence according to a well is led;According to the different type for pre-setting every Interlayer responds lower limit standard strategy, is explained using Logging Curves and leads type of a well every interlayer, and bent using dipmeter logging Line determines to lead interface depth of a well every interlayer;According to type and interface depth of a well every interlayer is led, a well is led in constantly adjustment Corresponding window is long, step-length and search angle so that determine to lead interface depth and basis of a well every interlayer using dipmeter logging curve What the core data determined lead, and a well is corresponding every interlayer occurrence, it is determined that leading tendency and inclination angle of a well every interlayer;Obtain water Horizontal well section log response abnormal data, and according to net horizontal section log response abnormal data determine horizontal well every interlayer scope and Type, the relative variation according to horizontal well every interlayer depth resistivity data, shallow resistivity data and resistivity reservoir determines water Horizontal well is every the position of interlayer relative level well well;A well will be led to be changed to horizontal well path place every the tendency of interlayer and inclination angle On section, with complete lead a well every interlayer the matching every interlayer with horizontal well;Predefine the production number of water injection well and producing well According to horizontal well every the scope of interlayer, type and horizontal well every the position of interlayer relative level well well relation, it is determined that not Same level well is every the corresponding production feature of Distribution Pattern of Interlayer;Obtain the production of the water injection well and producing well in purpose survey region Data, determine the production feature of the water injection well and producing well in purpose survey region, according to lead a well every interlayer and horizontal well every The matching result of interlayer, by the production feature of the water injection well in purpose survey region and producing well and varying level well every interlayer point The corresponding production feature of cloth pattern is matched, and determines the net horizontal section of water injection well in purpose survey region and producing well every folder Layer distribution.The present invention can be solved the problems, such as currently to lack and accurately explained every interlayer and the method for portraying its spatial distribution, Data utilizes unification, the problem of Forecasting Methodology simplification during prior art can be improved, between can overcoming Horizontal Well Every the low shortcoming of interbed distribution precision of prediction, the precision of Horizontal Well spacer interlayers forecast of distribution is improve.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also Other accompanying drawings are obtained with according to these accompanying drawings.
Fig. 1 is a kind of flow chart one of Horizontal Well spacer interlayers Forecasting Methodology provided in an embodiment of the present invention;
Fig. 2 is a kind of flowchart 2 of Horizontal Well spacer interlayers Forecasting Methodology provided in an embodiment of the present invention;
Fig. 3 is the HD4-52 wells in the embodiment of the present invention every interlayer dipmeter logging curve processing result map;
Fig. 4 is the HD112H wells horizontal segment in the embodiment of the present invention every interlayer interpretation results figure;
Fig. 5 leads a well and matches result map with Horizontal Well spacer interlayers for the HD112H wells in the embodiment of the present invention;
Fig. 6 is the HD112H water flooding curves and HD4-52H production curve figures in the embodiment of the present invention;
Fig. 7 is that the HD112-HD112H-HD4-52 wells in the embodiment of the present invention predict result map every interlayer;
Fig. 8 is a kind of structural representation of Horizontal Well spacer interlayers prediction meanss provided in an embodiment of the present invention.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.
Realize it is of the invention during, inventor find prior art development and prior art Horizontal Well Spacer interlayers research is as follows:
Existing level well well spacer interlayers Forecasting Methodology has followed the conventional method that clastic rock is predicted every interlayer, i.e., only with often Rule well log interpretation individual well is every interlayer, the continuous layered distribution of well spacer interlayers.But the well spacing in horizontal well production area is (generally higher than 500m) it is typically much deeper than the well spacing (generally less than 200m) of conventional straight well oil-producing region, the Forecasting Methodology and concept of this simplification The distribution pattern of change obviously cannot meet the actual demand in oil field.
Although having carried out the well spacer interlayers forecasting research of straight well both at home and abroad, achieve some impressive progresses and into Really.But also in the starting stage effective methods and techniques, and inclination angle are not formed still to the prediction of Horizontal Well spacer interlayers Log, horizontal wellbore logging curve, horizontal well creation data are not used also in the middle of the prediction of well spacer interlayers.How rock is utilized Many data prediction level well well spacer interlayers such as the heart, conventional logging, dipmeter logging, horizontal wellbore logging and oil-water well creation data point Cloth is the key problem for being badly in need of solving.
Prior art related to the present invention mainly includes every interlayer individual well well logging recognition combination Geological Mode and well Between two aspects of phased Geologic modeling:
Prior art one:With rock core-log well as foundation, to classifying every interlayer, the well logging recognition every interlayer is established Standard, and individual well longitudinal direction has been carried out every interlayer identification, while the essential characteristic and the origin cause of formation every interlayer are studied, for the well every interlayer Between forecast of distribution (Shu Qinglin petroleum journals, 2006;What fragrant perfume offshore oil, 2012;Han Rubing, waits oil explorations and opens Hair, 2014).By to prototype model of appearing, rock core and well-log information comprehensive analysis, based on sand body configuration, study in reservoir Portion is every sandwich type, formation mechenism and distribution pattern, form, rule using the different parfacies of layer interfaces analytic approach research every interlayer Mould and its stacked pattern, carry out forecast of distribution between well (Lynds R.AAPG Bulletin, 2006;Wei Zhang.Petroleum Science, 2008;Print forest, waits Petroleum finances, 2013;Li Haiyan, waits oil explorations With exploitation, 2015).
The shortcoming of prior art one:
In clastic rock every mainly conventional logging recognition methods is used in the identification of interlayer, i.e., added by core observation Log modeling explain individual well every interlayer, but due to some positions development it is smaller every thickness of interlayer, Logging Curves hang down Accuracy of identification is not reached to resolution ratio, the method is recognizing there is certain error every sandwich type, thickness.To appear, rock core With well logging based on study geology distribution pattern every interlayer, generally to every interlayer form, scale, the research precision of distribution pattern Not enough, the geology distribution pattern of summary has larger artificial subjectivity, and the prediction of well spacer interlayers has uncertainty higher.
Prior art two:Electrical, the physical property recognition mode every interlayer are set up, knowledge of the whole district every interlayer is carried out using the pattern Other and across comparison.On this basis, respectively with Decided modelling, the control of physical property lower limit, the prediction of the method such as condition simulation every Be distributed between the well of interlayer, set up every interlayer three-dimensional distribution model (Zhang Guoyi waits science and technology Leaders, 2010).Entered using gray theory Row individual well is recognized every interlayer, on this basis, with reference to made ground measurement, to carrying out quantitatively characterizing every interlayer.With fixed every interlayer Gauge mould is constraints, and using the phased method being combined with stochastic modeling, foundation can reflect be distributed spy every mezzanine space The three-dimensional geological model levied (Sun Tianjian waits Petroleum finances, 2014).
The shortcoming of prior art two:
Development polytype is every interlayer in clastic reservoir, and different type is big every the form of interlayer, scale difference, space Complex distribution is very big, it is necessary to substantial amounts of made ground investigation and dense well pattern data are set up every the Quantitative anatomy research difficulty of interlayer Every the quantitative relationship of interlayer.But material resource condition, hydrodynamic condition that made ground is deposited with Gu etc. has larger difference, directly will What made ground observation was obtained is used for there is larger uncertainty among ancient deposition every interlayer scale.Meanwhile, three-dimensional geological is built The method of mould, with multiple results, last selection has larger artificial subjectivity.During being predicted every interlayer, do not make With dipmeter logging curve, horizontal wellbore logging curve, horizontal well creation data, technical basis are more single, and lack to horizontal well Every interlayer identification and the effective ways of inter-wells predicting.
In order to overcome above-mentioned problem of the prior art, as shown in figure 1, the embodiment of the present invention provides a kind of Horizontal Well interval Interlayer Forecasting Methodology, including:
Step 101, the core data obtained in purpose survey region, and determined to lead a well every folder according to the core data The type of layer, distributing position and occurrence, and determine to lead a well every interlayer every the distributing position of interlayer and occurrence according to a well is led Level.
The different type that step 102, basis pre-set responds lower limit standard strategy every interlayer, using Logging Curves Type of a well every interlayer is led in explanation, and determines to lead interface depth of a well every interlayer using dipmeter logging curve.
Step 103, according to leading type and interface depth of a well every interlayer, constantly adjustment leads that the corresponding window of a well is long, step Long and search angle so that lead a well using the determination of dipmeter logging curve true with according to the core data every the interface depth of interlayer It is fixed to lead that a well is corresponding every interlayer occurrence, it is determined that leading tendency and inclination angle of a well every interlayer.
Step 104, acquisition net horizontal section log response abnormal data, and it is true according to net horizontal section log response abnormal data Scope and type of the horizontal well every interlayer are determined, according to horizontal well every interlayer depth resistivity data, shallow resistivity data and oil reservoir electricity The relative variation of resistance rate determines horizontal well every the position of interlayer relative level well well.
Step 105, a well will be led changed to section where horizontal well path every the tendency of interlayer and inclination angle, to complete to lead Eye well matching every interlayer with horizontal well every interlayer.
Step 106, the scope of creation data and horizontal well every interlayer for predefining water injection well and producing well, type and Horizontal well determines varying level well every the corresponding production of Distribution Pattern of Interlayer every the relation of the position of interlayer relative level well well Feature.
The creation data of step 107, the water injection well and producing well obtained in purpose survey region, determines purpose survey region The production feature of interior water injection well and producing well, according to lead a well every interlayer and horizontal well every interlayer matching result, by purpose The production feature of water injection well and producing well in survey region and varying level well are every the corresponding production feature of Distribution Pattern of Interlayer Matched, determine the net horizontal section of water injection well in purpose survey region and producing well every interbed distribution.
A kind of Horizontal Well spacer interlayers Forecasting Methodology provided in an embodiment of the present invention, obtains in purpose survey region first Core data, and determine to lead a well every the type of interlayer, distributing position and occurrence according to the core data, and according to leading eye Well determines to lead level of a well every interlayer every the distributing position of interlayer and occurrence;Rung every interlayer according to the different type for pre-setting Lower limit standard strategy is answered, is explained using Logging Curves and is led type of a well every interlayer, and determined using dipmeter logging curve Lead interface depth of a well every interlayer;According to type and interface depth of a well every interlayer is led, it is corresponding that a well is led in continuous adjustment Window is long, step-length and search angle so that using dipmeter logging curve determine to lead a well every interlayer interface depth with according to the rock Calculation leads that a well is corresponding every interlayer occurrence according to what is determined, it is determined that leading tendency and inclination angle of a well every interlayer;Obtain net horizontal section Log response abnormal data, and scope and type of the horizontal well every interlayer are determined according to net horizontal section log response abnormal data, According to horizontal well every the relative variation of interlayer depth resistivity data, shallow resistivity data and resistivity reservoir determine horizontal well every The position of interlayer relative level well well;Well section where the tendency of interlayer and inclination angle are changed to horizontal well path will be led On, with complete lead a well every interlayer the matching every interlayer with horizontal well;Predefine the creation data of water injection well and producing well with Horizontal well, every the relation of the position of interlayer relative level well well, determines different water every the scope of interlayer, type and horizontal well Horizontal well is every the corresponding production feature of Distribution Pattern of Interlayer;Obtain the production number of the water injection well and producing well in purpose survey region According to the production feature of the water injection well and producing well in purpose survey region being determined, according to leading a well every interlayer and horizontal well every folder The matching result of layer, by the production feature of the water injection well in purpose survey region and producing well and varying level well every interbed distribution The corresponding production feature of pattern is matched, and determines the net horizontal section of water injection well in purpose survey region and producing well every interlayer Distribution.The present invention can be solved the problems, such as currently to lack and accurately explained every interlayer and the method for portraying its spatial distribution, can With data in improving prior art using unification, the problem of Forecasting Methodology simplification, Horizontal Well can be overcome to be spaced The low shortcoming of interbed distribution precision of prediction, improves the precision of Horizontal Well spacer interlayers forecast of distribution.
In order that those skilled in the art is better understood by the present invention, a more detailed embodiment is set forth below,. It should be evident that drawings in the following description are only specific embodiment of the invention, described embodiment is only this hair Bright a part of embodiment, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art are not having The every other embodiment obtained under the premise of creative work is made, the scope of the present invention is belonged to.
As shown in Fig. 2 the embodiment of the present invention provides a kind of Horizontal Well spacer interlayers Forecasting Methodology, including:
Step 201, the core data obtained in purpose survey region, and determined to lead a well every folder according to the core data The type of layer, distributing position and occurrence, and determine to lead a well every interlayer every the distributing position of interlayer and occurrence according to a well is led Level.
Specifically, a well of leading includes shale, calcareous and calcium shale every the type of interlayer;Lead level of a well every interlayer Including in interlayer, layer, between sand body and in single sand body.
Herein, the principle of step 201 is directed to research area's geological condition, and different its lithology every interlayer are different, drop hydrochloric acid Reactions vary.
Step 202, lower limit standard strategy is responded every interlayer according to the different type that pre-sets, using leading a well every interlayer Gamma ray log response characteristic value, density log response characteristic value and formation resistivity logging response characteristic value explain and lead Type of the eye well every interlayer.
Wherein, the different type for pre-setting includes that different type is corresponding every interlayer every interlayer response lower limit standard strategy The span of gamma ray log response characteristic value, the span of density log response characteristic value and formation resistivity are surveyed The span of well response characteristic value.
Step 203, flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading a well every folder The interface depth of layer.
Logging Curves are obvious every interlayer response characteristic to different type.According to different type every interlayer log response, By the method for Core-Calibrated Logging, natural gamma GR, density DEN, formation resistivity RT of the class of rock core well three every interlayer are read Logging response character value, set up different type every interlayer log response lower limit standard.It is possible thereby to according to Logging Curves To recognize the type every interlayer, but most thin only have several centimeters every interlayer, and Logging Curves cannot be accurately judged every interlayer Upper and lower interface.Compared to the sampling interval of Logging Curves 0.125m, dipmeter logging curve is 0.003m, and longitudinal frame can Up to 1cm, at the top bottom interface of interlayer, there is obvious flex point (resistivity becomes suddenly big point) in dipmeter logging curve, by This can exactly recognize the upper and lower interface every interlayer.
Step 204, stratigraphic dip polar plot is determined according to dipmeter logging curve.
The corresponding window of a well long, step-length and search angle are led in step 205, constantly adjustment so that utilize stratigraphic dip polar plot It is determined that lead a well every interlayer interface depth with determined according to the core data to lead a well corresponding every interlayer occurrence.
Azimuth frequency diagram and Schmidt polar plot of a well every interlayer are led in step 206, determination, count the azimuth frequency diagram Middle main peak direction, as leading tendency and inclination angle of a well every interlayer.
Dipmeter logging curve can obtain stratigraphic dip polar plot by treatment, can reflect well subsurface formations and every Interlayer occurrence information.Large change is often there is in small depth bounds every the occurrence of interlayer, it is right using short window correlation long Than method, by constantly adjustment, preferred interpretation processing parameter can obtain the azimuth frequency diagram and Schmidt's vector every interlayer after treatment Figure, main peak direction is the tendency every interlayer and inclination angle in the azimuth frequency diagram of statistics.
Step 207, acquisition net horizontal section log response abnormal data, and it is true according to net horizontal section log response abnormal data Scope and type of the horizontal well every interlayer are determined, according to horizontal well every interlayer depth resistivity data, shallow resistivity data and oil reservoir electricity The relative variation of resistance rate determines horizontal well every the position of interlayer relative level well well.
The logging environment of horizontal well is different from straight well, but both have similar log response at the position developed every interlayer Feature.Oil reservoirs of horizontal well Duan Shen, shallow resistivity curve are essentially coincided, and sensing is surveyed in horizontal well well diverse location every interlayer The response difference of well depth, shallow resistivity.
Step 208, a well will be led changed to section where horizontal well path every the tendency of interlayer and inclination angle, to complete to lead Eye well matching every interlayer with horizontal well every interlayer.
It is in different directions inconsistent every the occurrence of interlayer, dipmeter logging pickup is truly inclined to and is inclined every interlayer Angle, in order to realize that lead a well is combined the distribution of pre- well logging spacer interlayers, it is necessary to straight well is changed every interlayer occurrence with horizontal well Onto section where horizontal well track, with reach straight well with horizontal well matching every interlayer occurrence.
Herein, specifically can be according to formula (1):A well will be led On the section where the tendency of interlayer and inclination angle are changed to horizontal well path, to complete to lead a well every interlayer and horizontal well every interlayer Matching;Wherein, α is apparent dip;θ is true dip angle;ω is the angle between line of dip and true line of dip.
Extend trend and horizontal well every interlayer interpretation results every interlayer occurrence according to a well is led, carry out continuous every between interlayer well Journal of Sex Research.
Step 209, the creation data according to water injection well and producing well, are analyzed by perforation correspondence and effective response for water flood, produced Curvilinear characteristic is studied, it is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous every the corresponding list of interlayer spread Well Production development response characteristic.
It should be noted that difference influences different every interlayer areal extent on Production Wells situation, by perforation correspondence With effective response for water flood analysis, production curve properties study, continuous, semicontinuous, discontinuous three class can be summarized every interlayer areal extent Corresponding individual well Production development response characteristic.It is possible thereby to use the anti-distribution characteristics pushed away every interlayer of oil-water well creation data.
The creation data of step 210, the water injection well and producing well obtained in purpose survey region, determines purpose survey region The production feature of interior water injection well and producing well, according to lead a well every interlayer and horizontal well every interlayer matching result, by purpose The production feature of water injection well and producing well in survey region and varying level well are every the corresponding production feature of Distribution Pattern of Interlayer Matched, determine the net horizontal section of water injection well in purpose survey region and producing well every interbed distribution.
Herein, it is right that the individual well Production development response characteristic enumerated in the production feature of analysis individual well, with step 209 is carried out Than, can find the production feature classification for meeting the most, the anti-continuity pushed away every interbed distribution, comprehensively lead a well and horizontal well every The result of interlayer matching, completes Horizontal Well spacer interlayers forecast of distribution.
In order that above-mentioned steps 201 are clearer to step 210, a specific example is set forth below:
By above-mentioned steps 201 to step 210 to the " straight well-level of Tarim Basin Hadexun oilfield Donghe Sandstone The short section HD112-HD112H-HD4-52 of well-straight well " carries out Horizontal Well spacer interlayers prediction, and wherein HD112H is water filling Well, HD112 is a well of leading of HD112H, HD4-52H is producing well, and HD4-52 is that HD4-52H leads a well.
1. Tarim Basin Hadexun oilfield Donghe Sandstone HD112 well 5089.69-5104.05m well sections be have chosen as quilt Process object, the well section is bored and meets 5 substratums, wherein 1 substratum is bored meets 1 sand body, 2 substratums are bored meets 2 sand bodies;While also HD4- 52 well 5081.5-5091.17m well sections, the well section is bored and meets 4 substratums, wherein 1 substratum is bored meets 1 sand body, 2 substratums are bored meets 2 sand Body.
For research area's geological condition, according to core observation, every interbed distribution position and occurrence, the class every interlayer is determined Type and level, research separate interlayer and there are shale, calcareous and calcium shale three types, and shale is celadon mud every interlayer lithology Rock, argillaceous siltstoue, shale packsand, fill a large amount of shales between rock particles;Calcareous is that canescence is calcareous every interlayer lithology Siltstone, calcareous packsand, its calcareous content are high, and part drop hydrochloric acid reaction is violent;Calcium shale is by calcareous glue every interlayer lithology The mud stone of knot, siltstone containing mud and packsand, rock arrangement are fine and close.Research separates interlayer to be included between interlayer, sand body and single sand body Interior three levels.Layer spacer interlayers are present in substratum interface and occurrence is closely parallel;Sand body spacer interlayers are present in substratum in layer The interface that interior single sand body is stacked with single sand body, occurrence is inclined;In single sand body single sand body inside and occurrence inclination are located at every interlayer.
2. by the method for Core-Calibrated Logging, different type is set up every interlayer log response lower limit standard.Shale is every folder Layer is presented natural gamma numerical value (being more than 38API) higher, density value (2.38~2.58g/ higher on Logging Curves Cm3), the feature of resistivity value of a relatively high (being more than 0.6 Ω m);It is calcareous that nature is presented on Logging Curves every interlayer Gamma value relatively low (being less than 38API), density value (being more than 2.48g/cm3) high, the feature of resistivity value (being more than 0.9 Ω m) high; Calcium shale is presented natural gamma value relatively low (less than 38API) every interlayer on Logging Curves, density value is relatively low (is less than 2.48g/cm3), the feature of resistivity (being more than 0.7 Ω m) higher.Using natural gamma GR, density DEN, formation resistivity RT Logging response character value and dipmeter logging curve in the flex point at the bottom interface of interlayer top, explain H112 and HD4-52 every interlayer Interface depth and type (as shown in following table 1HD112 and HD4-52 every interlayer interpretation results table).
Table 1:
3. by constantly adjustment, final preferred interpretation process parameter can, 0.02m step-lengths long for 0.2m windows, 30 degree of explorations Angle, treatment HD112 and the dipmeter logging curve of the well section of HD4-52 mesh, count the azimuth frequency diagram and Schmidt's vector every interlayer The main peak direction of figure determines the occurrence (table 1) every interlayer, and Fig. 3 is HD4-52 individual wells every interlayer dipmeter logging curve processing achievement Figure, the wherein less possibility continuously distributed every between interlayer well in inclination angle is larger, and larger continuously distributed every between interlayer well in inclination angle can Energy property is smaller, but specific continuity degree will also be studied with reference to horizontal well every interlayer.
4. Fig. 4 is HD112H wells horizontal segment every interlayer interpretation results figure.Shale is anti-on horizontal wellbore logging curve every interlayer Reflect significantly increase for gamma ray curve numerical value, interval transit time curve numerical value it is relatively low;It is calcareous to be then reflected as natural gamma every interlayer Curve numerical value is substantially reduced, interval transit time curve numerical value is low;Calcium shale is reflected as natural gamma between shale every interlayer every interlayer And it is calcareous between interlayer, interval transit time curve numerical value it is higher.When depth, shallow resistivity have increasing by a relatively large margin relative to oil reservoir When big and numerical value is close, developed at well every interlayer;When depth, shallow resistivity is with respect to oil reservoir increase and deep resistivity is more than shallow electricity During resistance rate, more nearby developed in well every interlayer;When deep resistivity is with respect to the increase of resistivity reservoir amplitude, but shallow resistivity phase When being basically unchanged to oil reservoir, remotely developed in well every interlayer.
According to horizontal segment log, HD112H well 5220-5260m natural gammas and resistivity curve have obvious exception Value, reflects that it bores the breccia section for meeting Donghe Sandstone top, and horizontal well enters the Donghe Sandstone of overall more homogeneous from 5260m Section, the exception on natural gamma and sound wave curve can reflect type and the position every interlayer, while deep resistivity and oil reservoir electricity Resistance rate relative value is larger, is responded on shallow resistivity not substantially, illustrates to be located at well remote position every interlayer, thus may be used A shale in the range of well 0.5m is picked up at 5260~5385m every interlayer, one is picked up at 5275~5310m Apart from the calcareous every interlayer of well 0.5m or so, 5385~5615m pick up one outside well 0.5m scopes it is calcareous every Interlayer, is connected with the shale of 5260~5385m every interlayer, and one is identified at 5440~5615m apart from well 0.5m scopes Outer shale is every interlayer.Due to the detection feature of horizontal wellbore logging curve, be only capable of recognizing within investigative range every interlayer, other Position is not developed not necessarily, and be above investigative range and fail identification, this be accomplished by and lead a well carry out every interlayer With studying.
5. Fig. 5 leads a well and matches result map with Horizontal Well spacer interlayers for HD112H wells.It is small that horizontal well HD112H bores chance 1 The 1-3 sand bodies of layer and 2-2 sand bodies, the 2-3 sand bodies of 2 substratums, a corresponding well HD112 projection orientation of leading is 22 °, using formula (1) occurrence conversion is carried out, first shale is 22 ° of 0.5 ° of ∠ (tendency ∠ inclination angles) every the occurrence that interlayer projects to horizontal segment, the Two shales every interlayer be 22 ° of 0.34 ° of ∠, the 3rd it is calcareous every interlayer be 22 ° of 0.26 ° of ∠, these three every the inclination angle of interlayer all Very little.Matched every interlayer by with net horizontal section, mud of first shale every 5440~5615m of interlayer and HD112H well sections Matter is connected every interlayer, and second shale sand body spacer interlayers is connected with the shale of 5260~5385m of HD112H well sections every interlayer, the Three calcareous calcareous being connected every interlayer every interlayer and 5275~5310m of HD112H well sections, but due to the part after 5310m Not in investigative range, do not identify corresponding every interlayer, but do not represent it and do not exist, this is accomplished by dynamic data to enter One step research.
6. Fig. 6 is HD112H water flooding curves and HD4-52H production curve figures.In research area, perforation correspondence continuous every interlayer When effective response for water flood, the speed that takes effect of tracer, now producing well be presented that water-free oil production period is short, after water breakthrough moisture content it is quick on Rise, produce the feature that liquid rises daily;Do not taken effect when perforation is not to correspondence, now producing well presentation water-free oil production period is more long, day Produce the feature that liquid is gradually reduced.Semicontinuous every interlayer, no matter whether perforation corresponds to, and effective response for water flood can be presented on tracer test And the slow feature that takes effect, producing well is presented has certain water-free oil production period, and moisture content grows steadily after water breakthrough, daily output The feature of liquid stabilization.Discontinuous every interlayer, no matter whether perforation corresponds to, and effective response for water flood can be presented on tracer test and take effect The very slow feature of speed.
The HD112H wells metaideophone of in April, 2008 turns into water injection well, penetrates out 5339.42-5624.22m well sections, the i.e. 2-2 of 2 substratums With 2-3 sand bodies, HD4-52H penetrates out 2-2 the and 2-3 sand bodies of 5195-5498m well sections, i.e. 2 substratums in May, 2009, two mouthfuls of wells it Between perforation layer position correspondence, location of water injection is less than oil recovery position, the feature of " low note height adopt " is presented.HD112H was in November 26 in 2010 Day tracer is delivered, HD4-52H show that take effect speed in there is tracer after 276 days, and up to 3.84m/d, well operation is Water breakthrough, and the rising of initial stage moisture content is quickly, and 55% is risen to from May, 2009 to August moisture content from 27.9%, produces the liquid initial stage daily Also the feature of rising is presented.
7. Fig. 7 is that HD112-HD112H-HD4-52 wells predict result map every interlayer.Production feature meets " perforation above During to correspondence, the speed that takes effect (be more than 3.5m/d) of tracer, producing well is presented that water-free oil production period is short, after water breakthrough The feature that moisture content rapid increase, daily output liquid rise " illustrates that injection water is advanced by leaps and bounds along layer, and sand body between wells are connective preferably.Comprehensively lead Eye well and horizontal well every the result that interlayer match can be seen that interlayer and sand body spacer interlayers continuity is good, expanded range greatly, singly It is discontinuous every interlayer inside sand body, injection water can be effectively held, it is not flowed downward under gravity, can also hinder Only profit slows down the speed that takes effect in the mutual channelling of small interlayer.
HD112H wells are in 2 months 2009 injection rates from 170m3Drop to 120m3, later in August, 2010 (18 altogether Month) go back up to 150m3, stabilization is kept afterwards;On HD4-52H wells then it is corresponding show as from the moisture content of in November, 2009 from 50% drops to 37%, and in April, 2011 (altogether 17 months), moisture content rose to 54% later, stablizes left 60% afterwards The right side, both time intervals are 9 months, and the time consistency taken effect with tracer illustrates that profit is smooth in interflow and does not have Generation channelling, the production measure change of water injection well can be reflected on producing well, and this represents interlayer and sand body spacer interlayers are continuous Property is good, discontinuous every interlayer inside single sand body.Predicting the outcome before demonstrating, between illustrating the Horizontal Well of many Information integrations Every interlayer Forecasting Methodology high precision, production needs are disclosure satisfy that.
It can be seen that, by above-mentioned instantiation, this have the advantage that:
Method of the characterized systematically every interlayer static state geologic feature is creatively proposed, using rock core, conventional logging, inclination angle Well-log information carries out type every interlayer, level, individual well and explains and occurrence pickup, effectively improves static geologic feature to level The directive function of well well spacer interlayers forecast of distribution.Based on rock core information, differ according to every interlayer lithology, drip the anti-of hydrochloric acid Reflect differ, the difference of distributing position and thickness, carried out leading division of a well every sandwich type and level, solve existing method In there was only type and rare level every interlayer, it is impossible to for instructing the shortcoming every interbed distribution prediction;Dipmeter logging curve longitudinal direction Resolution ratio determines the method that type, dipmeter logging curve delimit face up to 1cm using Logging Curves, and explanation leads a well every folder The type of layer and position, solve in existing method because conventional logging resolution ratio is too low, it is impossible to accurately identify thickness compared with septulum The shortcoming of interlayer upper and lower interface;Dipmeter logging curve can obtain stratigraphic dip polar plot by treatment, can reflect well Every the occurrence of interlayer, by continuous preferred process parameter, pickup leads tendency and inclination angle of a well every interlayer, compensate for existing method In characterize not enough shortcoming every mezzanine space occurrence.The present invention improves the static geology sign precision every interlayer, disclosure satisfy that Production needs.
Creatively propose horizontal well to be recognized every interlayer and lead the method that a well is matched with horizontal well every interlayer, can make The type every interlayer and position with horizontal wellbore logging curve Direct Recognition, effectively compensate for horizontal wellbore logging data and are not used for level Section is every interlayer identification and the deficiency of Horizontal Well spacer interlayers forecast of distribution.Horizontal well with lead a well at the position developed every interlayer With similar logging response character, it is possible to use natural gamma and sound wave curve response abnormality carry out interpretation level well every interlayer Scope and type, solve the shortcoming that horizontal well in existing method cannot carry out being explained every interlayer;Every interlayer in horizontal well well Diverse location is deep to induction logging, shallow resistivity response is different, using the phase every interlayer deep and shallow resistivity and resistivity reservoir To variable quantity identification every position of the interlayer with respect to well, method of the horizontal well every mezzanine space location recognition is established;By leading Eye well is changed every the occurrence of interlayer, and trend and horizontal well are extended every interlayer interpretation results every interlayer occurrence based on a well is led, and is carried out The continuous Journal of Sex Research every between interlayer well, it is proposed that lead the method that a well is matched with horizontal well every interlayer.The present invention proposes a set of water The method that horizontal well is explained every interlayer, the need for disclosure satisfy that production.
Creatively propose the dynamic characterization method every interlayer, using oil-water well creation data it is counter push away every interlayer point Cloth, improves the confidence level of Horizontal Well spacer interlayers prediction.The situation of taking effect is adopted using the corresponding relation and note of perforation layer position, is ground Study carefully the spread of water injection well and producing well well spacer interlayers, compensate for note in existing method and adopt correspondence being only used for sand body company with taking effect property The general character is evaluated, and is not used for the deficiency of well spacer interlayers forecast of distribution;Summarize moisture content, daily output liquid curvilinear motion and every interlayer exhibition The relation of cloth, establishes the different methods for producing feature every the corresponding oil-water well of Distribution Pattern of Interlayer, in compensate for existing method Production curve is difficult to use in the deficiency of prediction level well well spacer interlayers distribution.
The present invention proposes a set of pre- every the dynamic characterization thinking of interlayer and the Horizontal Well spacer interlayers of many Information integrations Survey method, can all kinds of data of integrated application to greatest extent, reduce the uncertainty in the middle of interlayer inter-wells predicting so that pre- The result of survey is more reasonable, can be used in the middle of production practices.
Corresponding to the embodiment of the method described in above-mentioned Fig. 1 and Fig. 2, as shown in figure 8, the embodiment of the present invention provides a kind of level Well well spacer interlayers prediction meanss, including:
Core analysis unit 301 is for obtaining the core data in purpose survey region and true according to the core data Surely a well is led every the type of interlayer, distributing position and occurrence, and determines to lead every the distributing position of interlayer and occurrence according to a well is led Level of the eye well every interlayer.
A well is led every sandwich type and interface depth determining unit 302, for according to the different type for pre-setting every folder Layer response lower limit standard strategy, is explained using Logging Curves and leads type of a well every interlayer, and utilize dipmeter logging curve It is determined that leading interface depth of a well every interlayer.
Tendency and inclination angle determining unit 303 of a well every interlayer are led, type and interface of a well every interlayer are led for basis The corresponding window of a well long, step-length and search angle are led in depth, constantly adjustment so that determine to lead a well every folder using dipmeter logging curve The interface depth of layer with determined according to the core data to lead a well corresponding every interlayer occurrence, it is determined that leading a well every interlayer Tendency and inclination angle.
Horizontal well every interlayer information determination unit 304, for obtaining net horizontal section log response abnormal data, and according to water Horizontal well section log response abnormal data determines scope and type of the horizontal well every interlayer, according to horizontal well every interlayer depth resistivity number Relative variation according to, shallow resistivity data and resistivity reservoir determines horizontal well every the position of interlayer relative level well well.
Tendency and inclination angle converting unit 305, change to horizontal well path for will lead a well every the tendency of interlayer and inclination angle On the section of place, with complete lead a well every interlayer the matching every interlayer with horizontal well.
Production characteristics determining unit 306, for predefining the creation data and horizontal well of water injection well and producing well every folder Scope, type and the horizontal well of layer determine varying level well every interlayer every the relation of the position of interlayer relative level well well The corresponding production feature of distribution pattern.
Net horizontal section every interbed distribution predicting unit 307, for obtaining water injection well and producing well in purpose survey region Creation data, the production feature of the water injection well and producing well in purpose survey region is determined, according to leading a well every interlayer and water Horizontal well every interlayer matching result, by the production feature of the water injection well in purpose survey region and producing well and varying level well every The corresponding production feature of Distribution Pattern of Interlayer is matched, and determines the horizontal well of the water injection well and producing well in purpose survey region Section is every interbed distribution.
Specifically, the well of leading in the core analysis unit 301 includes shale, calcareous and calcium mud every the type of interlayer Matter;A well is led every the level of interlayer including in interlayer, layer, between sand body and in single sand body.
Additionally, a well is led every sandwich type and interface depth determining unit 302, specifically for:
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, using leading natural gal of a well every interlayer Horse logging response character value, density log response characteristic value and formation resistivity logging response characteristic value are explained and lead a well every folder The type of layer;Wherein, the different type for pre-setting includes different type every interlayer correspondence every interlayer response lower limit standard strategy Gamma ray log response characteristic value span, the span and formation resistivity of density log response characteristic value The span of logging response character value.
Flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading a well every the interface of interlayer Depth.
Additionally, tendency and inclination angle determining unit 303 of a well every interlayer are led, specifically for:
Stratigraphic dip polar plot is determined according to dipmeter logging curve.
The corresponding window of a well long, step-length and search angle are led in constantly adjustment so that determine to lead eye using stratigraphic dip polar plot Well every interlayer interface depth with determined according to the core data to lead a well corresponding every interlayer occurrence.
It is determined that leading azimuth frequency diagram and Schmidt polar plot of a well every interlayer, main peak side in the azimuth frequency diagram is counted To as leading tendency and inclination angle of a well every interlayer.
In addition, tendency and inclination angle converting unit 305, specifically for:
According to formula:A well will be led every the tendency of interlayer and inclined Angle is changed to section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well;Wherein, α is Apparent dip;θ is true dip angle;ω is the angle between line of dip and true line of dip.
In addition, production characteristics determining unit 306, specifically for:
According to water injection well and the creation data of producing well, by perforation correspondence and effective response for water flood analysis, production curve feature Research, it is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous dynamic every the corresponding individual well production of interlayer spread State response characteristic.
A kind of specific implementation of Horizontal Well spacer interlayers prediction meanss provided in an embodiment of the present invention may refer to The corresponding embodiments of the method for above-mentioned Fig. 1 and Fig. 2, here is omitted.
A kind of Horizontal Well spacer interlayers prediction meanss provided in an embodiment of the present invention, obtain in purpose survey region first Core data, and determine to lead a well every the type of interlayer, distributing position and occurrence according to the core data, and according to leading eye Well determines to lead level of a well every interlayer every the distributing position of interlayer and occurrence;Rung every interlayer according to the different type for pre-setting Lower limit standard strategy is answered, is explained using Logging Curves and is led type of a well every interlayer, and determined using dipmeter logging curve Lead interface depth of a well every interlayer;According to type and interface depth of a well every interlayer is led, it is corresponding that a well is led in continuous adjustment Window is long, step-length and search angle so that using dipmeter logging curve determine to lead a well every interlayer interface depth with according to the rock Calculation leads that a well is corresponding every interlayer occurrence according to what is determined, it is determined that leading tendency and inclination angle of a well every interlayer;Obtain net horizontal section Log response abnormal data, and scope and type of the horizontal well every interlayer are determined according to net horizontal section log response abnormal data, According to horizontal well every the relative variation of interlayer depth resistivity data, shallow resistivity data and resistivity reservoir determine horizontal well every The position of interlayer relative level well well;Well section where the tendency of interlayer and inclination angle are changed to horizontal well path will be led On, with complete lead a well every interlayer the matching every interlayer with horizontal well;Predefine the creation data of water injection well and producing well with Horizontal well, every the relation of the position of interlayer relative level well well, determines different water every the scope of interlayer, type and horizontal well Horizontal well is every the corresponding production feature of Distribution Pattern of Interlayer;Obtain the production number of the water injection well and producing well in purpose survey region According to the production feature of the water injection well and producing well in purpose survey region being determined, according to leading a well every interlayer and horizontal well every folder The matching result of layer, by the production feature of the water injection well in purpose survey region and producing well and varying level well every interbed distribution The corresponding production feature of pattern is matched, and determines the net horizontal section of water injection well in purpose survey region and producing well every interlayer Distribution.The present invention can be solved the problems, such as currently to lack and accurately explained every interlayer and the method for portraying its spatial distribution, can With data in improving prior art using unification, the problem of Forecasting Methodology simplification, Horizontal Well can be overcome to be spaced The low shortcoming of interbed distribution precision of prediction, improves the precision of Horizontal Well spacer interlayers forecast of distribution.
It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Apply the form of example.And, the present invention can be used and wherein include the computer of computer usable program code at one or more The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) is produced The form of product.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product Figure and/or block diagram are described.It should be understood that every first-class during flow chart and/or block diagram can be realized by computer program instructions The combination of flow and/or square frame in journey and/or square frame and flow chart and/or block diagram.These computer programs can be provided The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce A raw machine so that produced for reality by the instruction of computer or the computing device of other programmable data processing devices The device of the function of being specified in present one flow of flow chart or multiple one square frame of flow and/or block diagram or multiple square frames.
These computer program instructions may be alternatively stored in can guide computer or other programmable data processing devices with spy In determining the computer-readable memory that mode works so that instruction of the storage in the computer-readable memory is produced and include finger Make the manufacture of device, the command device realize in one flow of flow chart or multiple one square frame of flow and/or block diagram or The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented treatment, so as in computer or The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in individual square frame or multiple square frames.
Apply specific embodiment in the present invention to be set forth principle of the invention and implementation method, above example Explanation be only intended to help and understand the method for the present invention and its core concept;Simultaneously for those of ordinary skill in the art, According to thought of the invention, will change in specific embodiments and applications, in sum, in this specification Appearance should not be construed as limiting the invention.

Claims (12)

1. a kind of Horizontal Well spacer interlayers Forecasting Methodology, it is characterised in that including:
Obtain purpose survey region in core data, and according to the core data determine lead a well every interlayer type, point Cloth position and occurrence, and determine to lead level of a well every interlayer every the distributing position of interlayer and occurrence according to a well is led;
According to the different type for pre-setting every interlayer respond lower limit standard strategy, using Logging Curves explain lead a well every The type of interlayer, and determine to lead interface depth of a well every interlayer using dipmeter logging curve;
According to type and interface depth of a well every interlayer is led, the corresponding window of a well long, step-length and search angle are led in constantly adjustment, are made Must utilize dipmeter logging curve determine lead a well every interlayer interface depth with according to the core data determine lead a well every Interlayer occurrence is corresponding, it is determined that leading tendency and inclination angle of a well every interlayer;
Net horizontal section log response abnormal data is obtained, and determines horizontal well every folder according to net horizontal section log response abnormal data The scope and type of layer, according to horizontal well every interlayer depth resistivity data, the relative change of shallow resistivity data and resistivity reservoir Change amount determines horizontal well every the position of interlayer relative level well well;
A well will be led to be changed to section where horizontal well path every the tendency of interlayer and inclination angle, with complete to lead a well every interlayer with Matching of the horizontal well every interlayer;
The creation data of water injection well and producing well and horizontal well are predefined every the scope of interlayer, type and horizontal well every interlayer The relation of the position of relative level well well, determines varying level well every the corresponding production feature of Distribution Pattern of Interlayer;
Obtain the creation data of the water injection well and producing well in purpose survey region, determine water injection well in purpose survey region and The production feature of producing well, according to lead a well every interlayer and horizontal well every interlayer matching result, by purpose survey region The production feature of water injection well and producing well and varying level well are matched every the corresponding production feature of Distribution Pattern of Interlayer, it is determined that The net horizontal section of water injection well and producing well in purpose survey region is every interbed distribution.
2. Horizontal Well spacer interlayers Forecasting Methodology according to claim 1 a, it is characterised in that well of leading is every interlayer Type include shale, calcareous and calcium shale;A well is led every the level of interlayer including in interlayer, layer, between sand body and in single sand body.
3. Horizontal Well spacer interlayers Forecasting Methodology according to claim 2, it is characterised in that according to pre-setting not Same type responds lower limit standard strategy every interlayer, is explained using Logging Curves and leads type of a well every interlayer, and utilization is inclined Angle log determines to lead interface depth of a well every interlayer, including:
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, is surveyed using natural gamma of a well every interlayer is led Well response characteristic value, density log response characteristic value and formation resistivity logging response characteristic value are explained and lead a well every interlayer Type;Wherein, the different type for pre-setting includes that different type is corresponding certainly every interlayer every interlayer response lower limit standard strategy The span of gemma ray logging response characteristic value, the span of density log response characteristic value and formation resistivity logging The span of response characteristic value;
Flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading interface depth of a well every interlayer Degree.
4. Horizontal Well spacer interlayers Forecasting Methodology according to claim 3, it is characterised in that according to leading a well every interlayer Type and interface depth, constantly adjustment lead the corresponding window of a well long, step-length and search angle so that using dipmeter logging curve it is true Surely lead a well every interlayer interface depth with determined according to the core data to lead a well corresponding every interlayer occurrence, it is determined that leading Tendency and inclination angle of the eye well every interlayer, including:
Stratigraphic dip polar plot is determined according to dipmeter logging curve;
Constantly adjustment lead the corresponding window of a well long, step-length and search angle so that using stratigraphic dip polar plot determine lead a well every The interface depth of interlayer with determined according to the core data to lead a well corresponding every interlayer occurrence;
It is determined that leading azimuth frequency diagram and Schmidt polar plot of a well every interlayer, main peak direction in the azimuth frequency diagram is counted, As leading tendency and inclination angle of a well every interlayer.
5. Horizontal Well spacer interlayers Forecasting Methodology according to claim 4 a, it is characterised in that well will be led every interlayer Tendency and inclination angle are changed to section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well, Including:
According to formula:A well will be led to turn every the tendency of interlayer and inclination angle Shift on section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well;Wherein, α is depending on inclining Angle;θ is true dip angle;ω is the angle between line of dip and true line of dip.
6. Horizontal Well spacer interlayers Forecasting Methodology according to claim 5, it is characterised in that predefine water injection well and The creation data of producing well and horizontal well are every the scope of interlayer, type and horizontal well every the position of interlayer relative level well well Relation, determine varying level well every the corresponding production feature of Distribution Pattern of Interlayer, including:
According to water injection well and the creation data of producing well, by perforation correspondence and effective response for water flood analysis, production curve properties study, It is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous every the corresponding individual well Production development response of interlayer spread Feature.
7. a kind of Horizontal Well spacer interlayers prediction meanss, it is characterised in that including:
Core analysis unit, for obtaining the core data in purpose survey region, and determines to lead eye according to the core data Well every the type of interlayer, distributing position and occurrence, and according to lead a well every the distributing position of interlayer and occurrence determine to lead a well every The level of interlayer;
A well is led every sandwich type and interface depth determining unit, for according to the different type for pre-setting every under interlayer response Limit standard strategy, is explained using Logging Curves and leads type of a well every interlayer, and determines to lead eye using dipmeter logging curve Interface depth of the well every interlayer;
Tendency and inclination angle determining unit of a well every interlayer are led, type and interface depth of a well every interlayer is led for basis, no The corresponding window of a well long, step-length and search angle are led in disconnected adjustment so that determine to lead a well every the boundary of interlayer using dipmeter logging curve A well is corresponding every interlayer occurrence with leading of being determined according to the core data for face depth, it is determined that lead a well every the tendency of interlayer and Inclination angle;
Horizontal well every interlayer information determination unit, for obtaining net horizontal section log response abnormal data, and according to net horizontal section Log response abnormal data determines scope and type of the horizontal well every interlayer, according to horizontal well every interlayer depth resistivity data, shallow Resistivity data determines horizontal well every the position of interlayer relative level well well with the relative variation of resistivity reservoir;
Tendency and inclination angle converting unit, for well section where the tendency of interlayer and inclination angle are changed to horizontal well path will to be led On, with complete lead a well every interlayer the matching every interlayer with horizontal well;
Production characteristics determining unit, the model of creation data and horizontal well every interlayer for predefining water injection well and producing well Enclose, type and horizontal well every the relation of the position of interlayer relative level well well, determine varying level well every interbed distribution mould The corresponding production feature of formula;
Net horizontal section is every interbed distribution predicting unit, the production number for obtaining the water injection well in purpose survey region and producing well According to the production feature of the water injection well and producing well in purpose survey region being determined, according to leading a well every interlayer and horizontal well every folder The matching result of layer, by the production feature of the water injection well in purpose survey region and producing well and varying level well every interbed distribution The corresponding production feature of pattern is matched, and determines the net horizontal section of water injection well in purpose survey region and producing well every interlayer Distribution.
8. Horizontal Well spacer interlayers prediction meanss according to claim 7, it is characterised in that the core analysis unit In lead a well every interlayer type include shale, calcareous and calcium shale;Lead a well every the level of interlayer include interlayer, in layer, Between sand body and in single sand body.
9. Horizontal Well spacer interlayers prediction meanss according to claim 8, it is characterised in that lead a well every sandwich type And interface depth determining unit, specifically for:
Lower limit standard strategy is responded every interlayer according to the different type for pre-setting, is surveyed using natural gamma of a well every interlayer is led Well response characteristic value, density log response characteristic value and formation resistivity logging response characteristic value are explained and lead a well every interlayer Type;Wherein, the different type for pre-setting includes that different type is corresponding certainly every interlayer every interlayer response lower limit standard strategy The span of gemma ray logging response characteristic value, the span of density log response characteristic value and formation resistivity logging The span of response characteristic value;
Flex point of a well at the bottom interface of interlayer top is being led according to dipmeter logging curve, it is determined that leading interface depth of a well every interlayer Degree.
10. Horizontal Well spacer interlayers prediction meanss according to claim 9, it is characterised in that lead a well every interlayer Tendency and inclination angle determining unit, specifically for:
Stratigraphic dip polar plot is determined according to dipmeter logging curve;
Constantly adjustment lead the corresponding window of a well long, step-length and search angle so that using stratigraphic dip polar plot determine lead a well every The interface depth of interlayer with determined according to the core data to lead a well corresponding every interlayer occurrence;
It is determined that leading azimuth frequency diagram and Schmidt polar plot of a well every interlayer, main peak direction in the azimuth frequency diagram is counted, As leading tendency and inclination angle of a well every interlayer.
11. Horizontal Well spacer interlayers prediction meanss according to claim 10, it is characterised in that tendency and inclination angle conversion Unit, specifically for:
According to formula:A well will be led to turn every the tendency of interlayer and inclination angle Shift on section where horizontal well path, with complete to lead a well every interlayer the matching every interlayer with horizontal well;Wherein, α is depending on inclining Angle;θ is true dip angle;ω is the angle between line of dip and true line of dip.
12. Horizontal Well spacer interlayers prediction meanss according to claim 11, it is characterised in that production feature determines single Unit, specifically for:
According to water injection well and the creation data of producing well, by perforation correspondence and effective response for water flood analysis, production curve properties study, It is determined that it is continuous every interlayer spread, it is semicontinuous every interlayer spread and discontinuous every the corresponding individual well Production development response of interlayer spread Feature.
CN201710149370.5A 2017-03-14 2017-03-14 Horizontal well interval interlayer prediction method and device Active CN106869915B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710149370.5A CN106869915B (en) 2017-03-14 2017-03-14 Horizontal well interval interlayer prediction method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710149370.5A CN106869915B (en) 2017-03-14 2017-03-14 Horizontal well interval interlayer prediction method and device

Publications (2)

Publication Number Publication Date
CN106869915A true CN106869915A (en) 2017-06-20
CN106869915B CN106869915B (en) 2020-05-08

Family

ID=59170938

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710149370.5A Active CN106869915B (en) 2017-03-14 2017-03-14 Horizontal well interval interlayer prediction method and device

Country Status (1)

Country Link
CN (1) CN106869915B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107589469A (en) * 2017-08-04 2018-01-16 中国石油天然气股份有限公司 The determination method and apparatus of oil-water interfaces
CN110910267A (en) * 2019-11-18 2020-03-24 中国石油天然气股份有限公司 Method and device for determining marine carbonate rock interlayer
CN113933479A (en) * 2020-07-14 2022-01-14 中国石油天然气股份有限公司 Method and device for judging cause of mudstone interlayer in fluvial facies stratum

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102893283A (en) * 2010-05-21 2013-01-23 兰德马克绘图国际公司 Systems and methods for horizontal well correlation and geosteering
CN103410504A (en) * 2013-07-22 2013-11-27 中国石油天然气股份有限公司 Method and device for determining true resistivity of horizontal well/highly-deviated well
CN103573250A (en) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 Method for calculating distances between horizontal well borehole and upper and lower interfaces of stratum
CN104895554A (en) * 2015-04-30 2015-09-09 中国石油天然气股份有限公司 Depth and time conversion method and device of horizontal well logging data
CN105317431A (en) * 2014-07-30 2016-02-10 中国石油化工股份有限公司 Method for explaining and evaluating logging parameters of horizontal well
WO2016053330A1 (en) * 2014-10-01 2016-04-07 Halliburton Energy Services, Inc. Transposition of logs onto horizontal wells
CN105507890A (en) * 2014-09-26 2016-04-20 中国石油化工股份有限公司 Logging evaluation method suitable for shale gas reservoir geological model
CN105551083A (en) * 2015-12-08 2016-05-04 中国石油天然气集团公司 Method for depicting relationship between horizontal well track and strata
CN105804730A (en) * 2015-01-19 2016-07-27 国际商业机器公司 Resource identification method and system using historic well data
CN106204725A (en) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 A kind of horizontal well path bearing calibration and device

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102893283A (en) * 2010-05-21 2013-01-23 兰德马克绘图国际公司 Systems and methods for horizontal well correlation and geosteering
CN103410504A (en) * 2013-07-22 2013-11-27 中国石油天然气股份有限公司 Method and device for determining true resistivity of horizontal well/highly-deviated well
CN103573250A (en) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 Method for calculating distances between horizontal well borehole and upper and lower interfaces of stratum
CN105317431A (en) * 2014-07-30 2016-02-10 中国石油化工股份有限公司 Method for explaining and evaluating logging parameters of horizontal well
CN105507890A (en) * 2014-09-26 2016-04-20 中国石油化工股份有限公司 Logging evaluation method suitable for shale gas reservoir geological model
WO2016053330A1 (en) * 2014-10-01 2016-04-07 Halliburton Energy Services, Inc. Transposition of logs onto horizontal wells
CN105804730A (en) * 2015-01-19 2016-07-27 国际商业机器公司 Resource identification method and system using historic well data
CN104895554A (en) * 2015-04-30 2015-09-09 中国石油天然气股份有限公司 Depth and time conversion method and device of horizontal well logging data
CN105551083A (en) * 2015-12-08 2016-05-04 中国石油天然气集团公司 Method for depicting relationship between horizontal well track and strata
CN106204725A (en) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 A kind of horizontal well path bearing calibration and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
傅涛: "HDD油田4井区测井储层评价", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
李维禄等: "基于动、静模式约束的隔、夹层空间分布综合预测——以塔里木盆地东河塘油田井区为例", 《西北大学学报(自然科学版)》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107589469A (en) * 2017-08-04 2018-01-16 中国石油天然气股份有限公司 The determination method and apparatus of oil-water interfaces
CN110910267A (en) * 2019-11-18 2020-03-24 中国石油天然气股份有限公司 Method and device for determining marine carbonate rock interlayer
CN110910267B (en) * 2019-11-18 2023-11-28 中国石油天然气股份有限公司 Method and device for determining sea-phase carbonate rock interlayer
CN113933479A (en) * 2020-07-14 2022-01-14 中国石油天然气股份有限公司 Method and device for judging cause of mudstone interlayer in fluvial facies stratum
CN113933479B (en) * 2020-07-14 2023-09-26 中国石油天然气股份有限公司 Method and device for distinguishing factors of mudstone interlayer in river-phase stratum

Also Published As

Publication number Publication date
CN106869915B (en) 2020-05-08

Similar Documents

Publication Publication Date Title
CN107121699B (en) A kind of sedimentary facies identification method under earthquake phase control
CN106951660A (en) A kind of marine clastics horizontal well reservoir log interpretation method and device
CN109388817A (en) A kind of Reservoir Fracture three-dimensional modeling method
CN106501854B (en) Oil gas drilling target geological quantifying risk evaluation method
CN109061765A (en) The evaluation of trap method of heterogeneous thin sandstone alternating layers oil reservoir
CN107807407A (en) A kind of petroleum zone efficiency evaluation method and apparatus
CN106597544B (en) The brittle prediction technique of fine and close oil-gas reservoir and device
CN110517794A (en) A method of that establishes shale gas reservoir buries-thermal evolution history figure
CN104898161B (en) Effective sandstone predicting method based on logging response simulator
CN105093313B (en) A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device
CN103437760B (en) A kind of method using array induction data Fast Evaluation oil-water-layer
CN105510993A (en) Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method
CN105242307B (en) Carbonate complex seismic reservoir porosity acquisition methods and device
CN106894761A (en) Using the geosteering method of the geosteering model of time-domain seismic volume
CN109783944B (en) Three-dimensional flow unit model building method, device and equipment
CN107515957A (en) Mud shale SEQUENCE STRATIGRAPHIC DIVISION method
CN106597543A (en) Stratigraphic sedimentary facies division method
CN109870719A (en) A kind of carbonate rock thin tight reservoir distribution determination method, apparatus and system
CN106249300B (en) It shakes to combine based on well and determines source rock TOC contents and anisotropism method and apparatus
CN107861152A (en) A kind of Mesozoic Strata fine correlation method
CN106869915A (en) A kind of Horizontal Well spacer interlayers Forecasting Methodology and device
CN106940450A (en) Geosteering method for establishing model based on time-domain seismic volume
CN109469477A (en) A kind of prediction technique and device of man-made fracture extending direction
CN106908838A (en) The method for building target area stratum inelastic attenuation quality factor three-dimensional model
CN106990448B (en) A kind of marine clastics horizontal well is every interlayer log interpretation method and device

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