WO2017120099A1 - Modélisation permettant de caractériser un réseau de fractures dans des gisements de pétrole homogènes - Google Patents

Modélisation permettant de caractériser un réseau de fractures dans des gisements de pétrole homogènes Download PDF

Info

Publication number
WO2017120099A1
WO2017120099A1 PCT/US2016/069210 US2016069210W WO2017120099A1 WO 2017120099 A1 WO2017120099 A1 WO 2017120099A1 US 2016069210 W US2016069210 W US 2016069210W WO 2017120099 A1 WO2017120099 A1 WO 2017120099A1
Authority
WO
WIPO (PCT)
Prior art keywords
formation
conductivity
model
pressure
fault
Prior art date
Application number
PCT/US2016/069210
Other languages
English (en)
Inventor
Faisal M. ALTHAWAD
Mahmoud Jamiolahmady
Original Assignee
Saudi Arabian Oil Company
Aramco Services Company
Heriot-Watt University
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 Saudi Arabian Oil Company, Aramco Services Company, Heriot-Watt University filed Critical Saudi Arabian Oil Company
Publication of WO2017120099A1 publication Critical patent/WO2017120099A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing 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
    • E21B49/008—Testing 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 by injection test; by analysing pressure variations in an injection or production test, e.g. for estimating the skin factor

