CN106869915B - Horizontal well interval interlayer prediction method and device - Google Patents

Horizontal well interval interlayer prediction method and device Download PDF

Info

Publication number
CN106869915B
CN106869915B CN201710149370.5A CN201710149370A CN106869915B CN 106869915 B CN106869915 B CN 106869915B CN 201710149370 A CN201710149370 A CN 201710149370A CN 106869915 B CN106869915 B CN 106869915B
Authority
CN
China
Prior art keywords
interlayer
well
pilot hole
horizontal well
production
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710149370.5A
Other languages
Chinese (zh)
Other versions
CN106869915A (en
Inventor
余义常
徐怀民
阳建平
徐朝晖
高兴军
宁超众
雷诚
韩如冰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to 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

Images

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 method and a device for predicting a horizontal well interval interlayer, and relates to the technical field of clastic rock oil and gas development. The method comprises the following steps: determining the type, distribution position and production state of the pilot hole well interlayer according to the core data, and determining the grade of the pilot hole well interlayer; explaining the type of the pilot hole well interlayer by using a conventional logging curve, and determining the interface depth of the pilot hole well interlayer by using an inclination logging curve; determining the inclination and dip angle of the pilot hole interlayer; determining the range and the type of a horizontal well interlayer, and determining the position of the horizontal well interlayer relative to a horizontal well borehole; converting the inclination and the inclination angle of the interlayer of the pilot hole well to the section where the horizontal well track is located; and determining the production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes, and determining the horizontal well section interlayer distribution of the water injection well and the oil production well in the target research area.

Description

Horizontal well interval interlayer prediction method and device
Technical Field
The invention relates to the technical field of oil and gas development of clastic rocks, in particular to a method and a device for predicting a horizontal well interval interlayer.
Background
The clastic rock reservoir mainly comprises clastic sedimentary rocks such as various sandstone, conglomerate, siltstone and the like, is the most important reservoir type in China, and is also one of the main reservoir types of the world oil and gas fields. At present, for example, oil fields in Daqing, Shengli, Liaohe and the like in China are mainly clastic rock reservoirs; the famous kowitt boolean oil field in the world, the russian samotero oil field, are also clastic rock reservoirs. In oil and gas development and acquisition, horizontal wells are widely applied; the horizontal well is an advanced technical process applied to oil extraction of a clastic rock reservoir, the horizontal well refers to a well with a well inclination angle reaching or approaching 90 degrees, a well body drills a well with a certain length along the horizontal direction, a well hole penetrates through an oil layer with a long length, and the contact surface area between the well and the oil layer is greatly increased, so that the single-well yield of the oil well is high, the production speed of the oil well is high, and the production time is shortened.
In the development process of clastic rock oil and gas, the distribution of the interlayer enhances the heterogeneity of clastic rock reservoirs, and is a key factor influencing the flow of oil and water in the reservoirs and controlling the distribution of residual oil. At present, single-well logging identification is mainly combined with a geological model at home and abroad, and an interwell phase-controlled geological modeling method is mainly adopted to represent the distribution rule of the interlayer, wherein the former has certain subjectivity and poor quantitative relation effect; the latter relies on the algorithm more, has the characteristics of realizing more. Meanwhile, because the horizontal well logging instruments are horizontally arranged and are greatly different from the vertical arrangement in a vertical well, the logging environments of the horizontal well logging instruments and the vertical well logging instruments are different, logging interpretation of the horizontal well interval interlayer is difficult, and research on a prediction method of the horizontal well interval interlayer is more lacking. It can be seen that there is currently a lack of a way to accurately interpret the interlayer and characterize its spatial distribution.
Disclosure of Invention
The embodiment of the invention provides a method and a device for predicting a horizontal well interval interlayer, which aim to solve the problem that a method for accurately explaining the interval interlayer and describing the spatial distribution of the interval interlayer is lacked at present.
In order to achieve the purpose, the invention adopts the following technical scheme:
a horizontal well interval interbed prediction method comprises the following steps:
obtaining core data in a target research area, determining the type, the distribution position and the production shape of a pilot hole well interlayer according to the core data, and determining the grade of the pilot hole well interlayer according to the distribution position and the production shape of the pilot hole well interlayer;
according to a preset lower limit standard strategy of response of different types of interlayer layers, the type of the guide eye well interlayer is explained by using a conventional logging curve, and the interface depth of the guide eye well interlayer is determined by using an inclination logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the rock core data;
acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the oil layer resistivity;
converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer;
the method comprises the steps of predetermining the relation between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes;
the method comprises the steps of obtaining production data of a water injection well and a production well in a target research area, determining production characteristics of the water injection well and the production well in the target research area, matching the production characteristics of the water injection well and the production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to matching results of a pilot hole interlayer and a horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the production well in the target research area.
Specifically, the types of the guide hole well interlayer comprise argillaceous substances, calcareous substances and calcium argillaceous substances; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
Specifically, according to a preset standard strategy of response lower limit of different types of isolation layers, the type of the pilot hole isolation layer is explained by using a conventional logging curve, and the interface depth of the pilot hole isolation layer is determined by using an inclination logging curve, and the method comprises the following steps:
according to preset lower limit standard strategies for response of different types of isolation interlayers, the types of the isolation interlayer of the pilot hole well are explained by using the natural gamma logging response characteristic value, the density logging response characteristic value and the formation resistivity logging response characteristic value of the isolation interlayer of the pilot hole well; the preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers;
and determining the interface depth of the pilot hole well interlayer according to the inflection point of the dip angle logging curve at the top and bottom interface of the pilot hole well interlayer.
Specifically, according to the type and the interface depth of the pilot hole interlayer, continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by using the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the core data, and the inclination angle of the pilot hole interlayer are determined, and the method comprises the following steps:
determining a formation dip angle vector diagram according to the dip angle logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the interface depth of the pilot hole interlayer determined by using the stratigraphic dip angle vector diagram corresponds to the pilot hole interlayer attitude determined according to the core data;
and determining an azimuth frequency diagram and a Schmidt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram as the inclination and dip angle of the pilot hole interlayer.
Specifically, the method comprises the following steps of converting the inclination and the inclination angle of the pilot hole interlayer to a section where a horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer, wherein the method comprises the following steps:
according to the formula:
Figure BDA0001245087870000031
and converting the inclination and the inclination angle of the pilot hole interlayer to the section where the horizontal well track is located to complete the matching of the pilot hole interlayer and the horizontal well interlayer, wherein α is an apparent inclination angle, theta is a true inclination angle, and omega is an included angle between an apparent inclination line and a true inclination line.
Specifically, the method comprises the steps of predetermining the relationship between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes, wherein the method comprises the following steps:
and according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
A horizontal well interval interbed prediction device, comprising:
the core analysis unit is used for acquiring core data in a target research area, determining the type, the distribution position and the production state of the pilot hole interlayer according to the core data, and determining the grade of the pilot hole interlayer according to the distribution position and the production state of the pilot hole interlayer;
the pilot hole interlayer type and interface depth determining unit is used for explaining the type of the pilot hole interlayer by using a conventional logging curve according to a preset lower limit standard strategy of response of different types of interlayers, and determining the interface depth of the pilot hole interlayer by using an inclination angle logging curve;
the inclination and inclination angle determining unit is used for continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination and inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the production state of the pilot hole interlayer determined according to the rock core data;
the horizontal well interlayer information determining unit is used for acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the resistivity of an oil layer;
the inclination and inclination angle conversion unit is used for converting the inclination and inclination angle of the pilot hole interlayer to a section where a horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer;
the production characteristic determining unit is used for determining the relation between the production data of the water injection well and the oil production well and the range and the type of the horizontal well interlayer and the position of the horizontal well interlayer relative to the horizontal well hole in advance and determining the production characteristics corresponding to the distribution modes of the horizontal well interlayers;
and the horizontal well section interlayer distribution prediction unit is used for acquiring the production data of the water injection well and the oil production well in the target research area, determining the production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to the matching result of the pilot hole interlayer and the horizontal well interlayer, and determining the horizontal well section interlayer distribution of the water injection well and the oil production well in the target research area.
Specifically, the types of the guide hole well interlayer in the core analysis unit comprise argillaceous substances, calcareous substances and calcium argillaceous substances; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
In addition, the pilot hole interlayer type and interface depth determining unit is specifically used for:
according to preset lower limit standard strategies for response of different types of isolation interlayers, the types of the isolation interlayer of the pilot hole well are explained by using the natural gamma logging response characteristic value, the density logging response characteristic value and the formation resistivity logging response characteristic value of the isolation interlayer of the pilot hole well; the preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers;
and determining the interface depth of the pilot hole well interlayer according to the inflection point of the dip angle logging curve at the top and bottom interface of the pilot hole well interlayer.
In addition, the inclination and dip angle determination unit of the pilot hole interlayer is specifically configured to:
determining a formation dip angle vector diagram according to the dip angle logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the interface depth of the pilot hole interlayer determined by using the stratigraphic dip angle vector diagram corresponds to the pilot hole interlayer attitude determined according to the core data;
and determining an azimuth frequency diagram and a Schmidt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram as the inclination and dip angle of the pilot hole interlayer.
In addition, the inclination and tilt conversion unit is specifically configured to:
according to the formula:
Figure BDA0001245087870000051
and converting the inclination and the inclination angle of the pilot hole interlayer to the section where the horizontal well track is located to complete the matching of the pilot hole interlayer and the horizontal well interlayer, wherein α is an apparent inclination angle, theta is a true inclination angle, and omega is an included angle between an apparent inclination line and a true inclination line.
In addition, the production characteristic determination unit is specifically configured to:
and according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
The method and the device for predicting the horizontal well interlayer provided by the embodiment of the invention comprise the steps of firstly obtaining core data in a target research area, determining the type, the distribution position and the production shape of the pilot hole interlayer according to the core data, and determining the level of the pilot hole interlayer according to the distribution position and the production shape of the pilot hole interlayer; according to a preset lower limit standard strategy of response of different types of interlayer layers, the type of the guide eye well interlayer is explained by using a conventional logging curve, and the interface depth of the guide eye well interlayer is determined by using an inclination logging curve; continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the rock core data; acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the oil layer resistivity; converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer; the method comprises the steps of predetermining the relation between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes; the method comprises the steps of obtaining production data of a water injection well and a production well in a target research area, determining production characteristics of the water injection well and the production well in the target research area, matching the production characteristics of the water injection well and the production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to matching results of a pilot hole interlayer and a horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the production well in the target research area. The method can solve the problem that the method for accurately explaining the interlayer and describing the spatial distribution of the interlayer is lacked at present, can improve the problems of single utilization of data and simplification of a prediction method in the prior art, can overcome the defect of low prediction precision of the distribution of the interlayer of the horizontal well, and improves the prediction precision of the distribution of the interlayer of the horizontal well.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a first flow chart of a horizontal well interval interbed prediction method provided by an embodiment of the invention;
fig. 2 is a second flowchart of a horizontal well interval interbed prediction method according to an embodiment of the present invention;
FIG. 3 is a diagram of the processing results of an HD4-52 well bay dip log in accordance with an embodiment of the present invention;
FIG. 4 is a diagram illustrating the results of an HD112H well horizontal interval layer in an embodiment of the present invention;
FIG. 5 is a diagram of the matching result of the interval interlayer between the HD112H pilot hole well and the horizontal well in the embodiment of the invention;
FIG. 6 is a graph showing the water injection curve of HD112H and the production curves of HD4-52H in the embodiment of the present invention;
FIG. 7 is a graph of predicted results for HD 112-HD 112H-HD 4-52 well spacers in accordance with an embodiment of the present invention;
fig. 8 is a schematic structural diagram of a horizontal well interval interbed prediction device provided in an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the process of implementing the invention, the inventor finds the development condition of the prior art and the interval interlayer research of the horizontal well of the prior art as follows:
the conventional method for predicting the clastic rock interlayer is adopted in the conventional method for predicting the horizontal well interlayer, namely, the continuous layered distribution of the single well interlayer and the well interlayer is explained only by using conventional well logging. However, the well spacing (generally greater than 500m) of the horizontal well production area is usually much larger than that (generally less than 200m) of the conventional vertical well production area, and the simplified prediction method and the conceptualized distribution mode obviously cannot meet the actual requirements of the oil field.
Although the prediction research of the well interval interlayer of the vertical well has been carried out at home and abroad, some important progresses and achievements are achieved. However, the interval prediction of the horizontal well is still in a starting stage, no effective method and technology is formed, and the dip angle well logging curve, the horizontal well logging curve and the horizontal well production data are not used in the interval prediction. How to predict the interval interlayer distribution of the horizontal well by using multiple data such as rock cores, conventional well logging, dip angle well logging, horizontal well logging, oil-water well production data and the like is a critical problem which needs to be solved urgently.
The prior art scheme related to the invention mainly comprises two aspects of interlayer single well logging identification combined with geological mode and interwell phased geological modeling:
the first prior art is as follows: the method comprises the steps of classifying the interlayer by taking a rock core-logging as a basis, establishing a logging identification standard of the interlayer, identifying the single-well longitudinal interlayer, researching the basic characteristics and the cause of the interlayer at the same time, and predicting the inter-well distribution of the interlayer (Duqinglin. Petroleum institute, 2006; joss essence. marine oil, 2012; Korea such as ice, and the like. oil exploration and development, 2014). Through comprehensive analysis of outcrop prototype models, cores and logging data, the types, formation mechanisms and distribution modes of interlayer layers in reservoirs are researched based on sand body configurations, and the forms, scales and stacking patterns of different sub-phase interlayer layers are researched by utilizing a hierarchical interface analysis method to predict well distribution (Lynds R.AAPG Bulletin, 2006; WeiZhang. petroleum Science, 2008; forest printing, et al. oil exploration and development, 2013; Lihai swallow, et al. oil exploration and development, 2015).
The first prior art has the following defects:
the method mainly adopts a conventional logging identification method on the identification of the clastic rock interlayer, namely, the single-well interlayer is explained by core observation and logging curve modeling, but because the thickness of the interlayer developed at certain positions is small, the vertical resolution of the conventional logging curve cannot reach the identification precision, and certain errors exist in the identification of the type and the thickness of the interlayer. The geological distribution mode of the interlayer is researched mainly by exposing, rock core and well logging, the research precision of the shape, scale and distribution pattern of the interlayer is not enough on the whole, the summarized geological distribution mode has large human subjectivity, and the prediction of the well interlayer has high uncertainty.
The second prior art is: and establishing an electrical property and physical property identification mode of the interlayer, and carrying out identification and transverse comparison of the whole interlayer by utilizing the mode. On the basis, the well distribution of the interlayer is predicted by using methods such as deterministic modeling, physical property lower limit value control, condition simulation and the like, and a three-dimensional interlayer spreading model is established (Zhangyi, et al. science and technology guidance, 2010). And (3) identifying the interlayer of the single well by adopting a grey theory, and quantitatively characterizing the interlayer by combining modern deposition measurement on the basis. A three-dimensional geological model (Sunzian, et al. oil exploration and development, 2014) capable of reflecting the spatial distribution characteristics of the interlayer is established by taking the quantitative scale of the interlayer as a constraint condition and adopting a method combining phase control and random modeling.
The second prior art has the following defects:
various types of interlayers are developed in the clastic rock reservoir, the different types of interlayers have large differences in form and scale, the space distribution is complex, the quantitative dissection research difficulty of the interlayers is high, and a large amount of modern deposition investigation and dense well pattern data are needed to establish the quantitative relationship of the interlayers. However, the source conditions, hydrodynamic conditions and the like of modern deposition and ancient deposition are greatly different, and great uncertainty exists when the interlayer scale obtained by modern deposition observation is directly used in the ancient deposition. Meanwhile, the three-dimensional geological modeling method has a plurality of results, and the final selection has larger artificial subjectivity. In the interlayer prediction process, an inclination angle well logging curve, a horizontal well logging curve and horizontal well production data are not used, the technical basis is single, and an effective method for identifying horizontal well interlayers and predicting wells is lacked.
In order to overcome the problems in the prior art, as shown in fig. 1, an embodiment of the present invention provides a method for predicting a horizontal well interval zone, including:
step 101, obtaining core data in a target research area, determining the type, distribution position and production shape of a pilot hole well interlayer according to the core data, and determining the grade of the pilot hole well interlayer according to the distribution position and the production shape of the pilot hole well interlayer.
And 102, according to preset response lower limit standard strategies of different types of isolation layers, explaining the type of the isolation layer of the pilot hole by using a conventional logging curve, and determining the interface depth of the isolation layer of the pilot hole by using an inclination logging curve.
And 103, continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the interface depth of the pilot hole interlayer determined by using the dip angle logging curve corresponds to the attitude of the pilot hole interlayer determined according to the core data, and determining the inclination and the dip angle of the pilot hole interlayer.
And 104, acquiring logging response abnormal data of the horizontal well section, determining the range and the type of the horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to the horizontal well borehole according to the deep resistivity data of the horizontal well interlayer, the shallow resistivity data and the relative variation of the oil layer resistivity.
And 105, converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located, so as to complete matching of the pilot hole interlayer and the horizontal well interlayer.
And 106, predetermining the relation between the production data of the water injection well and the oil production well and the range and type of the horizontal well interlayer and the position of the horizontal well interlayer relative to the horizontal well hole, and determining the production characteristics corresponding to the distribution modes of the horizontal well interlayers.
And 107, acquiring production data of the water injection well and the oil production well in the target research area, determining the production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to the matching result of the pilot hole interlayer and the horizontal well interlayer, and determining the horizontal well interval interlayer distribution of the water injection well and the oil production well in the target research area.
The horizontal well interlayer prediction method provided by the embodiment of the invention comprises the steps of firstly obtaining core data in a target research area, determining the type, the distribution position and the production shape of a pilot hole interlayer according to the core data, and determining the level of the pilot hole interlayer according to the distribution position and the production shape of the pilot hole interlayer; according to a preset lower limit standard strategy of response of different types of interlayer layers, the type of the guide eye well interlayer is explained by using a conventional logging curve, and the interface depth of the guide eye well interlayer is determined by using an inclination logging curve; continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the rock core data; acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the oil layer resistivity; converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer; the method comprises the steps of predetermining the relation between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes; the method comprises the steps of obtaining production data of a water injection well and a production well in a target research area, determining production characteristics of the water injection well and the production well in the target research area, matching the production characteristics of the water injection well and the production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to matching results of a pilot hole interlayer and a horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the production well in the target research area. The method can solve the problem that the method for accurately explaining the interlayer and describing the spatial distribution of the interlayer is lacked at present, can improve the problems of single utilization of data and simplification of a prediction method in the prior art, can overcome the defect of low prediction precision of the distribution of the interlayer of the horizontal well, and improves the prediction precision of the distribution of the interlayer of the horizontal well.
A more detailed example is set forth below to provide a better understanding of the present invention to those skilled in the art. It is to be understood that the drawings in the following description are merely exemplary embodiments of the invention, and that the embodiments described are only a few, and not all, embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the scope of the present invention.
As shown in fig. 2, an embodiment of the present invention provides a horizontal well interval interbed prediction method, including:
step 201, obtaining core data in a target research area, determining the type, distribution position and production state of a pilot hole well interlayer according to the core data, and determining the grade of the pilot hole well interlayer according to the distribution position and the production state of the pilot hole well interlayer.
Specifically, the types of the guide hole well interlayer comprise argillaceous substances, calcareous substances and calcium argillaceous substances; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
Here, the principle of step 201 is that for the geological conditions of the research area, the lithology of different interlayer is different, and the reaction of dripping hydrochloric acid is different.
Step 202, according to preset lower limit standard strategies of response of different types of isolation layers, the types of the isolation layers of the pilot hole well are explained by using the natural gamma logging response characteristic values, the density logging response characteristic values and the formation resistivity logging response characteristic values of the isolation layers of the pilot hole well.
The preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers.
And 203, determining the interface depth of the pilot hole interlayer according to the inflection point of the dip angle logging curve at the top and bottom interfaces of the pilot hole interlayer.
The response characteristics of the conventional well logging curve to different types of interlayer are obvious. According to logging responses of different types of interlayer, logging response characteristic values of natural gamma GR, density DEN and formation resistivity RT of three types of interlayer of the core well are read through a core scale logging method, and the lower limit standard of the logging response of the different types of interlayer is established. Therefore, the type of the interlayer can be identified according to the conventional well logging curve, but the thinnest interlayer is only a few centimeters, and the upper interface and the lower interface of the interlayer cannot be accurately judged according to the conventional well logging curve. Compared with a conventional logging curve with a sampling interval of 0.125m, the dip angle logging curve is 0.003m, the longitudinal resolution can reach 1cm, and an obvious inflection point (a point with suddenly increased resistivity) exists in the dip angle logging curve at the top and bottom interfaces of the interlayer, so that the upper and lower interfaces of the interlayer can be accurately identified.
And 204, determining a formation dip angle vector diagram according to the dip angle logging curve.
And step 205, continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole well, so that the interface depth of the pilot hole well interlayer determined by the stratigraphic dip angle vector diagram corresponds to the pilot hole well interlayer attitude determined according to the rock core data.
And step 206, determining an azimuth frequency diagram and a Schmitt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram to be used as the inclination and dip angle of the pilot hole interlayer.
The dip angle well logging curve can be processed to obtain a stratum dip angle vector diagram, and the information of the underground stratum and the attitude of the interlayer can be well reflected. The appearance of the interlayer is often greatly changed in a small depth range, a short-window long-correlation contrast method is adopted, processing parameters are preferably explained through continuous adjustment, an azimuth frequency diagram and a Schmidt vector diagram of the interlayer can be obtained after processing, and the main peak direction in the statistical azimuth frequency diagram is the inclination and dip angle of the interlayer.
And step 207, acquiring logging response abnormal data of the horizontal well section, determining the range and the type of the horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to the horizontal well borehole according to the deep resistivity data of the horizontal well interlayer, the shallow resistivity data and the relative variation of the oil layer resistivity.
The logging environment of a horizontal well is different from that of a vertical well, but the two wells have similar logging response characteristics at the position where the interlayer develops. The oil layer depth and shallow resistivity curves of the horizontal well are basically overlapped, and the response of the interlayer to the induction logging depth and shallow resistivity at different positions of the horizontal well is different.
And 208, converting the inclination and the inclination angle of the pilot hole interlayer to the section where the horizontal well track is located so as to complete the matching of the pilot hole interlayer and the horizontal well interlayer.
The attitude of the interlayer is inconsistent in different directions, the inclination angle logging picks up the real inclination and inclination angle of the interlayer, and in order to predict the distribution of the interlayer of the well by combining the pilot hole well and the horizontal well, the attitude of the interlayer of the vertical well needs to be converted to the section where the well track of the horizontal well is located, so that the attitude of the interlayer of the vertical well and the horizontal well is matched.
Here, the following formula (1) may be specifically used:
Figure BDA0001245087870000111
the inclination and the inclination angle of the pilot hole interlayer are converted to the section where the horizontal well track is located, so that the matching of the pilot hole interlayer and the horizontal well interlayer is completedThe formula is shown in the specification, wherein α is an apparent inclination angle, theta is a true inclination angle, and omega is an included angle between an apparent inclination line and a true inclination line.
And carrying out interlayer inter-well continuity research according to the attitude extension trend of the interlayer of the pilot hole and the interpretation result of the interlayer of the horizontal well.
And 209, according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
It should be noted that the production conditions of the oil production well are affected differently by the different interlayer spreading ranges, and the single-well production dynamic response characteristics corresponding to the continuous, semi-continuous and discontinuous interlayer spreading ranges can be summarized through perforation correspondence, water injection effect analysis and production curve characteristic research. The distribution characteristics of the interlayer can thus be inferred using oil and water well production data.
Step 210, obtaining production data of the water injection well and the oil production well in the target research area, determining production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to the matching result of the pilot hole interlayer and the horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the oil production well in the target research area.
Here, the production characteristics of the single well are analyzed and compared with the dynamic response characteristics of the single well production listed in step 209, so that the best-fit production characteristic category can be found, the continuity of the interlayer distribution is reversely deduced, and the interlayer matching result of the pilot hole well and the horizontal well is synthesized, so that the interlayer distribution prediction of the horizontal well is completed.
In order to make the above steps 201 to 210 more clear, a specific example is listed below:
horizontal well interval interlayer prediction is carried out on short sections HD 112-HD 112H-HD 4-52 of 'vertical well-horizontal well-vertical well' of Tarilwood basin Hadison east river sandstone, wherein HD112H is a water injection well, HD112 is a pilot hole of HD112H, HD4-52H is a production well, and HD4-52 is a pilot hole of HD4-52H through the steps 201 to 210.
① A well section of Turkey Hadison oil field Donghe sandstone HD112 well 5089.69-5104.05m is selected as the object to be processed, the well section encounters 5 small layers, wherein 1 small layer encounters 1 sand body, 2 small layers encounter 2 sand bodies, and the well section also has HD4-52 well 5081.5-5091.17m well section, the well section encounters 4 small layers, wherein 1 small layer encounters 1 sand body, and 2 small layers encounter 2 sand bodies.
Aiming at the geological condition of a research area, determining the type and the grade of an interlayer according to core observation, the distribution position and the production shape of the interlayer, wherein the interlayer of the research area has three types of argillaceous substances, calcareous substances and calcareous argillaceous substances, the lithology of the argillaceous interlayer is gray green argillaceous rock, argillaceous siltstone and argillaceous fine sandstone, and a large amount of argillaceous substances are filled among rock particles; the calcareous interlayer lithology is gray calcareous siltstone and calcareous fine sandstone, the calcareous content is high, and partial hydrochloric acid drops react violently; the lithology of the calcium argillaceous interlayer is argillaceous rock cemented by calcium, sandstone containing silt and fine sandstone, and the arrangement of the rock is compact. The research area separating interlayer comprises three levels of interlayer, sand body and single sand body. The interlayer exists at the interface of the small layer and is approximately parallel in appearance; the interlayer of the sand body in the layer exists at the superposed interface of the single sand body and the single sand body in the small layer, and the shape is inclined; the interlayer in the single sand body is positioned in the single sand body and is inclined in shape.
② the logging response lower limit standards of different types of interlayer are established by a core calibration logging method, the shale interlayer presents the characteristics of high natural gamma value (more than 38API), high density value (2.38-2.58 g/cm3) and relatively high resistivity value (more than 0.6 omega.m) on a conventional logging curve, the calcareous interlayer presents the characteristics of low natural gamma value (less than 38API), high density value (more than 2.48g/cm3) and high resistivity value (more than 0.9 omega.m) on the conventional logging curve, the calcareous interlayer presents the characteristics of low natural gamma value (less than 38API), low density value (less than 2.48g/cm3) and high resistivity value (more than 0.7 omega.m) on the conventional logging curve, the logging response characteristic values of natural gamma GR, density DEN, formation RT and dip angle logging curves at the interlayer top and bottom interfaces, and inflection points of H112, H4-52 interlayer depth types and HD4 interpretive depth interpretation results are shown in an HD-52 table 112-52.
Table 1:
Figure BDA0001245087870000131
③ through continuous adjustment, the final preferred explanation processing parameters can be 0.2m window length, 0.02m step length, 30 degrees exploration angle, processing dip logging curves of HD112 and HD4-52 mesh well sections, counting the azimuth frequency diagram of the interlayer and the main peak direction of the Schmidt vector diagram to determine the production state of the interlayer (Table 1), and FIG. 3 is a processing result diagram of the dip logging curves of HD4-52 single well interlayer, wherein the interlayer with smaller dip angle has higher possibility of continuous distribution, the interlayer with larger dip angle has lower possibility of continuous distribution, but the specific continuity degree needs to be researched in combination with the horizontal well interlayer.
④ fig. 4 is an explanatory result diagram of HD112H well horizontal section interlayer, the mud interlayer reflects the obvious increase of natural gamma curve value and the lower of sound wave time difference curve value on the horizontal well logging curve, the calcium interlayer reflects the obvious decrease of natural gamma curve value and the lower of sound wave time difference curve value, the calcium mud interlayer reflects the natural gamma between the mud interlayer and the calcium interlayer and the higher of sound wave time difference curve value, when the deep and shallow resistivity is both greatly increased and the value is similar to the oil layer, the interlayer develops at the well hole, when the deep and shallow resistivity is increased relative to the oil layer and the deep resistivity is larger than the shallow resistivity, the interlayer develops at the near well hole, when the deep resistivity is increased relative to the oil layer resistivity, but the shallow resistivity is not changed relative to the oil layer, the interlayer develops at the far well hole.
According to the horizontal section logging curve, the natural gamma and resistivity curve of the HD112H well 5220-5260m has obvious abnormal values, which reflect that the curve drills into the breccia section at the upper part of the east river sandstone, the horizontal well from 5260m enters into the integral more homogeneous east river sandstone section, the abnormality on the natural gamma and sonic curve can reflect the type and the position of the interlayer, meanwhile, the relative value of the deep resistivity and the oil layer resistivity is large, the response on the shallow resistivity is not obvious, which shows that the interlayer is positioned at a position far away from the well hole, therefore, a argillaceous interlayer within 0.5m away from the well hole can be picked up at 5260-5385 m, picking up a calcareous interlayer which is about 0.5m away from a borehole at a position of 5275-5310 m, picking up a calcareous interlayer which is 0.5m away from the borehole at 5385-5615 m, connecting with a muddy interlayer of 5260-5385 m, and identifying a argillaceous interlayer which is 0.5m away from the borehole at the position of 5440-5615 m. Due to the detection characteristic of the horizontal well logging curve, only the interlayer within the detection range can be identified, other positions are not necessarily developed, but exceed the detection range and cannot be identified, and therefore the interlayer matching with the pilot hole is needed to be studied.
⑤ fig. 5 is a diagram of matching results of HD112H well pilot hole wells and horizontal well interval interlayers, a horizontal well HD112H drill encounters 1 small layer of 1-3 sand bodies, 2-2 sand bodies and 2-3 sand bodies of 2 small layers, a corresponding pilot hole HD112 projection orientation is 22 °, a formula (1) is used for performing attitude transformation, a first argillaceous interlayer projects to a horizontal section with an attitude of 22 ° ∠.5 ° (inclined to ∠ inclination angle), a second argillaceous interlayer is 22 ° ∠.34 °, a third calcareous interlayer is 22 ° 5630.26 °, inclination angles of the three interlayers are all small, the first argillaceous interlayer is matched with the horizontal well interval interlayer, the first argillaceous interlayer is connected with the argillaceous interlayers of HD112 well sections 5440-5615 m, the second argillaceous sand interval interlayer is connected with the argillaceous interlayer of HD112 well sections 5260-5385 m, the third calcareous interlayer is not connected with the argillaceous interlayer of HD H well sections 5260-5315 m, and no dynamic identification data of the argillaceous interlayer is needed in a dynamic range of HD112, and no identification data of HD 5375 is needed.
⑥ fig. 6 is HD112H water injection curve and HD4-52H production curve, in the research area, the interlayer is continuous, the water injection takes effect when the perforation corresponds, the effective speed of the tracer is faster, at this time, the oil production well presents the characteristics of short anhydrous oil recovery period, fast water content rise after water breakthrough, and daily liquid production rise, when the perforation does not correspond, the water injection does not take effect, at this time, the oil production well presents the characteristics of long anhydrous oil recovery period and gradual liquid production fall.
The HD112H well is transferred to be a water injection well in 2008 in 4 months, a 5339.42-5624.22m well section, namely 2-2 and 2-3 sand bodies in 2 small layers, is jetted, the HD4-52H is jetted, in 2009 in 5 months, a 5195-5498m well section, namely 2-2 and 2-3 sand bodies in 2 small layers, the perforating layer position between the two wells corresponds, the water injection position is lower than the oil production position, and the characteristic of low injection and high production is presented. The tracer is put into the HD112H at 26.11.2010 months, the tracer appears after 276 days in the HD4-52H, the effect taking speed is high and reaches 3.84m/d, the well is put into production, the water content at the initial stage rises quickly, the water content rises from 27.9% to 55% from 5 to 8 months in 2009, and the characteristic of rising is also presented at the initial stage of daily fluid production.
⑦ fig. 7 is HD 112-HD 112H-HD 4-52 well interlayer prediction result diagram, the above production characteristics conform to the 'water injection effect at perforation corresponding time, the effect velocity of the tracer is fast (more than 3.5m/d), the oil production well presents the characteristics of short anhydrous oil extraction period, fast water content rise after water breakthrough and daily production fluid rise' which indicate that the injected water enters along the layer suddenly and the inter-well sand body connectivity is good.
The injection amount of the HD112H well in 2 months of 2009 is 170m3Down to 120m3Later 8 months (18 months in total) in 2010 back to 150m3And then remains stable; the water content of the HD4-52H well correspondingly decreases from 50% to 37% in 11 months in 2009, and then increases to 54% in 4 months (17 months in total) in 2011, and then stabilizes at about 60%, and the time interval between the two is 9 monthsAnd the time is consistent with the time for taking effect of the tracer, so that the oil-water flows smoothly between layers and no channeling occurs, and the production measure change of the water injection well can be reflected on the oil production well, which means that the continuity between layers and the interlayer between sand bodies is good and the interlayer inside a single sand body is discontinuous. The previous prediction result is verified, and the multi-data comprehensive horizontal well interval interlayer prediction method is high in precision and can meet production requirements.
It can be seen that, by the above specific examples, the embodiments of the present invention have the following advantages:
the method for systematically representing the static geological characteristics of the interlayer is creatively provided, the type, the grade, the single well explanation and the attitude pickup of the interlayer are carried out by utilizing rock cores, conventional well logging and dip angle well logging data, and the guiding effect of the static geological characteristics on the interlayer distribution prediction of the horizontal well is effectively improved. Based on rock core data, the types and the grades of the interlayer of the pilot hole well are divided according to the different lithological characters of the interlayers, the different reflection of the dripping hydrochloric acid and the difference of the distribution positions and the thicknesses, so that the defect that the interlayer in the prior art is only of the type but less of the grade and cannot be used for guiding the distribution prediction of the interlayers is overcome; the longitudinal resolution of the dip angle well logging curve can reach 1cm, the type and the position of the guide eye well interlayer are explained by adopting a conventional well logging curve type determination method and a dip angle well logging curve interface determination method, and the defect that the upper interface and the lower interface of the interlayer with small thickness cannot be accurately identified due to the fact that the conventional well logging resolution is too low in the conventional method is overcome; the dip angle well logging curve can be processed to obtain a stratum dip angle vector diagram, the attitude of the interlayer can be well reflected, the inclination and the dip angle of the interlayer of the pilot hole well are picked up by continuously optimizing processing parameters, and the defect that the spatial attitude characterization of the interlayer in the existing method is insufficient is overcome. The method improves the static geological characterization precision of the interlayer and can meet the production requirement.
The horizontal well interlayer identification and leading hole well and horizontal well interlayer matching method is creatively provided, the type and the position of the interlayer can be directly identified by using a horizontal well logging curve, and the defect that horizontal well logging information is not used for horizontal section interlayer identification and horizontal well interlayer distribution prediction is effectively overcome. The horizontal well and the pilot hole have similar logging response characteristics at the position where the interlayer develops, the range and the type of the interlayer of the horizontal well can be explained by utilizing response abnormity of natural gamma and acoustic curves, and the defect that the interlayer explanation of the horizontal well cannot be carried out in the existing method is overcome; the response of the interlayer to the induction logging depth and shallow resistivity at different positions of a horizontal well bore is different, the position of the interlayer relative to the well bore is identified by adopting the relative variation of the interlayer depth and shallow resistivity and the oil layer resistivity, and a method for identifying the space position of the horizontal well interlayer is established; through the attitude conversion of the pilot hole interlayer, based on the attitude extension trend of the pilot hole interlayer and the interpretation achievement of the horizontal well interlayer, the continuity research of the interlayer wells is carried out, and a method for matching the pilot hole with the horizontal well interlayer is provided. The invention provides a method for explaining the interlayer of the horizontal well, which can meet the production requirement.
The dynamic characterization method of the interlayer is creatively provided, the distribution of the interlayer is reversely deduced by using the production data of the oil-water well, and the reliability of the prediction of the interlayer of the horizontal well is improved. By utilizing the corresponding relation of the perforation layer positions and the injection and production effect conditions, the spread of the interval interlayer of the water injection well and the oil production well is researched, and the defect that the injection and production correspondence and effect are only used for evaluating the sand body connectivity but not used for predicting the well interval interlayer distribution in the existing method is overcome; the relation between the water content, the daily fluid production curve change and the interlayer distribution is summarized, the method for establishing the oil-water well production characteristics corresponding to different interlayer distribution modes is established, and the defect that the production curve is difficult to be used for predicting the interlayer distribution of the horizontal well in the conventional method is overcome.
The invention provides a dynamic characterization idea of an interlayer and a multi-data comprehensive horizontal well interlayer prediction method, which can comprehensively apply various data to the maximum extent, reduce the uncertainty in interlayer inter-well prediction, ensure that the prediction result is more reasonable and can be used in production practice.
Corresponding to the method embodiment described in fig. 1 and 2, as shown in fig. 8, an embodiment of the present invention provides a horizontal well interval interbed prediction apparatus, including:
the core analysis unit 301 is configured to obtain core data in a target research area, determine the type, distribution position, and attitude of the pilot hole interlayer according to the core data, and determine the level of the pilot hole interlayer according to the distribution position and attitude of the pilot hole interlayer.
And the pilot hole interlayer type and interface depth determining unit 302 is used for explaining the type of the pilot hole interlayer by using a conventional logging curve according to a preset response lower limit standard strategy of different types of interlayers, and determining the interface depth of the pilot hole interlayer by using an inclination angle logging curve.
And the inclination and inclination angle determining unit 303 is used for continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination and inclination angle of the pilot hole interlayer is determined by determining the interface depth of the pilot hole interlayer by using the inclination angle logging curve to correspond to the production state of the pilot hole interlayer determined according to the core data.
The horizontal well interlayer information determining unit 304 is used for acquiring logging response abnormal data of a horizontal well section, determining the range and the type of the horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data of the horizontal well interlayer, the shallow resistivity data and the relative variation of the resistivity of an oil layer.
And the inclination and inclination angle conversion unit 305 is used for converting the inclination and inclination angle of the pilot hole interlayer to the section where the horizontal well track is located so as to complete the matching of the pilot hole interlayer and the horizontal well interlayer.
And the production characteristic determining unit 306 is used for determining the relationship between the production data of the water injection well and the oil production well and the range and the type of the horizontal well interlayer and the position of the horizontal well interlayer relative to the horizontal well borehole in advance, and determining the production characteristics corresponding to the distribution modes of the different horizontal well interlayers.
And the horizontal well section interlayer distribution prediction unit 307 is used for acquiring production data of the water injection well and the oil production well in the target research area, determining the production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to the matching result of the pilot hole well interlayer and the horizontal well interlayer, and determining the horizontal well section interlayer distribution of the water injection well and the oil production well in the target research area.
Specifically, the types of the guide hole well interlayer in the core analysis unit 301 include argillaceous substances, calcareous substances, and calcareous argillaceous substances; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
In addition, the guidance hole interlayer type and interface depth determination unit 302 is specifically configured to:
according to preset lower limit standard strategies for response of different types of isolation interlayers, the types of the isolation interlayer of the pilot hole well are explained by using the natural gamma logging response characteristic value, the density logging response characteristic value and the formation resistivity logging response characteristic value of the isolation interlayer of the pilot hole well; the preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers.
And determining the interface depth of the pilot hole well interlayer according to the inflection point of the dip angle logging curve at the top and bottom interface of the pilot hole well interlayer.
In addition, the lead hole compartment dip and inclination determination unit 303 is specifically configured to:
and determining a formation dip angle vector diagram according to the dip angle logging curve.
And continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the interface depth of the pilot hole interlayer determined by using the stratigraphic dip angle vector diagram corresponds to the pilot hole interlayer attitude determined according to the rock core data.
And determining an azimuth frequency diagram and a Schmidt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram as the inclination and dip angle of the pilot hole interlayer.
In addition, the inclination and tilt angle conversion unit 305 is specifically configured to:
according toFormula (II):
Figure BDA0001245087870000181
and converting the inclination and the inclination angle of the pilot hole interlayer to the section where the horizontal well track is located to complete the matching of the pilot hole interlayer and the horizontal well interlayer, wherein α is an apparent inclination angle, theta is a true inclination angle, and omega is an included angle between an apparent inclination line and a true inclination line.
In addition, the production characteristic determining unit 306 is specifically configured to:
and according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
The specific implementation manner of the horizontal well interval interlayer prediction device provided by the embodiment of the invention can be seen in the method embodiments corresponding to fig. 1 and fig. 2, and details are not repeated here.
The horizontal well interlayer prediction device provided by the embodiment of the invention comprises the steps of firstly obtaining core data in a target research area, determining the type, the distribution position and the production shape of a pilot hole interlayer according to the core data, and determining the level of the pilot hole interlayer according to the distribution position and the production shape of the pilot hole interlayer; according to a preset lower limit standard strategy of response of different types of interlayer layers, the type of the guide eye well interlayer is explained by using a conventional logging curve, and the interface depth of the guide eye well interlayer is determined by using an inclination logging curve; continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the rock core data; acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the oil layer resistivity; converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer; the method comprises the steps of predetermining the relation between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes; the method comprises the steps of obtaining production data of a water injection well and a production well in a target research area, determining production characteristics of the water injection well and the production well in the target research area, matching the production characteristics of the water injection well and the production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to matching results of a pilot hole interlayer and a horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the production well in the target research area. The method can solve the problem that the method for accurately explaining the interlayer and describing the spatial distribution of the interlayer is lacked at present, can improve the problems of single utilization of data and simplification of a prediction method in the prior art, can overcome the defect of low prediction precision of the distribution of the interlayer of the horizontal well, and improves the prediction precision of the distribution of the interlayer of the horizontal well.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (12)

1. A horizontal well interval interlayer prediction method is characterized by comprising the following steps:
obtaining core data in a target research area, determining the type, the distribution position and the production shape of a pilot hole well interlayer according to the core data, and determining the grade of the pilot hole well interlayer according to the distribution position and the production shape of the pilot hole well interlayer;
according to preset lower limit standard strategies for response of different types of isolation interlayers, the type of the isolation interlayer of the pilot hole well is explained by using a natural gamma logging response characteristic value, a density logging response characteristic value and a formation resistivity logging response characteristic value of the isolation interlayer of the pilot hole well, and the interface depth of the isolation interlayer of the pilot hole well is determined by using an inclination angle logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination angle and the inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the attitude of the pilot hole interlayer determined according to the rock core data;
acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the oil layer resistivity;
converting the inclination and the inclination angle of the pilot hole interlayer to a section where the horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer;
the method comprises the steps of predetermining the relation between the production data of a water injection well and a production well and the range and type of a horizontal well interlayer and the position of the horizontal well interlayer relative to a horizontal well hole, and determining the production characteristics corresponding to different horizontal well interlayer distribution modes;
the method comprises the steps of obtaining production data of a water injection well and a production well in a target research area, determining production characteristics of the water injection well and the production well in the target research area, matching the production characteristics of the water injection well and the production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to matching results of a pilot hole interlayer and a horizontal well interlayer, and determining horizontal well section interlayer distribution of the water injection well and the production well in the target research area.
2. The horizontal well interval interbed prediction method of claim 1, wherein the types of pilot well interval interbed comprise argillaceous, calcareous, and calcic argillaceous; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
3. The method for predicting the interval interlayer of the horizontal well according to claim 2, wherein the method for predicting the interval interlayer of the pilot hole well is characterized in that the type of the interval interlayer of the pilot hole well is explained by using a natural gamma logging response characteristic value, a density logging response characteristic value and a formation resistivity logging response characteristic value of the interval interlayer of the pilot hole well according to a preset response lower limit standard strategy of different types of interval interlayers, and the interface depth of the interval interlayer of the pilot hole well is determined by using an inclination angle logging curve, and comprises the following steps:
determining the interface depth of the pilot hole well interlayer according to the inflection point of the dip angle logging curve at the top-bottom interface of the pilot hole well interlayer; the preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers.
4. The method for predicting the interval interlayer of the horizontal well according to claim 3, wherein the window length, the step length and the exploration angle corresponding to the pilot well are continuously adjusted according to the type and the interface depth of the interval interlayer of the pilot well, so that the inclination angle and the dip angle of the interval interlayer of the pilot well are determined by using the dip angle logging curve, the interface depth of the interval interlayer of the pilot well corresponds to the production state of the interval interlayer of the pilot well determined according to the core data, and the method comprises the following steps:
determining a formation dip angle vector diagram according to the dip angle logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the interface depth of the pilot hole interlayer determined by using the stratigraphic dip angle vector diagram corresponds to the pilot hole interlayer attitude determined according to the core data;
and determining an azimuth frequency diagram and a Schmidt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram as the inclination and dip angle of the pilot hole interlayer.
5. The horizontal well interval interlayer prediction method according to claim 4, wherein the method for converting the inclination and the dip angle of the pilot hole interval interlayer to the section where the horizontal well track is located so as to complete the matching of the pilot hole interval interlayer and the horizontal well interval interlayer comprises the following steps:
according to the formula:
Figure FDA0002232697450000021
and converting the inclination and the inclination angle of the pilot hole interlayer to the section where the horizontal well track is located to complete the matching of the pilot hole interlayer and the horizontal well interlayer, wherein α is an apparent inclination angle, theta is a true inclination angle, and omega is an included angle between an apparent inclination line and a true inclination line.
6. The method for predicting the horizontal well interval interlayer according to claim 5, wherein the relation between the production data of the water injection well and the oil production well and the range and the type of the horizontal well interval interlayer and the position of the horizontal well interval interlayer relative to a horizontal well hole are determined in advance, and the production characteristics corresponding to the distribution modes of the different horizontal well interval interlayers are determined, and the method comprises the following steps:
and according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
7. A horizontal well interval interlayer prediction device is characterized by comprising:
the core analysis unit is used for acquiring core data in a target research area, determining the type, the distribution position and the production state of the pilot hole interlayer according to the core data, and determining the grade of the pilot hole interlayer according to the distribution position and the production state of the pilot hole interlayer;
the pilot hole interlayer type and interface depth determining unit is used for explaining the type of the pilot hole interlayer by utilizing the natural gamma logging response characteristic value, the density logging response characteristic value and the formation resistivity logging response characteristic value of the pilot hole interlayer according to preset different types of interlayer response lower limit standard strategies and determining the interface depth of the pilot hole interlayer by utilizing an inclination angle logging curve;
the inclination and inclination angle determining unit is used for continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole according to the type and the interface depth of the pilot hole interlayer, so that the inclination and inclination angle of the pilot hole interlayer are determined by utilizing the inclination angle logging curve to determine that the interface depth of the pilot hole interlayer corresponds to the production state of the pilot hole interlayer determined according to the rock core data;
the horizontal well interlayer information determining unit is used for acquiring logging response abnormal data of a horizontal well section, determining the range and the type of a horizontal well interlayer according to the logging response abnormal data of the horizontal well section, and determining the position of the horizontal well interlayer relative to a horizontal well hole according to the deep resistivity data and the shallow resistivity data of the horizontal well interlayer and the relative variation of the resistivity of an oil layer;
the inclination and inclination angle conversion unit is used for converting the inclination and inclination angle of the pilot hole interlayer to a section where a horizontal well track is located so as to complete matching of the pilot hole interlayer and the horizontal well interlayer;
the production characteristic determining unit is used for determining the relation between the production data of the water injection well and the oil production well and the range and the type of the horizontal well interlayer and the position of the horizontal well interlayer relative to the horizontal well hole in advance and determining the production characteristics corresponding to the distribution modes of the horizontal well interlayers;
and the horizontal well section interlayer distribution prediction unit is used for acquiring the production data of the water injection well and the oil production well in the target research area, determining the production characteristics of the water injection well and the oil production well in the target research area, matching the production characteristics of the water injection well and the oil production well in the target research area with the production characteristics corresponding to different horizontal well interlayer distribution modes according to the matching result of the pilot hole interlayer and the horizontal well interlayer, and determining the horizontal well section interlayer distribution of the water injection well and the oil production well in the target research area.
8. The horizontal well interbed prediction device of claim 7, wherein the types of leading hole interbed in the core analysis unit comprise argillaceous, calcareous, and calcic argillaceous; the grades of the guide hole interlayer comprise interlayer, in-layer, sand body and single sand body.
9. The horizontal well interval interlayer prediction device according to claim 8, wherein the pilot well interval interlayer type and interface depth determination unit is specifically configured to:
determining the interface depth of the pilot hole well interlayer according to the inflection point of the dip angle logging curve at the top-bottom interface of the pilot hole well interlayer; the preset different types of interlayer response lower limit standard strategies comprise value ranges of natural gamma logging response characteristic values, value ranges of density logging response characteristic values and value ranges of formation resistivity logging response characteristic values corresponding to the different types of interlayers.
10. The horizontal well interval interbed prediction device of claim 9, wherein the leading hole interval interbed dip and inclination determination unit is specifically configured to:
determining a formation dip angle vector diagram according to the dip angle logging curve;
continuously adjusting the window length, the step length and the exploration angle corresponding to the pilot hole, so that the interface depth of the pilot hole interlayer determined by using the stratigraphic dip angle vector diagram corresponds to the pilot hole interlayer attitude determined according to the core data;
and determining an azimuth frequency diagram and a Schmidt vector diagram of the pilot hole interlayer, and counting the main peak direction in the azimuth frequency diagram as the inclination and dip angle of the pilot hole interlayer.
11. The horizontal well interval interbed prediction device of claim 10, wherein the dip and dip conversion unit is specifically configured to:
according to the formula:
Figure FDA0002232697450000041
inclination and inclination of the lead eye well interlayerAnd converting the angle to a section where a horizontal well track is located so as to complete the matching of the pilot hole interlayer and the horizontal well interlayer, wherein α is an apparent dip angle, theta is a true dip angle, and omega is an included angle between an apparent dip line and a true dip line.
12. The horizontal well interval interbed prediction device of claim 11, wherein the production characteristic determination unit is specifically configured to:
and according to the production data of the water injection well and the oil production well, determining the single-well production dynamic response characteristics corresponding to the continuous interlayer spreading, the semi-continuous interlayer spreading and the discontinuous interlayer spreading through perforation correspondence, water injection effect analysis and production curve characteristic research.
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 CN106869915A (en) 2017-06-20
CN106869915B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107589469B (en) * 2017-08-04 2019-07-09 中国石油天然气股份有限公司 Method and device for determining oil-water interface
CN110910267B (en) * 2019-11-18 2023-11-28 中国石油天然气股份有限公司 Method and device for determining sea-phase carbonate rock interlayer
CN113933479B (en) * 2020-07-14 2023-09-26 中国石油天然气股份有限公司 Method and device for distinguishing factors of mudstone interlayer in river-phase stratum

Citations (9)

* 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 stratum in horizontal well/highly deviated well
CN103573250A (en) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 Method for calculating distance from horizontal well hole to upper and lower interfaces of stratum
CN104895554A (en) * 2015-04-30 2015-09-09 中国石油天然气股份有限公司 Depth time conversion method and device for horizontal well logging data
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
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 中国石油天然气股份有限公司 Horizontal well trajectory correction method and device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016053330A1 (en) * 2014-10-01 2016-04-07 Halliburton Energy Services, Inc. Transposition of logs onto horizontal wells

Patent Citations (9)

* 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 stratum in horizontal well/highly deviated well
CN103573250A (en) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 Method for calculating distance from horizontal well hole to 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
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 time conversion method and device for 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 中国石油天然气股份有限公司 Horizontal well trajectory correction method and device

Non-Patent Citations (2)

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

Also Published As

Publication number Publication date
CN106869915A (en) 2017-06-20

Similar Documents

Publication Publication Date Title
CA2895549C (en) Fracturing and reactivated fracture volumes
US11209561B2 (en) Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function
Close et al. Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada
CN106842301B (en) Quantitative identification and prediction method for favorable reservoir of tufaceous sandstone
CA2823710C (en) Methods and systems regarding models of underground formations
Kuanzhi et al. A reserve calculation method for fracture-cavity carbonate reservoirs in Tarim Basin, NW China
CN109870719B (en) Well position layout method, device and system for carbonate rock tight thin reservoir
CN104914465A (en) Volcanic rock crack quantitative prediction method and device
AU2012388240A1 (en) Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models
CN106869915B (en) Horizontal well interval interlayer prediction method and device
Price et al. Influence of stratigraphy on barriers to fracture growth and completion optimization in the Meramec Stack Play, Anadarko Basin, Oklahoma
CN107831540A (en) The direct new method for extracting of reservoir physical parameter
CN109752761B (en) Method and device for evaluating reservoir characteristics of dissolved hydrocarbon reservoir
CN102288996A (en) Method for quantitatively predicting gas containing probability by AFI (prestack inversion)
Trippetta et al. Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy
CN111190228B (en) Submarine fan water channel reservoir configuration identification method
CN106443772B (en) One kind going diapir original stratigraphic thickness restoration methods
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
Jambayev Discrete fracture network modeling for a carbonate reservoir
CN103217718B (en) A kind of method of hidden layer under additional well
Lyu et al. Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China
CN113376692B (en) Method and device for optimizing fracturing modification scheme of tight sandstone gas horizontal well
Yao et al. Using seismic methods to detect connectivity of fracture networks controlled by strike-slip faults in ultra-deep carbonate reservoirs: A case study in northern tarim basin, China
Li et al. Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China
Podnebesnykh et al. New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field

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