Patent Grant
12360278
The present disclosure relates to the field of well logging to determine characteristics of formations surrounding an earth borehole and, more particularly, to methods for determining earth formation characteristics using dielectric permittivity and conductivity measurements.
The determination of formation water saturation (Sw) is a challenging petrophysical calculation and it is used to quantify information useful for reservoir evaluation and analysis, such as formation hydrocarbon saturation (1−Sw) and then volume of original oil in place (i.e., stock-tank original oil in place) and/or volume of original gas in place (i.e., stock-tank original gas in place).
Resistivity logs have traditionally been used to estimate formation water saturation for reservoir evaluation and analysis. There are also applications of using high frequency permittivity tools for estimation of formation water saturation. Among them the following patents are relevant: U.S. Pat. No. 4,786,873, to Sherman entitled “Method for Evaluating Water Saturation in Petroleum Reservoirs from Dielectric Permittivity” and U.S. Pat. No. 7,363,160 to Seleznev at al. entitled “Technique for Determining Properties of Earth Formations Using Dielectric Permittivity Measurements.”
In these applications, formation water saturation is based on Archie's model or equation as follows:
where Sw is formation water saturation, Rw is formation water resistivity which is a function of water salinity and measurement temperature, RT is true formation resistivity, ϕ is formation porosity, and m and n are formation rock electrical properties usually determined in laboratory core analysis.
The underlying assumption of using Archie's model of Eqn. (1) is that the formation water has lower resistivity than the hydrocarbons (oil and/or gas) of the formation. However, in a fresh-water environment, the formation water can have a higher resistivity, closer to that of hydrocarbons (oil and/or gas) of the formation. In this case, the Archie's equation is not capable to differentiate water from hydrocarbon, thus estimating water saturation using Archie's equation is not possible in such environments. In addition, since formation water salinity or resistivity is a direct input to the Archie's equation, formation water saturation calculated from Archie's equation can carry very large uncertainties in unknown mixed formation water salinity environments.
This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
Method and systems are provided that characterize a geological formation, which involve collecting or obtaining electrical signal data that results from interaction of electrical signals with the formation. The electrical signal data is processed to determine data characterizing relative permittivity of the formation. The data characterizing relative permittivity of the formation is used as an input to at least one model to determine data characterizing water saturation of the formation.
In embodiments, the at least one model is based on a first equation and second equation, wherein the first equation relates a parameter representing relative permittivity of the formation to a parameter representing water volume of the formation, and wherein the second equation relates the parameter representing water volume of the formation and formation porosity to a parameter representing permittivity-based water saturation of the formation.
In one embodiment, the first equation is of the form
VW=Aεr B, where εr is the parameter representing relative permittivity of the formation, Vw is the parameter representing water volume of the formation, and A and B are calibration constants.
In another embodiment, the first equation is of the form Vw2+aVw=bεr+c, where εr is the parameter representing relative permittivity of the formation, Vw is the parameter representing water volume of the formation, and a, b, and c are calibration constants.
In one embodiment, the second equation is of the form where
SW (εr) is the parameter representing permittivity-based water saturation of the formation, Vw is the parameter representing water volume of the formation, and ϕ is formation porosity.
In embodiments, the parameter representing permittivity-based water saturation of the formation can be used to generate the data characterizing water saturation of the formation. For example, in a freshwater environment, the parameter representing permittivity-based water saturation of the formation can be equated to the data characterizing water saturation of the formation or given weight in generating the data characterizing water saturation of the formation.
In embodiments, the at least one model can further include at least one additional equation that relates formation resistivity or conductivity to a parameter representing resistivity-based water saturation of the formation resistivity.
In embodiments, the parameter representing permittivity-based water saturation of the formation and the parameter representing resistivity-based water saturation of the formation can be combined or reconciled to generate the data characterizing water saturation of the formation. For example, in a high salinity environment, the data characterizing water saturation of the formation can be equated to the parameter representing resistivity-based water saturation of the formation, or more weight can be given to the parameter representing resistivity-based water saturation of the formation relative to the parameter representing permittivity-based water saturation of the formation in generating the data characterizing water saturation of the formation. In another example, in a freshwater (low salinity) or unknown or mixed salinity environment, the data characterizing water saturation of the formation can be equated to the parameter representing permittivity-based water saturation of the formation, or more weight is given to the parameter representing permittivity-based water saturation of the formation relative to the parameter representing resistance-based water saturation of the formation in generating the data characterizing water saturation of the formation.
In one embodiment, the at least one additional equation is of the form
where Sw is resistivity-based water saturation of the formation, Rw is formation water resistivity which is a function of water salinity and measurement temperature, RT is true formation resistivity, ϕ is formation porosity, and m and n are formation rock electrical properties.
In embodiments, the processing of the electrical signal data to determine data characterizing relative permittivity of the formation can involve an inversion process, such as a pixel-based inversion process.
In embodiments, the electrical signal data can include in-phase signal component data and out-of-phase signal component data acquired by an induction logging tool. For example, the in-phase signal component data and out-of-phase signal component data are acquired at a lower frequency within the range of a few hundred KHz to a few MHz.
In other embodiments, the electrical signal data can include attenuation signal component data and phase-shift signal component data acquired by a propagation logging tool at frequency within the range of a few hundred KHz to a few MHz.
In embodiments, the methods and systems further involve the operation of a downhole logging tool to acquire the electrical signal data of a) at a depth of investigation of at least 2 ft to 10 ft into the formation from the borehole wall.
In embodiments, the operations of the methods and systems can be performed by a processor.
In embodiments, the data characterizing water saturation of the formation is used for reservoir saturation monitoring (such as identifying by-passed hydrocarbon, monitoring water-hydrocarbon contact movement, and evaluating sweep efficiency by invading water through such as waterflooding in fresh water or mixed salinity environments).
In embodiments, the operations of the methods and systems can be repeated for electrical signal data measured at different depths and different azimuth directions in a wellbore in order to investigate different parts of the formation that is traversed by the wellbore.
In another aspect, methods and systems are provided for drilling a wellbore that traverses a formation using a drill bit, which involve collecting or obtaining electrical signal data that results from interaction of electrical signals with the formation. The electrical signal data is processed to determine data characterizing relative permittivity of the formation. The data characterizing relative permittivity of the formation is used as an input to at least one model to determine data characterizing water saturation of the formation. These operations are repeated for electrical signal data measured at different depths and different azimuth directions in the wellbore in order to investigate different parts of the formation that is traversed by the wellbore. The drill bit is steered to adjust and control trajectory of the drilled wellbore based on the data characterizing water saturation of the different parts of the formation.
The subject disclosure is further described in the detailed description which follows, in reference to the noted plurality of drawings by way of non-limiting examples of the subject disclosure, in which like reference numerals represent similar parts throughout the several views of the drawings, and wherein:
The particulars shown herein are by way of example and for purposes of illustrative discussion of the embodiments of the subject disclosure only and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the subject disclosure. In this regard, no attempt is made to show structural details in more detail than is necessary for the fundamental understanding of the subject disclosure, the description taken with the drawings making apparent to those skilled in the art how the several forms of the subject disclosure may be embodied in practice. Furthermore, like reference numbers and designations in the various drawings indicate like elements.
Dielectric permittivity can be estimated from electromagnetic measurements, such as wireline induction logs and/or propagation logging while drilling tools. Water has a high dielectric permittivity in comparison with other common materials commonly encountered in subsurface, because it can store charges in the dipole-dipole associated with the O—H covalent bond. The relative dielectric permittivity (εr) of the different materials of an earth formation vary widely (for example, 2.2 for oil, 7.5 for limestone, and 80 for water). This strong contrast between water and other materials provides an opportunity for using dielectric permittivity to estimate formation water saturation.
In a subsurface reservoir environment, dielectric permittivity is mainly caused by interfacial polarization (IP). This effect is widely observed from electromagnetic (EM) surveys conducted on porous media. The commonly measured downhole EM measurements have components that are also impacted by IP effects; thus, they can be used to estimate dielectric permittivity. More specifically, in a reservoir formation containing porous media (reservoir rock) and formation fluids (water, oil, and/or gas), application of an external electric field results in electric conduction (migration of charge carries) and interfacial polarization (IP). The conduction current and IP effects can be characterized by formation conductivity and dielectric permittivity, which are best described by a complex conductivity:
σ=σR+i2πfεrε0, Eqn. (2)
where σR is effective formation conductivity, εr relative dielectric permittivity, ε0 the permittivity of vacuum (8.854×10−12) and f measurement frequency.
A pixel inversion-based processing method can be used to simultaneous determine the conductivity σR (or resistivity) and relative dielectric permittivity εr of the formation, by measured R- and X- components of resistivity logs obtained from induction-type logging tools as well as phase shift and attenuation data from Logging While Drilling (LWD) tools. An example of pixel inversion-based processing method is described in detail in U.S. Pat. No. 10,036,826, commonly assigned to assignee of the subject disclosure. In this example pixel inversion-based processing method, there is no assumption as to a maximal number of layers (an assumption usually present in model-based inversion approaches). Furthermore, instead of inverting for distance to boundaries layer thicknesses, bed resistivities, and dip, the example pixel inversion-based processing method inverts for distribution of resistivities using predefined thin layers referred to as “pixels.” The boundary positions are inferred from changes in resistivities.
In other embodiments, model-based inversion processing methodology can be used. The inversion process can employ a model based on Maxwell's equations to solve for the data characterizing resistivity (or conductivity) and complex permittivity (or dielectric constant) of the formation. The stopping criteria for the inversion process can also be based on cost function. The cost function may include one or more parameters related to physical properties of the formation, such as resistivity (or conductivity), permittivity (or dielectric constant), and position within the formation. For example, the cost function may include a misfit term, an entropy term, and a smoothing term. In some embodiments, one or more of the terms (e.g., the misfit term, the entropy term, and the smoothing term) may be parameterized based on the resistivity (or conductivity) and/or permittivity (or dielectric constant) as defined by the model.
Through model simulations, it is found that the dielectric permittivity has strong dependence on porosity and water saturation and little dependence on water salinity. Furthermore, the dielectric permittivity has almost a linear relationship with the volume of water, i.e., the product of porosity and saturation, which makes it possible to estimate water saturation using dielectric permittivity.
Vw=Swϕ. Eqn. (3)
Specifically,
Vw=Aεr+B, Eqn. (4a)
where A and B are calibration constants, which can be determined from linear regression analysis on laboratory measurements of core samples or other suitable measurements.
Specifically,
Vw2+aVw=bεr+C. Eqn. (4b)
where a, b, and c are calibration constants, which can be determined from non-linear regression analysis on laboratory measurements of core samples or other suitable measurements.
Specifically,
Equations 4a and 4b allow the determination of Vw from the relative dielectric permittivity εr. With independently determined formation porosity (ϕ), a permittivity-based water saturation, Sw(εr), can then be calculated based on the definition of Vw (Eq. 3) as follows:
Equations 5 and 4a and/or 4b provide a model for permittivity-based formation water saturation Sw(εr) as a function of relative dielectric permittivity εr which is referred to as a “Deep Reading Relative Permittivity model” or “DRRP model” herein.
The workflow begins in block 501 which determines or quantifies an initial estimate of formation water salinity, which can be obtained from field testing or lab measurements of formation water samples.
In block 503, formation porosity can be determined independently from the estimate of formation water salinity of block 501. For example, formation porosity can be determined from porosity logging operations by acoustic, nuclear (e.g., combination of density and neutron), NMR or other suitable porosity logging operations. The formation porosity of block 503 can be used as input to the Archie-type model (block 507) and the DRRP model (block 509).
In block 505, rock properties or other calibration constants that are used as inputs to the Archie-type model (block 507) and the DRRP model (block 509) are obtained, for example from laboratory core analysis or other suitable methods.
In blocks 511A, 511B and 511C, deep reading electromagnetic log measurements are obtained and processed to extract formation conductivity (and/or resistivity) and formation relative permittivity. The electromagnetic log measurements can be obtained by operation of a low frequency (few tens KHz) induction logging tool (wireline or LWD) or other suitable electromagnetic logging tool (wireline or LWD) with a deep depth of investigation (>2 ft). The formation contains porous media (reservoir rock). Formation fluids, such as connate water, oil, and/or hydrocarbon gas, can be contained within the porous media (reservoir rock). The processing of blocks 511B and 511C can involve pixel inversion-based processing to simultaneous determine the conductivity σR (or resistivity) and relative permittivity εr of the formation based on the measured R- and X- components of resistivity logs obtained from induction-type logging tools (or based on measured phase shift and attenuation data from Logging While Drilling (LWD) tools) as described above. In other embodiments, model-based inversion processing can also possibly be used.
In block 513, the Archie-type model (block 507) is used to determine a resistivity-based formation water saturation, Sw(RT). For example, the Archie-type model of block 507 can be a formation resistivity-based Archie model, such as the model of Eqn. (1) described above where the formation water saturation Sw is determined as function of true formation resistivity RT as determined in block 511B.
In block 515, the DRRP model (block 509) is used to determine a permittivity-based formation water saturation, Sw(εr). For example, the DRRP model of block 507 can be a permittivity-based model, such as the model of Eqns. (5) and 4(a) and/or 4(b) described above where the formation water saturation Sw is determined as function of relative dielectric permittivity εr as determined in block 511C.
In block 517, the resistivity-based formation water saturation Sw(RT) of block 513 and the permittivity-based formation water saturation Sw(εr) of block 515 are integrated and/or reconciled with one another to determine a final formation water saturation Sw based thereon. As part of the reconciliation to determine the final formation water saturation Sw, in a high salinity environment, the final formation water saturation Sw can be equated to Sw(RT), or more weight may be given to Sw(RT) relative to Sw(εr). On the other hand, in a freshwater (low salinity) environment or in an unknown or mixed salinity environment, the final formation water saturation Sw can be equated to Sw(εr) or more weight may be given to Sw(εr) relative to Sw(RT). The weighting factor is a function of water salinity.
The final formation water saturation Sw as determined in block 517 can be used to determine useful formation properties, such as hydrocarbon saturation, original oil in place (i.e., stock-tank original oil in place) and/or volume of original gas in place formation. In another example shown in Eqn. (6) below, by substituting Sw(RT) with Sw(εr) in the Archie-type model (e.g., Eqn. (1)), formation water salinity at a given temperature may also be determined as of block 517. This calculation is useful in unknown mixed salinity environments and estimated qualitatively in freshwater environments.
RW=Swn(εr)ϕmRT. Eqn. (6)
In other embodiments, the final formation water saturation Sw as determined in block 517 can be used for formation evaluation and saturation monitoring, such as identifying by-passed hydrocarbon, monitoring water movement, evaluating water sweep efficiency, identifying an invasion profile in fresh water or mixed salinity environments, identifying fluid and core sampling depths in fresh water or mixed salinity environments, or identifying completion perforation locations in fresh water or mixed salinity environments. For example, completion perforation locations that would interface to parts of the formation that are invaded by relatively high amount of water (as determined from the final formation water saturation Sw) can be avoided, while completion perforation locations that interface to parts of the formation that include a relatively high amount of hydrocarbons (e.g., oil or gas) can be used for completion.
In still other embodiments, the final formation water saturation Sw as determined in block 517 can be used to evaluate formations through horizontal and multilateral wells where effects such as polarization horns, formation anisotropy, etc. on measured formation resistivity are significant.
It will be appreciated that the workflow as described herein can employ a downhole logging tool (such as a wireline or logging-while-drilling tool) that is conveyable in a wellbore that traverses the reservoir rock under investigation, such as the logging tool described below with respect to
It should be appreciated that in an alternative embodiment such processing circuitry is capable of being located downhole in or near the logging tool TL itself. Such processing circuitry being capable of handling all the processing functionality pertaining to the various measurements and models described herein.
Moreover, while
In other embodiments, the induction logging tool can be an array induction logging tool, such as the AIT or Thrubit Induction tools of Schlumberger Technology Corporation. In still other embodiments, the workflow can be based on attenuation and phase shift measurements acquired by a propagation-type downhole tool as a substitute for the R and X signals acquired by the induction logging tool.
Note that the downhole measurement of the workflow can be performed at different depths and different azimuth directions in the wellbore in order to investigate different parts of the geological formation that is traversed by the wellbore, and the processing of the workflow can be repeated for the corresponding downhole measurements at the respective depths and azimuth directions in order to accurately measure water saturation or other formation parameters at the measured depths and azimuthal directions. Such analysis is particularly useful for formations in freshwater and/or unknown mixed salinity environments.
Device 2500 is one example of a computing device or programmable device and is not intended to suggest any limitation as to scope of use or functionality of device 2500 and/or its possible architectures. For example, device 2500 can comprise one or more computing devices, programmable logic controllers (PLCs), etc.
Further, device 2500 should not be interpreted as having any dependency relating to one or a combination of components illustrated in device 2500. For example, device 2500 may include one or more of computers, such as a laptop computer, a desktop computer, a mainframe computer, etc., or any combination or accumulation thereof.
Device 2500 can also include a bus 2508 configured to allow various components and devices, such as processors 2502, memory 2504, and local data storage 2510, among other components, to communicate with each other.
Bus 2508 can include one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. Bus 2508 can also include wired and/or wireless buses.
Local data storage 2510 can include fixed media (e.g., RAM, ROM, a fixed hard drive, etc.) as well as removable media (e.g., a flash memory drive, a removable hard drive, optical disks, magnetic disks, and so forth).
One or more input/output (I/O) device(s) 2512 may also communicate via a user interface (UI) controller 2514, which may connect with I/O device(s) 2512 either directly or through bus 2508.
In one possible implementation, a network interface 2516 may communicate outside of device 2500 via a connected network.
A media drive/interface 2518 can accept removable tangible media 2520, such as flash drives, optical disks, removable hard drives, software products, etc. In one possible implementation, logic, computing instructions, and/or software programs comprising elements of module 2506 may reside on removable media 2520 readable by media drive/interface 2518. Various processes of the present disclosure or parts thereof can be implemented by instructions and/or software programs that are elements of module 2506. Such instructions and/or software programs may reside on removable media 2520 readable by media drive/interface 2518 as is well known in the computing arts.
In one possible embodiment, input/output device(s) 2512 can allow a user (such as a human annotator) to enter commands and information to device 2500, and also allow information to be presented to the user and/or other components or devices. Examples of input device(s) 2512 include, for example, sensors, a keyboard, a cursor control device (e.g., a mouse), a microphone, a scanner, and any other input devices known in the art. Examples of output devices include a display device (e.g., a monitor or projector), speakers, a printer, a network card, and so on.
Various processes or parts of the workflow of the present disclosure may be described herein in the general context of software or program modules, or the techniques and modules may be implemented in pure computing hardware. Software generally includes routines, programs, objects, components, data structures, and so forth that perform particular tasks or implement particular abstract data types. An implementation of these modules and techniques may be stored on or transmitted across some form of tangible computer-readable media. Computer-readable media can be any available data storage medium or media that is tangible and can be accessed by a computing device. Computer readable media may thus comprise computer storage media. “Computer storage media” designates tangible media, and includes volatile and non-volatile, removable and non-removable tangible media implemented for storage of information such as computer readable instructions, data structures, program modules, or other data. Computer storage media include, but are not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other tangible medium which can be used to store the desired information, and which can be accessed by a computer.
In embodiments, any one or any portion or all of the steps or operations of the workflow as described above can be performed by a processor. The term “processor” should not be construed to limit the embodiments disclosed herein to any particular device type or system. The processor may include a computer system. The computer system may also include a computer processor (e.g., a microprocessor, microcontroller, digital signal processor, or general purpose computer) for executing any of the methods and processes described above.
The computer system may further include a memory such as a semiconductor memory device (e.g., a RAM, ROM, PROM, EEPROM, or Flash-Programmable RAM), a magnetic memory device (e.g., a diskette or fixed disk), an optical memory device (e.g., a CD-ROM), a PC card (e.g., PCMCIA card), or other memory device.
Some of the methods and processes described above, can be implemented as computer program logic for use with the computer processor. The computer program logic may be embodied in various forms, including a source code form or a computer executable form. Source code may include a series of computer program instructions in a variety of programming languages (e.g., an object code, an assembly language, or a high-level language such as C, C++, or JAVA). Such computer instructions can be stored in a non-transitory computer readable medium (e.g., memory) and executed by the computer processor. The computer instructions may be distributed in any form as a removable storage medium with accompanying printed or electronic documentation (e.g., shrink wrapped software), preloaded with a computer system (e.g., on system ROM or fixed disk), or distributed from a server or electronic bulletin board over a communication system (e.g., the Internet or World Wide Web).
Alternatively or additionally, the processor may include discrete electronic components coupled to a printed circuit board, integrated circuitry (e.g., Application Specific Integrated Circuits (ASIC)), and/or programmable logic devices (e.g., a Field Programmable Gate Arrays (FPGA)). Any of the methods and processes described above can be implemented using such logic devices.
While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein. Accordingly, the scope of the invention should be limited only by the attached claims. Moreover, embodiments described herein may be practiced in the absence of any element that is not specifically disclosed herein.
In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together, whereas a screw employs a helical surface, in the environment of fastening wooden parts, a nail and a screw may be equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. § 112, paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words ‘means for’ together with an associated function.
There have been described and illustrated herein one or more embodiments of methods and systems that extract deep formation conductivity (or resistivity) and formation relative permittivity from electromagnetic log measurements and process the deep formation conductivity (or resistivity) and formation relative permittivity to determine water saturation of formation rock and possibly other useful formation parameters and fluids parameters, such as reservoir wettability. It will therefore be appreciated by those skilled in the art that yet other modifications could be made to the provided invention without deviating from its spirit and scope as claimed.
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/US2020/020339 | 2/28/2020 | WO |
| Publishing Document | Publishing Date | Country | Kind |
|---|---|---|---|
| WO2021/173150 | 9/2/2021 | WO | A |
| Number | Name | Date | Kind |
|---|---|---|---|
| 4786873 | Sherman | Nov 1988 | A |
| 7363160 | Seleznev et al. | Apr 2008 | B2 |
| 10036826 | Thiel et al. | Jul 2018 | B2 |
| 11119239 | Kouchmeshky | Sep 2021 | B2 |
| 20100078165 | Seleznev | Apr 2010 | A1 |
| 20140320126 | Heaton | Oct 2014 | A1 |
| 20160245952 | Dupuis | Aug 2016 | A1 |
| 20170276832 | Kouchmeshky | Sep 2017 | A1 |
| 20180164466 | Zhang | Jun 2018 | A1 |
| 20190204466 | Zhang | Jul 2019 | A1 |
| 20190265382 | Itskovich et al. | Aug 2019 | A1 |
| Number | Date | Country |
|---|---|---|
| 2018183223 | Oct 2018 | WO |
| Entry |
|---|
| Substantive Exam issued in Saudi Arabia Patent Application No. 522440342 dated May 28, 2023, 24 pages with English translation. |
| International Preliminary Report on Patentability issued in PCT Application PCT/US2020/020339 dated Sep. 9, 2022, 9 pages. |
| Misra et al., 2016a. “Interfacial polarization of disseminated conductive minerals in absence of redox-active species—Part 1: Mechanistic model and validation.” Geophysics 81 (2): E139-E157. |
| Misra et al., 2016b. “Interfacial polarization of disseminated conductive minerals in absence of redox-active species—Part 2: Effective electrical conductivity and dielectric permittivity Interfacial polarization due to inclusions.” Geophysics 81 (2): E159-E176. |
| International Search Report and Written Opinion issued in PCT Application PCT/US2020/020339, dated Nov. 11, 2020 (12 pages). |
| Han et al., Hydrocarbon saturation in a Lower-Paleozoic organic-rich shale gas formation based on Markov-chain Monte Carlo stochastic inversion of broadband and electromagnetic dispersion logs, Fuel. May 1, 2019, vol. 234, No. 1, pp. 645-658. |
| Second Substantive Exam issued in Saudi Arabia Patent Application No. 522440342 dated Nov. 12, 2023, 9 pages with English translation. |
| Number | Date | Country | |
|---|---|---|---|
| 20230109189 A1 | Apr 2023 | US |