Definitions

  • the present invention relates to modeling the structure of subsurface reservoirs, and more particularly to forming models of fracture networks in a reservoir based on pressure transient test measurements obtained from a formation layer of interest in the reservoir.
  • Direct methods use direct measuring tools such as well logging tools. However, the ability of such tools to obtain data measurements as a function of depth into the reservoir from the tools is shallow, typically on the order of a few inches.
  • indirect measurements tools such as pressure gauges are used to record pressure changes due to well rate variations. Indirect measurements involve flowing the well and recording the pressure changes with time. The pressure data obtained are then processed in a number of different ways to describe the reservoir and model the fluid flow processes.
  • Reservoir modeling is, to a great extent, an art and has its benefits and restraints.
  • There are two main methods to model the reservoirs namely; numerical and analytical.
  • Numerical modeling is flexible, however, it can be inaccurate due to instability of computer processing to solve multiple, multi-variable non-linear differential equations expressing the physical relationships of reservoir rock and fluid phenomena and characteristics.
  • numerical models of a reservoir are organized into a large number of individual cells. The number of cells can be from tens to hundreds of millions for typical reservoirs. Instability in the modeling and the gridding effect make numerical modeling unsuitable to address the more complex cases.
  • the present invention provides a new and improved computer implemented method of determining a model of a subsurface earth formation, having a well intersecting a complex geological network of a fault and a fracture in the formation, based on a pressure transient test of the formation.
  • the computer implemented method obtains a test measure of bottom-hole pressure and also obtains a test pressure derivative at sampled instants of measurement during a pressure transient test of the formation.
  • An estimated reference type- curve set of a plurality of type-curves of well pressures in the well and pressure derivatives for selected test values of reservoir capacity, fracture conductivity, fault conductivity value and distance from the tested well and the formation capacity are received.
  • a model well pressure of the formation is determined based on the test measure of well pressure and an estimated value of fracture conductivity, fault conductivity, distance and the formation capacity.
  • a model pressure derivative is then determined based on the test measure of well pressure and the estimated values in the formation.
  • a model type-curve is then formed of the determined model well pressure of the formation and the model pressure derivative. The model type-curve of the determined model well pressure of the formation and the model pressure derivative are then determined with the plurality of type-curves of the estimated reference type-curve set.
  • the estimated value of fracture conductivity and the estimated value of fault conductivity in the formation in addition to its proximity and reservoir capacity of the matched type-curve of the estimated reference type-curve set are stored as models of the fracture conductivity and the fault conductivity of the formation. If not, the estimated values are adjusted, and the steps of determining a model well pressure, determining a model pressure derivative, forming a model type-curve and comparing are repeated based on the adjusted estimated value of fracture conductivity, fault conductivity and its remoteness from the tested well along with the formation capacity.
  • the present invention also provides a new and improved data processing system for determining a model of a subsurface earth formation, having a well intersecting a complex geological network of a fault and a fracture in the formation, based on a pressure transient test of the formation.
  • the data processing system includes a processor, which obtains a test measure of bottom-hole pressure in the well and also a test pressure derivative at sampled instants of measurement during the pressure transient test of the formation.
  • the processor receives an estimated reference type-curve set of a plurality of type-curves of well pressures in the well and pressure derivatives for selected test values of fracture conductivity fault conductivity in the formation.
  • the processor determines a model well pressure of the formation based on the test measure of well pressure and an estimated value of fracture conductivity and an estimated value of fault conductivity and its proximity along with the reservoir capacity in the formation, and also determines a model pressure derivative based on the test measure of well pressure and the estimated value of fracture conductivity and the estimated value of fault conductivity in the formation.
  • the processor then forms a model type-curve of the determined model well pressure of the formation and the model pressure derivative.
  • the processor next compares the model type-curve of the determined model well pressure of the formation and the model pressure derivative with the plurality of type-curves of the estimated reference type-curve set.
  • the processor stores the estimated value of fracture conductivity and the estimated value of fault conductivity in the formation of the matched type-curve of the estimated reference type-curve set as models of the fracture conductivity and the fault conductivity and distance to fault in the formation.
  • the processor adjusts one or all of the estimated value of fracture conductivity, estimated value of fault conductivity and distance to the fault in the formation from the tested well and the quality of formation, and repeats the steps of determining a model well pressure, determining a model pressure derivative, forming a model type-curve and comparing based on the adjusted estimated value of fracture conductivity, fault conductivity and distance of the formation and other values.
  • a memory of the data processing system stores the estimated value of fracture conductivity and the estimated value of fault conductivity in the formation of the matched type-curve of the estimated reference type-curve set as models of the fracture conductivity and the fault conductivity of the formation.
  • the present invention further provides a new and improved data storage device having stored in a non-transitory computer readable medium computer operable instructions for causing a data processing system to determine a model of a subsurface earth formation, having a well intersecting a complex geological network of a fault and a fracture in the formation, based on a pressure transient test of the formation, the instructions stored in the data storage device causing the data processing system to perform a sequence of processing steps.
  • a test measure of bottom-hole pressure is obtained, and a test pressure derivative is also obtained.
  • An estimated reference type-curve set of a plurality of type-curves of well pressures in the well and pressure derivatives for selected test values of fracture conductivity, fault conductivity and distance with reservoir capacity in the formation are received.
  • a model well pressure of the formation is determined based on the test measure of well pressure and an estimated value of fracture conductivity, fault conductivity, distance and reservoir capacity in the formation.
  • a model pressure derivative is then determined based on the test measure of well pressure and the estimated value of fracture conductivity, the estimated value of fault conductivity and distance with reservoir capacity in the formation.
  • a model type- curve is then formed of the determined model well pressure of the formation and the model pressure derivative based on the determined model well pressure of the formation.
  • the model type-curve of the determined model well pressure of the formation and the model pressure derivative are then determined with the plurality of type-curves of the estimated reference type-curve set.
  • the estimated value of fracture conductivity and the estimated value of fault conductivity in the formation of the matched type-curve of the estimated reference type-curve set are stored as models of the fracture conductivity and the fault conductivity of the formation. If not, one or both of the estimated value of fracture conductivity, an estimated value of fault conductivity and distance with reservoir capacity in the formation are adjusted, and the steps of determining a model well pressure, determining a model pressure derivative, forming a model type-curve and comparing are repeated based on the adjusted estimated value of fracture conductivity, fault conductivity and distance with reservoir capacity of the formation.
  • Figure 1 is a schematic view, taken in cross-section, of a producing fractured well in a reservoir formation in the earth having a geologic fault nearby.
  • Figure 2 is a plan view of the flow geometry of the well and reservoir formation of Figure 1.
  • Figure 3 is a functional block diagram of a flow chart of data processing steps for developing type-curves of pressure and pressure derivatives as functions of time for different fracture and fault conductivities according to the present invention.
  • Figure 4 is a functional block diagram of a flow chart of data processing steps for modeling to characterize fractures networks according to the present invention.
  • Figure 5 is a schematic diagram of a data processing system for modeling to characterize fractures networks according to the present invention.
  • Figure 6 is a plot of a type-curve of dimensionless pressure and its log pressure derivative as functions of dimensionless time for a range of different fracture and fault conductivities according to the present invention.
  • Figure 7 is a plot of a type-curve of dimensionless log pressure derivative as functions of dimensionless time for a range of different fracture and fault conductivities and reflecting an early fracture linear flow according to the present invention.
  • Figure 8 is a plot of a type-curve of dimensionless pressure and its log pressure derivative as functions of dimensionless time obtained according to the present invention for a well intersecting a fracture network as compared to a synthetic numerical model of the same well and fracture network.
  • Figure 9 is a plot of a type-curve of dimensionless pressure and its log pressure derivative as functions of dimensionless time obtained according to the present invention for a well intersecting a fracture network as compared to a field example data for such a well.
  • Figures 1 and 2 illustrate a reservoir formation layer R of interest in an example subsurface hydrocarbon producing reservoir with complex flow geometry.
  • a hydrocarbon producing well 10 in a well bore 12 has been drilled into and through the reservoir as well as through subsurface formations above the depth of the reservoir.
  • the producing well is located in a fracture or fracture matrix 14 in the layer R.
  • the fracture 14 is one component of the complex flow geometry. As indicated at 16, the fracture 14 has a fracture width Wf.
  • a fault indicated schematically at 18 is also present in the layer R composed of Region 1, Region 2, and Region 3 nearby the well 10 and the fracture 14.
  • the fault 18 juxtaposes different geometry across a fault between regions identified as Region 2 and Region 3, which are same quality zones in the reservoir.
  • the fault 18 is another component of the complex flow geometry.
  • the fault 18 has a fault width WF.
  • FIG. 2 illustrates schematically the complex flow geometry in the reservoir layer R.
  • the flow domain of the reservoir layer R is composed of three areas: (i) the reservoir layer R with its specific rock and fluid parameters and properties), which defines flow as indicated schematically by arrows 22 into and around the fracture network 14, the fractured well 10 and the nearby fault 18; (ii) the fracture 14 (with its specific fracture properties), which allows fluids to flow as indicated by arrows 24 along and into the fracture 14 and towards the well 10; and (iii) the nearby fault 18 (with its specific properties), which allows fluids to flow as indicated by arrows 26 along, across and towards the fractured well 10.
  • the flow per unit area is defined to be positive from above and negative from below. It should be noted that the present invention provides indications to reservoir engineers and analysts information with regard to the reservoir pressure conditions as functions of both time and space in general, and wellbore pressure conditions with regard to time in particular.
  • the present invention provides a computer implemented methodology of modeling of subsurface reservoirs and formations, and reservoir simulation of such models.
  • the present invention provides improvements to the existing technological processes of characterizing and modeling of subsurface hydrocarbon reservoirs, where complex flow geometry with fractures and faults are present in order to evaluate their development.
  • the present invention is also potentially capable of improving the functioning of computers in performance of reservoir simulation, by reducing the processing time lost due to instability in the simulator processing of the reservoir model.
  • n l, 2, 3 & F
  • FcDf is the dimensionless fracture conductivity described by
  • k r is the matrix dimensionless permeability, and the dimensionless pressure is: k rf . h[p t - p f ]
  • FIG. 3 illustrates a flow chart T setting forth the methodology of the present invention for developing type-curves of pressure and pressure derivatives as functions of time for different fracture and fault conductivities.
  • Figure 4 illustrates a flow chart F setting forth the methodology of the present invention for modeling to characterize fractures networks.
  • the flow chart T ( Figure 3) and the flow chart F ( Figure 4) illustrate the structure of the logic of the present invention as embodied in computer program software.
  • the flow charts illustrate the structures of computer program code elements including logic circuits on an integrated circuit that function according to this invention.
  • the invention is practiced in its essential embodiment by a machine component that renders the program code elements in a form that instructs a digital processing apparatus (that is, a computer) to perform a sequence of data transformation or processing steps corresponding to those shown.
  • the flow chart T of Figure 3 illustrates schematically a preferred sequence of steps of a computer implemented process for developing type-curves of pressure and pressure derivatives as functions of dimensionless time for different fracture and fault conductivities for a subsurface formation or layer of interest.
  • processing according to the present invention begins with data regarding the formation rock, fluid and geometric properties of the layer R and well 10 being stored in memory of a data processing system D ( Figure 5).
  • the data regarding formation rock, fluid and geometric properties are obtained from well logs, core samples, fluid analysis reports, pressure transient tests and other conventional sources of such data.
  • a starting value for each of dimensional fracture conductivity Fcf and dimensional fault conductivity FCF are initialized.
  • model values of the pressure derivative are obtained by the data processing system D based on a specified input value of dimensional fracture conductivity Fcf and dimensional fault conductivity FC F and the pressure transient test data obtained for the region of interest.
  • the derivative is calculated using a computer code that multiplies the dimensionless time (to), in terms of Laplace transform parameter (s), by the change in well pressure with respect to time to produce the well test derivative and plot it in a log-log scale.
  • the model values of the pressure derivative for the specified input value of dimensional fracture conductivity Fcf and dimensional fault conductivity FC F during step 44 are stored in memory of the data processing system D, together with the input value of dimensional fracture conductivity Fcf and dimensional fault conductivity FC F ⁇
  • the values of fracture conductivity and fault conductivity are adjusted as required for a range of postulated values deemed likely to be present based on the pressure transient test data, and additional model values of pressure derivative obtained as described above in step 44 and stored in memory of the data processing system D.
  • a set of model values for type-curves are stored in the data processing system D and are available as indicated at step 50 for presentation as output displays from the data processing system D for analysis by reservoir engineers and analysts.
  • the flow chart F of Figure 4 illustrates schematically a preferred sequence of steps of a computer implemented process for modeling to characterize fractures networks.
  • a time range is selected from the pressure and time data obtained during pressure transient test of layer R.
  • the model and its structure have been described above in terms of Equation 1 in the Laplace domain.
  • model values of present pressure and pressure derivative with time are obtained according to the methodology of flowchart T of Figure 3 and formatted in a form display in log-log plots, and made available for comparison with actual test data and for output display as diagnostic plots by data processing system D ( Figure 5) in such format.
  • the petrophysical and reservoir data of both the well 10, fracture 14 and fault 18 are read in from storage memory for processing in the data processing system D.
  • step 66 actual values for well pressure and pressure derivative are obtained according to actual measured well pressure transient tests data according to Equation (1) above in the data processing system D.
  • step 68 model pressure and derivative plots based on actual pressure transient testing are generated and then made ready to compare with the model pressure and derivative of the data obtained during step 62.
  • the well pressure and pressure derivative values determined during step 68 are also formatted in a form for storage and subsequent display in log-log plots, and are available in that format for output display by data processing system D.
  • step 70 the values for well pressure and pressure derivative determined from actual pressure transient test data results during step 66 are compared with the model values of well pressure and pressure derivative in the log-plot format resulting from step 62. This comparison is done by super-imposing of pressure data from the actual pressure transient test data results on the proposed type curve or reservoir model.
  • Step 72 involves an evaluation of the results of comparison step 72. If the well pressure and pressure derivative values obtained during step 68, which are compared with model values during step 70 indicate that the generated actual values being compared do correspond within a specified acceptable degree to the model data, an acceptable value of well pressure and pressure derivative is indicated.
  • step 74 the fracture and fault parameters indicated by the model are reported as those for the layer being analyzed.
  • step 76 the fracture and fault parameters are displayed as results from the data processing system D.
  • step 72 If the results of step 72 indicate an unacceptable match between pressure or pressure derivative, or both of them, in the measured data and that of the model values being compared, the value of either or both of the dimensional fracture conductivity Fcf and dimensional fault conductivity FCF are adjusted during step 78.
  • the distance to fault dp and matrix permeability k may also be adjusted during step 78.
  • Processing returns to step 66 for processing of the actual well data based on the adjusted values of fracture and/or fault conductivity. Processing continues for further iterations until during step 72 an acceptable agreement is achieved between the measured data and the model data. This indicates, as noted, that the dimensional fracture conductivity and dimensional fault conductivity values of the current iteration are proper indications of the complex flow geometry.
  • the data processing system D includes a computer 100 having a processor 102 and memory 104 coupled to the processor 102 to store operating instructions, control information and database records therein.
  • the data processing system D may be a multicore processor with nodes such as those from Intel Corporation or Advanced Micro Devices (AMD), an HPC Linux cluster computer or a mainframe computer of any conventional type of suitable processing capacity such as those available from International Business Machines (IBM) of Armonk, N.Y. or other source.
  • the data processing system D may also be a computer of any conventional type of suitable processing capacity, such as a personal computer, laptop computer, or any other suitable processing apparatus. It should thus be understood that a number of commercially available data processing systems and types of computers may be used for this purpose.
  • the processor 102 is, however, typically in the form of a personal computer having a user interface 106 and an output display 108 for displaying output data or records of processing of force measurements performed according to the present invention.
  • the output display 108 includes components such as a printer and an output display screen capable of providing printed output information or visible displays in the form of graphs, data sheets, graphical images, data plots and the like as output records or images.
  • the user interface 106 of computer 100 also includes a suitable user input device or input/output control unit 110 to provide a user access to control or access information and database records and operate the computer 100.
  • Data processing system D further includes a database 114 stored in memory, which may be internal memory 114, or an external, networked, or non-networked memory as indicated at 116 in an associated database server 118.
  • the database 114 also contains various data including the time and pressure data obtained during pressure transient testing of the layer under analysis, as well as the rock, fluid and geometric properties of layer R and well 10, and other formation properties, physical constants, parameters, data measurements identified above with respect to Figure 1 and 2 and the Nomenclature table.
  • the data processing system D includes program code 120 stored in a data storage device, such as memory 104 of the computer 100.
  • the program code 120 is in the form of computer operable instructions causing the data processor 102 to perform the methodology of modeling to characterize fractures network in homogeneous petroleum reservoirs as shown in Figures 3 and 4.
  • program code 120 may be in the form of microcode, programs, routines, or symbolic computer operable languages that provide a specific set of ordered operations that control the functioning of the data processing system D and direct its operation.
  • the instructions of program code 120 may be stored in non-transitory memory 104 of the computer 100, or on computer diskette, magnetic tape, conventional hard disk drive, electronic read-only memory, optical storage device, or other appropriate data storage device having a computer usable medium stored thereon.
  • Program code 120 may also be contained on a data storage device such as server 118 as a non-transitory computer readable medium, as shown.
  • the processor 102 of the computer 100 accesses the pressure transient testing data and other input data measurements as described above to perform the logic of the present invention, which may be executed by the processor 102 as a series of computer-executable instructions.
  • the stored computer operable instructions cause the data processor computer 100 to develop type-curves of pressure and pressure derivatives as functions of time for different fracture and fault conductivities according to the methodology of Figure 3 and to develop models to characterize fractures networks according to Figure 4. Results of such processing are then available on output display 108.
  • Figures 6 through 9 are example displays of such results.
  • Figure 6 is a display of a model type-curve set of dimensionless time versus dimensionless pressure and its log-derivative for a selected set of different fracture and fault conductivities which have been determined in a data processing system according to the process of Figure 3 described above.
  • Figure 6 shows determined model pressure at 200 from well data and the determined model pressure derivative is shown at 202. Radial flow is indicated at 204 in Figure 6.
  • Model pressure type-curves for a set of selected values of dimensional fracture conductivity F Cf values of O. lel, 1.5e2 and 1.5e4 are shown at 210, 212, and 214, respectively, in Figure 6.
  • Model pressure derivative type curves for a set of selected values of dimensional fault conductivity F CF values of le8, le9 and lelO are shown at 216, 218, and 220, respectively.
  • the type-curves displayed in Figure 6 shows some distinctive features of flow conditions of subsurface wells where complex fracture networks may be occurring.
  • fractured well pressure behavior at early times such as at 222, where the type-curve indicates a 1 ⁇ 4 slope is of interest.
  • This segment of the type-curves is indicative of a bilinear flow well behavior reflecting two linear flow regimes along and into a fracture such as shown schematically at 14 in Figures 1 and 2.
  • the type-curves of Figure 6 indicate at 224 well behavior of a radial flow demonstrating transient flow in the matrix around a fracture, such as that shown schematically at 22 in Figures 1 and 2.
  • the subsequent increase or up-turn in Figure 6 indicates a bilinear flow regime and also indicating the finite nature of fault 18.
  • the type-curves of Figure 6 further indicate at 204 for subsequent time a radial flow regime of the bounding blocks of Regions 1 and Region 2.
  • Figure 7 is a display of a model type-curve set of dimensionless time versus dimensionless pressure log-derivative for a different selected set of different dimensional fracture conductivities F cf and dimensional fault conductivities F CF which have also been determined in a data processing system according to the process of Figure 3 described above.
  • the pressure derivative type-curves of Figure 7 are for dimensional fracture conductivity F cf of le3, le4, le5, and le6 as indicated and for dimensional fault conductivity FCF of le8, le9, and lelO, as indicated. It is noted that the pressure derivative type-curves of Figure 7 exhibit at 230 a distinctive feature of an early fracture linear flow regime at very early times until such a fracture linear flow ends as shown schematically at 232, with bilinear flow starting at 234 in the manner also indicated in the type-curves of Figure 6.
  • the early fracture linear flow feature shown in the type-curves of Figure 7 reflects the first fluid flow into a well from a fracture alone as shown schematically at 10 and 14, and confirms the stability of the solution even at a very early time.
  • the pressure derivative type-curves of Figure 7 also indicate at 236 beginning of a conductive fault and at 238, a bilinear flow regime of the type described in Figure 6.
  • a synthetic numerically-built model of simulated flow geometry with a well intersecting a fracture network was constructed and the pressure data were generated by a backward modelling of the given well rate, fluid, reservoir, fracture and fault parameters and properties.
  • the pressure data for the simulated flow geometry of the model were then analyzed in a commercial well-test package (i. e. ECRIN of KAPPA Associates).
  • the results obtained for the numerical model are shown at 250 in Figure 8 for pressure and at 252 for pressure derivative. Processing results obtained according to the methods of Figures 3 and 4 indicated a dimensional fracture conductivity Fc f value of le5 and a dimensional fault conductivity FCF value of le8 as indicated in Figure 8.
  • the results are plotted in a log-log plot format as shown at 262 for pressure and 264 for pressure derivative in Figure 8, along with the type-curves as indicated at 250 and 252 from the model data.
  • Initial bi-linear flow is indicated in Figure 8 at 266 in both the model data and the results according to the present invention, followed by beginning of a conductive fault as shown at 268, with a subsequent increase or up-turn shown at 270 indicating bilinear flow, and followed by radial flow as indicated at 272.
  • Heavy lines 280 and 282 have been added in Figure 8 in those portions of the data plots, where substantial conformity exists between the type-curves 250 and 260 and the pressure derivative type-curves 252 and 264 for the fracture conductivity of le5, and for a fault conductivity of le8. As is evident from Figure 8, a good agreement between the curves plotted is noted in those areas.
  • Table 1 contains a comparison of the fracture conductivities and other parameter values both for the synthetic flow geometry network model and according to the present invention. The value of the permeability is determined by dividing the flow capacity (kh) by the layer thickness. The flow capacity is an output of the model selected. Again, good agreement between the values is also indicated.
  • a well model case example form an actual well in a producing field provided a data set for comparison with processing results according to the present invention.
  • the well model corresponds to a vertical well intersecting a fracture network in a tight homogenous reservoir.
  • the case example was to evaluate results according to the present invention in comparison with an existing example of data from an actual well.
  • Pressure transient testing of the actual well has determined the well to exhibit a flow which is dominated by a fracture bi-linear flow regime for both pressure type-curve at 298 and pressure derivative type-curve at 256 in Figure 9, followed by beginning of a conductive fault in pressure derivative type- curve with a negative half-slope as indicated at 290, followed by an increase at 292 indicating a bi-linear flow regime.
  • the present invention provides a new methodology where pressure transient data is processed so that a complex flowing geometry with flow from fractures and faults is rigorously described based on values of fracture conductivities and fault conductivities which are determined.
  • the present invention provides models of the complex flow geometry which conforms to both numerical models and actual field data.
  • the present invention provides reliable reservoir models based on the pressure transient testing of a reservoir.
  • Type-curves such as those shown in Figure 6 indicate how a complex network of fractures appears when recording pressure transient data from oil and gas wells in a reservoir. Accordingly, data that show similar behavior can be matched using such type-curves. Once the match is obtained, the values of the fracture nature and dimensions can be determined with a high degree of accuracy.
  • the present invention thus provides accurate semi-analytical solutions for a well intersecting fractures network in homogenous reservoir(s). This is of considerable value in view of increasing activities in production from naturally faulted geological settings and unconventional reservoirs.
  • the developed present invention offers more flexible schemes to easily carry out modelling with increasing certainty and larger positive impact on the management decisions of such reservoirs.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Mining & Mineral Resources (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)

Abstract

L'invention concerne la création de modèles de systèmes de gisements complexes comprenant des réseaux de fractures et des défauts, à partir de données d'essai transitoire de pression provenant d'un puits situé dans une région d'intérêt du gisement. Une méthodologie de solution analytique est prévue, permettant d'interpréter la signature des données de test de puits à partir des données d'essai transitoire de pression sur la base de données de modèle regroupant la géométrie de flux simulée et le comportement des données de pression.
PCT/US2016/069210 2016-01-04 2016-12-29 Modélisation permettant de caractériser un réseau de fractures dans des gisements de pétrole homogènes WO2017120099A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US14/987,120 US20170191348A1 (en) 2016-01-04 2016-01-04 Modeling to characterize fractures network in homogeneous petroleum reservoirs
US14/987,120 2016-01-04

Publications (1)

Publication Number Publication Date
WO2017120099A1 true WO2017120099A1 (fr) 2017-07-13

Family

ID=57915077

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2016/069210 WO2017120099A1 (fr) 2016-01-04 2016-12-29 Modélisation permettant de caractériser un réseau de fractures dans des gisements de pétrole homogènes

Country Status (2)

Country Link
US (1) US20170191348A1 (fr)
WO (1) WO2017120099A1 (fr)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2018410470B2 (en) * 2018-02-20 2022-03-10 Chevron U.S.A. Inc. Systems and methods for generating permeability scaling functions to estimate permeability
US11035210B2 (en) * 2018-10-22 2021-06-15 Halliburton Energy Services, Inc. Optimized foam application for hydrocarbon well stimulation
US11401803B2 (en) 2019-03-15 2022-08-02 Saudi Arabian Oil Company Determining fracture surface area in a well
CN111950112B (zh) * 2019-05-16 2024-06-11 中国石油化工股份有限公司 一种适用于底部封闭的碳酸盐岩储层动态分析方法
CN114086934B (zh) * 2020-08-03 2024-03-19 中国石油化工股份有限公司 一种试井解释方法、装置、存储介质以及计算机设备
CN112377186B (zh) * 2020-11-30 2024-01-30 中国石油天然气股份有限公司 哑铃状缝洞型碳酸盐岩储层产量分析模型构建方法及装置
CN112682033B (zh) * 2020-11-30 2023-08-22 中国石油天然气股份有限公司 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
US11933165B2 (en) 2021-03-15 2024-03-19 Saudi Arabian Oil Company Hydraulic fracture conductivity modeling

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015191927A1 (fr) * 2014-06-11 2015-12-17 Advantek International Corporation Quantification du volume d'un réservoir et de la limite de pression d'une pompe

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015191927A1 (fr) * 2014-06-11 2015-12-17 Advantek International Corporation Quantification du volume d'un réservoir et de la limite de pression d'une pompe

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ANONYMOUS: "Fluid flow in hydraulically fractured wells", 4 June 2015 (2015-06-04), XP055357892, Retrieved from the Internet <URL:http://petrowiki.org/Fluid_flow_in_hydraulically_fractured_wells> [retrieved on 20170323] *
FIKRI KUCHUK ET AL: "Pressure-Transient Tests and Flow Regimes in Fractured Reservoirs", SPE RESERVOIR EVALUATION AND ENGINEERING, vol. 18, no. 02, 1 May 2015 (2015-05-01), US, pages 187 - 204, XP055357791, ISSN: 1094-6470, DOI: 10.2118/166296-PA *
FIKRI KUCHUK ET AL: "Transient Pressure Test Interpretation for Continuously and Discretely Fractured Reservoirs", 2012 SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION, 8 October 2012 (2012-10-08), San Antonio, Texas, USA, pages 1 - 35, XP055357813 *
M ZEYBEK ET AL: "Fault and Fracture Characterization Using 3D Interval Pressure Transient Tests", 10TH ABU DHABI INTERNATIONAL PETROLEUM EXHIBITION AND CONFERENCE, 13 October 2002 (2002-10-13), Abu Dhabi, pages 13 - 16, XP055357796 *

Also Published As

Publication number Publication date
US20170191348A1 (en) 2017-07-06

Similar Documents

Publication Publication Date Title
WO2017120099A1 (fr) Modélisation permettant de caractériser un réseau de fractures dans des gisements de pétrole homogènes
US10557333B2 (en) Estimating measures of formation flow capacity and phase mobility from pressure transient data under segregated oil and water flow conditions
US10392922B2 (en) Measuring inter-reservoir cross flow rate between adjacent reservoir layers from transient pressure tests
US9896930B2 (en) Three-dimensional reservoir pressure determination using real time pressure data from downhole gauges
Mirzaei-Paiaman et al. A new methodology for grouping and averaging capillary pressure curves for reservoir models
EP3248031B1 (fr) Mesure de débit d&#39;écoulement croisé à travers des fuites non intentionnelles dans des gaines en ciment d&#39;isolation zonale dans des puits décalés
US10119396B2 (en) Measuring behind casing hydraulic conductivity between reservoir layers
EP3047096A1 (fr) Identification d&#39;une structure de profondeur de formation géologique au moyen de données de diagraphie de puits
US9892366B2 (en) Facies definition using unsupervised classification procedures
EP4042211B1 (fr) Modélisation de la perméabilité de réservoir par estimation de la distribution et des propriétés de fracture naturelle
Morton et al. Global sensitivity analysis for natural fracture geological modeling parameters from pressure transient tests
US11487032B2 (en) Characterizing low-permeability reservoirs by using numerical models of short-time well test data
US20220027616A1 (en) Systems and methods for reservoir history matching quality assessment and visualization
EP3714132B1 (fr) Mesure de demi-longueur de fracture efficace et quantification de distribution de flux dans des fractures et autour de ces dernières dans des réservoirs de pétrole
US11493654B2 (en) Construction of a high-resolution advanced 3D transient model with multiple wells by integrating pressure transient data into static geological model
Li et al. Petrophysical Rock Typing Based on Pore Geometry Improves Permeability and Bound Fluid Volume Estimation in Heterogeneous Sandstone Formations
US20180179891A1 (en) Measuring effective fracture half-length and quantifying flux distribution in and around fractures in petroleum reservoirs
Mabrouk Quantifying Formation Damage Due to Drilling Through Constructing Electro-Facies Model from WFT Data and Suite of Well Logs
Ateeq et al. An integrated study to characterize and model natural fracture networks of gas condensate carbonate reservoirs, onshore Abu Dhabi, UAE
Anton et al. Integrated Approach to Investigating the Effects of Structural Fault Uncertainty on Well Test Design
Morton et al. Integrated Interpretation for Pressure Transient Tests in Discretely Fractured Reservoirs (SPE 154531)
Zaafran Importance of Checking Quality Control of Data in Describing Oil and Gas Reservoirs and Its Effect in Calculating Important Reservoir Parameters and Thus, Oil and Gas Reserves

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16831805

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 11.09.2018)

122 Ep: pct application non-entry in european phase

Ref document number: 16831805

Country of ref document: EP

Kind code of ref document: A1