Permeability anisotropy assessment in subsurface anisotropic formations

Information

  • Patent Grant
  • 11156738
  • Patent Number
    11,156,738
  • Date Filed
    Monday, January 25, 2016
    8 years ago
  • Date Issued
    Tuesday, October 26, 2021
    3 years ago
Abstract
Example embodiments are described for a method and system for logging data processing in determining permeability anisotropy effects. A permeability anisotropy model is used to derive a relationship between formation permeability anisotropy and resistivity anisotropy in both TI and BA formations. Implementations can provide the permeability anisotropy plus the true reservoir (or sand) permeability by using an integrated interpretation of the MCI resistivity anisotropy measurements with conventional permeability logs from other sensors (e.g., NMR or sonic). Biaxial and triaxial permeability components of the permeability anisotropy tensor can be determined for application to synthetic and field log interpretations.
Description
PRIORITY

The present application is a U.S. National Stage patent application of International Patent Application No. PCT/US2016/014679, filed on Jan. 25, 2016, the benefit of which is claimed and the disclosure of which is incorporated herein by reference in its entirety.


BACKGROUND

Modern operations for the exploration and production of oil and gas rely on access to a variety of information regarding subsurface geological parameters and conditions. Such information typically includes characteristics of Earth formations traversed by a borehole, as well as data relating to the size and mud of the borehole itself. The collection of information relating to subsurface conditions, which is commonly referred to as “logging,” can be performed by several methods, including wireline logging and logging while drilling (LWD).


In wireline logging, a sonde is lowered into the borehole after some or all of the well has been drilled. The sonde hangs at the end of a wireline cable that provides mechanical support to the sonde and also provides an electrical connection between the sonde and electrical equipment located at the surface. In accordance with existing logging techniques, various parameters of the Earth's formations are measured and correlated with the position of the sonde in the borehole as the sonde is pulled uphole. In LWD, a drilling assembly includes sensing instruments that measure various parameters as the formation is penetrated, thereby enabling measurement of the formation during the drilling operation.


Among the available wireline and LWD tools are a variety of logging tools including devices configured for taking multicomponent induction (MCI), nuclear magnetic resonance (NMR), and multipole sonic logging (MSL) measurements. Formation permeability can be characterized using such measurements. Current models can be used for the determination of permeability in isotropic formations. However, formation permeability can exhibit anisotropy in anisotropic formations, and is often dependent upon the measuring direction of logging tools.





BRIEF DESCRIPTION OF THE DRAWINGS

Some embodiments are illustrated by way of example and not limitation in the figures of the accompanying drawings in which:



FIG. 1 is a schematic view of a system for capturing subsurface measurement data in a logging while drilling operation, according to one or more example embodiments.



FIG. 2 is a schematic view of a system for capturing subsurface measurement data in a wireline logging operation, according to one or more example embodiments.



FIG. 3 is a schematic view of a coiled tubing logging system for capturing subsurface measurement data, according to one or more example embodiments.



FIG. 4 is a pair of plot diagrams illustrating simulated results for permeability components, kh and ky, and anisotropic ratio khv as a function of the resistivity anisotropy ratio, according to one or more example embodiments.



FIG. 5 is a diagram of a multi-component permeability model, according to one or more example embodiments.



FIG. 6 is a diagram illustrating a bimodal permeability model consisting of isotropic sand and transversely isotropic shale, according to one or more example embodiments.



FIG. 7 is a plot diagram of an example graphical solution for sand permeability ksd and laminar shale volume Vlam in laminar formations, according to one or more example embodiments.



FIG. 8 is plot diagram of another example graphical solution for sand permeability ksd and laminar shale volume Vlam in laminar formations, according to one or more example embodiments.



FIG. 9 is a pair of plots illustrating horizontal and vertical permeability, kh and kv, effective permeability ke, and anisotropic ratio khv=kh/kv as a function of laminar shale volume Vlam, according to one or more example embodiments.



FIG. 10 is a plot diagram illustrating a graphical solutions for ksd and ksh to at a laminar shale volume of Vlam=20%, according to one or more example embodiments.



FIG. 11 is a plot diagram illustrating simulated results for permeability anisotropy ratio khv as a function of resistivity anisotropy ratio Rvh, according to one or more example embodiments.



FIG. 12 is a pair of plot diagrams illustrating simulated results for triaxial permeability components, kx, ky, and kz, as a function of the resistivity anisotropy ratio, Rxy, according to one or more example embodiments.



FIG. 13 is flow chart illustrating a method for evaluation of formation permeability anisotropy, according to one or more example embodiments.



FIG. 14 is a set of plots that provides a synthetic data set without errors for predicting the permeability anisotropy ratio (khv) and the horizontal/vertical permeability components (kh and kv), according to one or more example embodiments.



FIG. 15 is a set of plots that provides a synthetic data set with a maximum random error of ±20% for predicting the permeability anisotropy ratio (khv) and the horizontal/vertical permeability components (kh and kv), according to one or more example embodiments.



FIG. 16 is a cross plot for kh vs kv vs Vlam, according to one or more example embodiments.



FIG. 17 is a schematic block diagram of a system for permeability determination, according to one or more example embodiments.



FIG. 18 is a diagrammatic representation of a machine in the example form of a computer system within which a set of instructions for causing the machine to perform a method for permeability determination may be executed, according to one or more example embodiments.





DETAILED DESCRIPTION

The following detailed description refers to the accompanying drawings that depict various details of examples selected to show how aspects of this disclosure may be practiced. The discussion addresses various examples of the disclosed subject matter at least partially in reference to these drawings, and describes the depicted embodiments in sufficient detail to enable those skilled in the art to practice the subject matter disclosed herein. Many other embodiments may be utilized for practicing the disclosed subject matter other than the illustrative examples discussed herein, and structural and operational changes in addition to the alternatives specifically discussed herein may be made without departing from the scope of the disclosed subject matter.


INTRODUCTION

Subsurface formations can exhibit the ability to allow formations fluids (e.g., oil, water, or gas) or multi-phase fluids to pass through it, referred to as formation permeability. Formation permeability is an important parameter in formation evaluation and reservoir characterization. For example, permeability information can be utilized for reservoir simulation, enhanced oil recovery, well completion design, and field exploitation/development strategies. Unlike isotropic petro-physical parameters such as porosity and saturation, formation permeability can be anisotropic in, for example, three common types of subsurface formations (e.g., laminated sand-shale sequences, sands with different grain sizes, and sand with thin resistive/conductive streaks). Thus, permeability is often strongly dependent on the measuring direction in the anisotropic formations.


Various interpretation models for deriving formation permeability from logging data (e.g., multicomponent induction [MCI], nuclear magnetic resonance [NMR], and multipole sonic logging [MSL]) have been developed. However, such models can typically only be used for determining permeability in isotropic formations. Current logging tools and their measuring data can only be used for determination of scalar or isotropic permeability; however, these conventional log-derived permeabilities cannot deliver permeability anisotropy and its components (such as horizontal and vertical permeability). Accurate determination of permeability that accounts for anisotropy can enhance formation evaluation and reservoir characterization.


Permeability can be determined using a permeability anisotropy model that is based at least in part on relationship between formation permeability and resistivity in transversely isotropic (TI) and biaxially anisotropic (BA) formations. In MCI logging, formation resistivity (or conductivity, which is inversely related to resistivity) shows azimuthal anisotropy of horizontal resistivity in the bedding plane of formations. Unless the text or context clearly indicates otherwise, “horizontal” or “transverse” means a direction or plane substantially coinciding with a bedding plane of the relevant formation, and “vertical” means a direction of plane substantially orthogonal to the bedding plane of the relevant formation.


Processing of MCI logging data may be performed based at least in part on BA parameterization, resulting in more accurate description of complex anisotropic formations than is the case based on TI models. The TI formation model is a model that represents simulated formation resistivity characteristics accounting for transverse formation isotropy to resistivity. The TI model can account for resistivity differences between orthogonal axes lying in a formation or bedding plane (e.g., the horizontal or transverse plane) and an axis perpendicular to the formation or bedding plane (e.g., the vertical axis). Thus, the TI model accounts for anisotropy between the vertical axis and the horizontal plane, but assumes isotropy between different axes in the horizontal plane. The BA model additionally accounts for anisotropy between orthogonal axes in the transverse plane.


One or more example embodiments described below provide a method and system for logging data processing in determining permeability anisotropy. The description of example embodiments that follow describe using the permeability anisotropy model to derive a relationship between formation permeability anisotropy and resistivity anisotropy in both TI and BA formations. As is known, MCI tools are capable of measuring resistivity anisotropy in formations. An effective permeability or a permeability component can be used to calibrate conventional log-derived permeabilities (e.g., that only account for scalar or isotropic permeability) for calculating the other permeability components. The limited vertical resolution of logging tools are sometimes unable to capture the reservoir (or sand) permeability in laminar formations. Thus, in some embodiments, the sand permeability is determined from calculating horizontal and vertical permeability components based on a multi-modal permeability model, as discussed further below.


Implementations of the disclosed example embodiments can provide the permeability anisotropy by using an integrated interpretation of the MCI resistivity anisotropy measurements with conventional permeability logs from other sensors (e.g., resistivity, NMR, or sonic). Thus, biaxial (for TI formations) or triaxial (for BA formations) permeability components of the permeability anisotropy tensor can be determined and applied to field log interpretations.


Example Measurement Environments



FIG. 1 is a schematic illustration of an example logging while drilling (LWD) environment. A drilling platform 102 is equipped with a derrick 104 that supports a hoist 106 for raising and lowering a drill string 108. The hoist 106 suspends a top drive 110 suitable for rotating the drill string 108 and lowering the drill string 108 through the well head 112. Connected to the lower end of the drill string 108 is a drill bit 114. As the drill bit 114 rotates, it creates a borehole 116 that passes through various formations 118. A pump 120 circulates drilling fluid through a supply pipe 122 to top drive 110, down through the interior of drill string 108, through orifices in drill bit 114, back to the surface via an annulus around drill string 108, and into a retention pit 124. The drilling fluid transports cuttings from the borehole 116 into the pit 124 and aids in maintaining the integrity of the borehole 116. Various materials can be used for drilling fluid, including a salt-water based conductive mud.


An assembly of LWD tools 126 is integrated into a bottom-hole assembly (BHA) near the bit 114. As the bit 114 extends the borehole 116 through the formations 118, LWD tools 126 collect measurements relating to various formation properties as well as the tool orientation and various other drilling conditions. The LWD tools 126 may take the form of a drill collar, i.e., a thick-walled tubular that provides weight and rigidity to aid the drilling process. In various example embodiments, the LWD tools 126 include a multi-array triaxial induction tool to measure formation resistivity and deliver MCI measurement data (also referred to as “log data”), such as described further herein. Further, the LWD tools 126 can include a NMR and/or a sonic logging tool to measure and deliver logging data. A telemetry sub 128 may be included to transfer images and measurement data to a surface receiver 200 and to receive commands from the surface. In some embodiments, the telemetry sub 128 does not communicate with the surface, but rather stores logging data for later retrieval at the surface when the logging assembly is recovered.


At various times during (or after) the drilling process, the drill string 108 may be removed from the borehole 116 as shown in FIG. 2. Once the drill string 108 has been removed, logging operations can be conducted using a wireline logging sonde 234, i.e., a probe suspended by a cable 242 having conductors for conducting power to the sonde 234, and for transmitting telemetry data from the sonde 234 to the surface. The example wireline logging sonde 234 may have pads and/or centralizing springs to maintain the sonde 234 near the central axis of the borehole 116 as sonde 234 is pulled uphole. The logging sonde 234 can include a variety of sensors including a multi-array triaxial induction tool for measuring formation resistivity and providing MCI measurement data. Further, the logging sonde 234 can include a NMR and/or sonic logging tool for measuring and providing measurement data. A logging facility 244 collects measurements from the logging sonde 234, and includes a computer system 245 for processing and storing the measurements gathered by the sensors.


Yet a further alternative logging technique is schematically illustrated in FIG. 3, which shows an example embodiment of a coiled tubing logging system 300. In system 300, coiled tubing 354 is pulled from a spool 352 by a tubing injector 356 and injected through a packer 358 and a blowout preventer 360 into the borehole 116. In the borehole 116, a supervisory sub 364 and one or more logging tools 365 are coupled to the coiled tubing 354 and configured to communicate to a surface computer system 366 via information conduits or other telemetry channels. An uphole interface 367 may be provided to exchange communications with the supervisory sub 364 and receive data to be conveyed to the surface computer system 366.


Surface computer system 366 is configured to communicate with supervisory sub 364 to set logging parameters and collect logging information from the one or more logging tools 365. Surface computer system 366 is configured by software (shown in FIG. 3 as being stored on example embodiments of removable storage media 372) to monitor and control downhole instruments 364, 365. The surface computer system 366 may be a computer system such as that described further herein.


Forward Models and Their Libraries


Forward modeling includes a numerical solution of Maxwell's equation in a mathematical boundary value problem, where the relevant formation or model specifies boundaries and shapes of regions of different resistivity. Processes for deriving formation parameters from a set of given field logs is known as inverse modeling, and typically comprises iteratively adjusting selected formation parameters in one or more layers of a formation model, and repeating forward modeling (e.g., by and the calculation or with reference to pre-calculated library data), until the observed field logs are satisfactorily replicated, e.g. until a set of variable borehole-formation parameters is found for a best fit with the observed log data based on the applicable formation model.


Some of the example processing schemes disclosed herein are based at least in part on a biaxial anisotropy (BA) model (otherwise known as orthorhombic anisotropy) and/or a transversely isotropic (TI) model. The TI model can account for resistivity differences between, on the one hand, orthogonal axes lying in a formation or bedding plane (sometimes referred to as the horizontal or transverse plane), and, on the other hand, an axis perpendicular to the formation or bedding plane (sometimes referred to as the vertical axis). The TI model thus can account for anisotropy between the “vertical” axis and the “horizontal” plane, but assumes isotropy between different axes in the “horizontal” or transverse plane. For this reason, the TI model is also referred to being TI anisotropic.


The BA model additionally accounts for anisotropy between orthogonal axes in the transverse plane, and is therefore also referred to as accounting for triaxial anisotropy. Note that, unless otherwise specified, “biaxial anisotropy” and its derivations refer to transverse biaxial anisotropy. Consistent with this terminology, a TI model does not account for biaxial anisotropy, even though it accounts for anisotropy between two axes (e.g., between the horizontal plane and the vertical axis).


We describe formation resistivity/conductivity of BA anisotropy in the formation's principal coordinate system. This system is chosen so that the x- or y-axis direction coincides with a conductivity tensor's principal axis having the largest conductivity (or smallest resistivity) component in the bedding plane; the z axis is parallel to the conductivity principal axis having the smallest conductivity component. In this principal axis coordinate system, the formation conductivity can be expressed as a diagonal tensor:

σt=diag(Cx,Cy,Cz).  (1)

where σ is the formation conductivity tensor, its elements Cx and Cy are the two conductivities in the two principal-axis (e.g., x- and y-axes) in the bedding plane, and Cz is the conductivity in the principal-axis direction perpendicular to the bedding plane.


If equation (1) is expressed in resistivity terms, the formation resistivity can be expressed as the following diagonal tensor:

Rt=diag(Rx,Ry,Rz),  (2)

where Rt is the formation resistivity tensors, and its elements Rx, Ry, and Rz are the triaxial resistivity components in the three principal-axis directions (xf, yf, and zf), respectively. It is noted that in a TI formation, resistivity can be represented as a diagonal tensor that is only described using two resistivity components: Rh (wherein Rh=Rx=Ry) and Rv (wherein Rv=Rz) in the principal axis system (e.g., x-y-z coordinate system), however, the resistivity tensor of the BA model is best described using the triaxial resistivity components: Rx, Ry, and Rz.


Based on equations (1) and (2), one can derive the relations among resistivity and conductivity components Rx=1/Cx, Ry=1/Cy, and Rz=1/Cz. If Rx=Ry=Rz, then the formation resistivity is isotropic. If only Rx=Ry≠Rz, then the formation resistivity is transversely isotropic; therefore, only one resistivity component is necessary in the bedding plane. In this case, both Rx and Ry are usually called horizontal resistivity and are often denoted as Rh (=Rx=Ry); and, Rz is called the vertical resistivity and is usually denoted as Rv. If Rx≠Ry≠Rz, the formation resistivity is of BA anisotropy. It is can be seen that the isotropy and transverse isotropy are only two special cases of the biaxial anisotropy. For the practical applications, different ratio notations are used. For example, ratios Rxy=Rx/Ry, or Rzx=Rz/Rx and Rzy=Rz/Ry.


Permeability Model in Anisotropic Formations


As discussed herein, when permeability is anisotropic, a permeability tensor (as opposed to a scalar) is used for expressing the formation permeability. Pressure in formations can be applied in three directions, and for each direction, the permeability can be measured (via Darcy's law that describes the flow of a fluid through a porous medium) in three directions, thus leading to a 3 by 3 tensor. In a three-dimensional (3D) coordinate system, the tensor is realized using a 3 by 3 matrix being both symmetric and positive definite, and can be represented using the following equation:










K
=


(




k
xx




k
xy




k
xz






k
xy




k
yy




k
yz






k
xz




k
yz




k
zz




)

=


(

k
ij

)


3
×
3




.




(
3
)







In equation (3), K is the true formation permeability tensor, i denotes the pressure direction, and j is the permeability measurement direction. Therefore, component kij represents the j-th directional permeability component at the i-th pressure direction. This permeability tensor K is diagonalizable (as it is both symmetric and positive definite). The eigenvectors will yield the principal directions of flow, representing the directions where flow is parallel to the pressure drop, and the three eigenvalues represent the three principal permeability components.


In the 3D principal coordinate system, the permeability tensor K is expressed as a diagonal tensor after the matrix diagonalization, which is equivalent to finding the matrix's eigenvalues, and can be represented using the following equation:

K=diag(kx,ky,kz)  (4)


In equation (4), kx and ky are the two permeability components in the two principal axis (X and Y) directions of the bedding plane, and kz is the permeability component in the Z principal axis direction perpendicular to the bedding plane. Following the formation resistivity descriptions, the formation permeability is isotropic if kx=ky=kz. If kx=ky≠kz, then the permeability is transversely isotropic. In the TI formation, only two components are needed for describing the formation permeability anisotropy. In the bedding plane, the permeability component can be denoted as kx=ky=kh and is referred to as the horizontal permeability. For the Z-directional component, kz, can be denoted as kz=kv and is referred to as the vertical permeability. Thus, the permeability tensor can also be represented using the equation: K=diag(kh, kh, kv). When kx≠ky≠kz, the permeability is of biaxial anisotropy (BA).


In electrically anisotropic formations, the resistivity anisotropy can be obtained from the 3D induction resistivity logs, such as MCI logs. If formations are both permeability and resistivity anisotropic, the permeability anisotropy can be evaluated from the resistivity anisotropy. It is assumed that the pore space is represented as a bundle of independent and tortuous tubes of different radii in isotropic formations. If the flow rate is low enough that it is laminar and not turbulent, then the permeability scalar can be determined by using the following Kozeny-Carman equation:









k
=


A
·


d
2


4

F



=

A
·

d
2

·


I
·

R
w



4


R
t









(
5
)








where k represents the permeability in an isotropic formation, A represents a shape factor for the pore tubes with a diameter of d, and F represents the formation factor.


The formation factor can be represented using the following equation:










F
=



R
o


R
w


=


R
t


I
·

R
w





,




(
6
)








where Ro represents the 100% water-bearing formation resistivity, Rw represents the formation water resistivity, Rt represents the true formation resistivity, and I represents the resistivity index.


The resistivity index can be represented using the following equation:









I
=


R
t


R
o






(
7
)







In TI anisotropic formations, the permeability, formation factor, and resistivity (or conductivity) index are not scalars. Instead, the permeability, formation factor, and resistivity (or conductivity) index are represented using tensors, namely: K, F, and I, respectively. F and I can be represented using the following two diagonal tensor equations:

F=diag(Fh,Fh,Fv)  (8)
I=diag(Ih,Ih,Iv)  (9)


A generalization of all the values to anisotropic formation yields the tensor expressions. For example, the following equations represent permeability scalars in TI formations:










k
h

=



A
h

·

d
h
2

·

1

4


F
h




=


A
h

·

d
h
2

·



I
h

·

R
w



4


R
h









(
10
)







k
v

=



A
v

·

d
v
2

·

1

4


F
v




=


A
v

·

d
v
2

·



I
v

·

R
w



4


R
v









(
11
)








where Fh and Fv represent the formation factor components measured along the horizontal and vertical direction, respectively, which can be represented by







F
h

=




R
h



I
h

·

R
w








and






F
v


=



R
v



I
v

·

R
w



.







The two components of the resistivity index tensor, I, can be represented by








I
h

=


R
h


R
O
h



,


and






I
v


=


R
v


R
O
v



,





where Rh and Rv represent the formation horizontal and vertical resistivity, respectively.


Relationship Between Permeability Anisotropy and Resistivity Anisotropy


In TI permeability-anisotropic formations, the relationship between permeability and formation resistivity factor obtained along two perpendicular principal-axis directions can be represented using the following equation:











k
h


k
v


=




A
h

·

d
h
2




A
v

·

d
v
2



·

(


F
v


F
h


)






(
12
)







Alternatively, the relationship between permeability and anisotropy and resistivity (or conductivity) anisotropy can be represented using the following equation:











k
h


k
v


=





A
h

·

d
h
2




A
v

·

d
v
2







I
h


I
v


·

(


R
v


R
h


)



=





A
h

·

d
h
2




A
v

·
d






I
h


I
v


·

(


C
h


C
v


)



=


C

R

P




(


R
v


R
h


)








(
13
)








where








c

R

P


=




A
h

·

d
h
2




A
v

·

d
v
2






I
h


I
v




,





with Ch and Cv representing the horizontal and vertical conductivity, respectively. A resistivity (or conductivity) anisotropy ratio can be defined as:








R

v

h


=



R
v


R
h


=



C
h


C
v


=

C
vh




,





and a permeability anisotropy ratio khv can be defined as:







k

h

v


=



k
h


k
v


.





This permeability-resistivity anisotropy relationship shows that the permeability anisotropy can be evaluated based on resistivity anisotropy measurements (such as from MCI logs) and the pore shape/diameter. A simplified relationship between the permeability and resistivity anisotropy of equation (13) can be determined if it is assumed that AhIh=AvIv, and can be represented using the following equation:












k
h


k
v


=



d
h
2


d
v
2


·

(


R
v


R
h


)



,


or






k

h

v



=


C
RP

·

R

v

h




,


or






R

v

h



=


1

C

R

P



·

k

h

v








(
14
)







Further, if it is assumed that








C

R

P


=



d
h
2


d
v
2


=
1


,





then equation (14) can further be reduced and represented using the following equation:











k
h


k
v


=

(


R
v


R
h


)





(
15
)








From equation (15), the permeability anisotropy can be estimated if the resistivity anisotropy is known, such as from MCI log inversion processing.


Determination of Permeability Tensor Components in TI Formations


An effective (or geometric mean) permeability as measured by a conventional formation logging tool can be represented using the following equation in TI formations:










k
e

=



k
h
2

·

k
v


3





(
16
)







The effective permeability, ke, should fall between kh and kv. It is noted that kh=kv=ke for isotropic media. For core measurements, it is possible to measure one permeability component or multi-component permeability, such as kh or kv or both components. In one case, if it is assumed that the log-derived permeabilities (e.g., using resistivity logs, sonic and/or NMR) are calibrated by ke, then they are approximately equal to this effective permeability, ke. In another case, if it is assumed that the log-derived permeabilities are calibrated by kh or kv, then they are approximately equal to the horizontal or vertical permeability, kh or kv.


For example, if it is assumed that the resistivity anisotropy can be obtained from the MCI resistivity measurements, then the horizontal and vertical permeabilities are determined from equations (15) and (16) if the log-derived permeabilities are equal to this effective permeability, ke. Thus, the following equations can be used for determining kh and kv:










k
h

=



k
e

·


k

h

v


3


=


k
e

·


R

v

h


3







(
17
)







k
v

=



k
e


k
hv

2
/
3



=


k
e


R
vh

2
/
3








(
18
)








FIG. 4 is a pair of plot diagrams illustrating simulated results for permeability components, kh and kv, and anisotropic ratio khv as a function of the resistivity anisotropy ratio, Rvh using equations (17) and (18). The effective permeability is not equal to the horizontal and vertical permeabilities if Rvh>1. Both kh and kv are two log-linear functions of Rvh. In the left plot, the effective permeability is 50 md. In the right plot, the effective permeability is 500 md. Therefore, after kh and kv are available, the permeability tensor K=diag(kh, kh, kv) can be determined.


True Permeability Determination in Laminar Reservoirs


It is known that the laminations in laminar reservoirs are often much thinner than a logging tool's vertical resolution. Therefore, the log-derived permeabilities can be represented as weighted averages of true reservoir permeability and shale permeability. FIG. 5 is a diagram of a multi-component permeability model consisting of M-number types of isotropic sands and one TI shale. In the model, ksd(i) is the i-th type true or sand permeability, and kshh and kshv represent the horizontal and vertical permeabilities of the pure shale formation. According to the multi-component permeability model, the following equations can be used for determination of true reservoir permeability in a laminar formation, such as sand-shale formations:


For permeability in a direction parallel (e.g., horizontal) to the bedding plane, the horizontal permeability can be represented using the following equation:

khi=1M(Vsd(i)·ksd(i))+kshh·Σi=1MVlam(i)i=1M(Vsd(i)·ksd(i))+kshh·Vlam  (19)


For permeability in a direction perpendicular (e.g., vertical) to the bedding plane, the vertical permeability can be represented using the following equation:










1

k
v


=






i
=
1

M



(


V
sd

(
i
)



k
sd

(
i
)



)


+



Σ

i
=
1

M



V
lam

(
i
)




k

s

h

v



=





i
=
1

M



(


V

s

d


(
i
)



k

s

d


(
i
)



)


+


V
lam


k

s

h

v








(
20
)







In equations (19) and (20), it is assumed that all sands are isotropic and the shale is TI anisotropic. M represents the total number of the sand types. kshh and kshv represent the horizontal and vertical permeabilities of the pure shale formation, respectively. Vlami=1M Vlam(i) represents the total laminated shale volume fraction, Vsd(i) represents the volume percentage for the i-th type sand or reservoir, ksd(i) represents the i-th type sand permeability, and Σi=1M (Vsd(i))+Vlam=1.


If it is assumed that M=1 or ksd(1)=k=sd(2)= . . . =ksd(M)=ksd, the above equations can be reduced to a bimodal permeability model for determination of reservoir permeability. FIG. 6 is a diagram illustrating a bimodal permeability model consisting of isotropic sand and TI shale. Here, ksd is the true sand (or reservoir) permeability, and kshh and kshv are the horizontal and vertical permeabilities of the pure shale formation. As illustrated, a multi-layer model can be reduced to an equivalent two-layer bimodal model.


The horizontal permeability component can be represented using the following reduced equation:

kh=(1−Vlamksd+Vlam·kshh  (21)


The vertical permeability component can be represented using the following reduced equation:











1

k
v


=



(

1
-

V
lam


)


k

s

d



+


V
lam


k

s

h

v




,


or






k
v


=

1



(

1
-

V
lam


)


k
sd


+


V
lam


k

s

h

v









(
22
)







In equations (21) and (22), if kh, kv, kshh and kshv are known, then the sand permeability ksd and Vlam can be solved. For example, the graphical solution of sand permeability ksd and laminar shale volume Vlam to equations (21) and (22) for laminar shale at kshh=10 md and kshv=10 md are shown in FIG. 7. In FIG. 7, the dashed-line contours correspond to values of constant Vlam, and the solid-line contours represent values of constant ksd. Similarly, FIG. 8 illustrates a graphical solution for kshh=10 md and kshv=5 md. In FIG. 8, the dashed-line contours correspond to values of constant Vlam, and the solid-line contours represent values of constant ksd.


However, if it is assumed that both sand and shale are isotropic (e.g., kshh=kshv=ksh), equations (21) and (22) can be expressed using the following equations:










k
h

=



(

1
-

V
lam


)

·

k
sd


+


V
lam

·

k

s

h








(
23
)








1

k
v


=



(

1
-

V
lam


)


k

s

d



+


V
lam


k

s

h





,


or






k
v


=

1



(

1
-

V
lam


)


k
sd


+


V
lam


k
sh









(
24
)







However, even though it is assumed that both the sand and shale are isotropic, the macro anisotropy of formation permeability can observed. For example, FIG. 9 is a pair of plots illustrating horizontal and vertical permeability, kh and kv, effective permeability ke, and anisotropic ratio khv=kh/kv as a function of laminar shale volume Vlam. If Vlam≠0, ksd and ksh represent the sandstone and shale permeability, respectively. If it is assumed that the two permeability components, kh, kv, and laminar shale volume Vlam are known, then ksd and ksh can be resolved from equations (23) and (24). After the lamination effects on the measured/calculated permeabilities, kh, kv, are corrected, the more accurate permeability ksd for sand can be obtained. In FIG. 9, the sandstone permeability ksd=100 md and kshh=10 md. The left plot has shale permeability kshHV=kshv=1 or the shale is isotropic. The right plot has shale permeability kshHV=kshv=5 or the shale is TI anisotropic. These plots show that the khv approaches the maximum value at approximately 50% shale by volume. FIG. 10 is a plot diagram illustrating a graphical solutions for ksd and ksh to Equations (23) and (24) at a laminar shale volume of Vlam=20%. In FIG. 10, the solid-line contours correspond to values of constant ksh, and the dashed-line contours represent values of constant ksd.


Using equations (23) and (24), the permeability anisotropy ratio







k

h

v


=


k
h


k
v







can be represented using the following equation:










k
hv

=



k
h


k
v


=




[



(

1
-

V
lam


)

·

k

s

d



+


V
lam

·

k

s

h




]




(

1
-

V
lam


)


k

s

d




+


V
lam


k

s

h




=

1
+



V
lam



(

1
-

V
lam


)


·

[



k

s

d



k

s

h



+


k

s

h



k

s

d



-
2

]









(
25
)







Further, the following equations can be used for calculating Rh and Rv:










1

R
h


=



(

1
-

V
lam


)


R

s

d



+


V
lam


R

s

h

h







(
26
)







R
v

=



(

1
-

V
lam


)

·

R
sd


+


V
lam

·

R

s

h

v







(
27
)







If Rshh=Rshv=Rsh, the following resistivity anisotropy ratio occurs:










R

v

h


=



R
v


R
h


=

1
+



V
lam



(

1
-

V
lam


)


·

[



R

s

d



R

s

h



+


R

s

h



R

s

d



-
2

]








(
28
)







From the above, it can be observed that: if









k

s

d



k

s

h



=


R

s

d



R

s

h




,





then khv=Rvh (or CRP=1). Otherwise, khv=CRP·Rvh, and CRP is a non-linear function of Vlam, ksd/ksh, and Rsd/Rsh, and can be represented using the following equation:










C

R

P


=


1
+



V
lam



(

1
-

V
lam


)


·

[



k

s

d



k

s

h



+


k

s

h



k

s

d



-
2

]




1
+



V
lam



(

1
-

V
lam


)


·

[



R

s

d



R

s

h



+


R

s

h



R

s

d



-
2

]








(
29
)







Once all of Vlam,








k

s

d



k

s

h



,

and







R

s

d



R

s

h









are known, equation (29) can be used to estimate CRP. FIG. 11 is a plot diagram illustrating simulated results for permeability anisotropy ratio khv as a function of resistivity anisotropy ratio Rvh using equations (25) and (28). In FIG. 11,









k

s

d



k

s

h



=
100

,





the dashed lines above the solid line have values of:










R

s

d



R

s

h



=
1

,
5

,
10
,
20
,





and 50 from left to right; blue solid line has a value of:









R

s

d



R

s

h



=
100





;





the dashed lines below the solid line have values of:









R

s

d



R

s

h



=

2

0

0


,

3

0

0

,

5

00

,





and 1000 from left to right; and Vlam=0−0.8. FIG. 11 provides a graphical solution to CRP, and illustrates a log-linear correlation between khv and Rvh if both ratios are less than 10, or CRP is not sensitive to Vlam.


Permeability Anisotropy Determination in BA Formations


Equations similar to those previously discussed with regard to equations (15) and (16) for TI formations can be inferred for BA formations and represented using the following equations:










k

x

z


=



k
x


k
z


=


(


R
z


R
x


)

=

R
zx







(
30
)







k

y

z


=



k
y


k
z


=


(


R
z


R
y


)

=

R
zy







(
31
)







k
yx

=



k
y


k
x


=


(


R
x


R
y


)

=

R

x

y








(
32
)







k
e

B

A


=



k
x

·

k
y

·

k
z


3





(
33
)








where Rzy=Rzx·Rxy. Once the triaxial resistivity components (Rx, Ry, and Rz) and the effective permeability keBA or one permeability component such as kx are available, then the tri-axial permeability components (kx, ky, and kz) and permeability anisotropy






(


e
.
g
.

,


k
x


k
y


,

and







k
y


k
z




)





can be obtained.


For example, if keBA and the resistivity anisotropy ratios are known, then equation (33) can be solved for the following:










k
x

=


k
e

B

A


·


(


R

z

x



R

x

y



)


1
/
3







(
34
)







k
y

=


k
e

B

A


·


(


R

x

y






·

R

z

y



)


1
/
3







(
35
)







k
z

=


k
e

B

A




(


R

x

y


·

R

z

y



)


1
/
3







(
36
)








FIG. 12 is a pair of plot diagrams illustrating simulated results for triaxial permeability components, kx, ky, and kz, as a function of the resistivity anisotropy ratio, Rxy, by using equations (34)-(36). In FIG. 12, Rzy=6, and two effective permeability plots are displayed. As illustrated, kx, ky, and kz, are three log-linear functions of horizontal resistivity anisotropy ratio, Rxy. The left plot diagram has effective permeability ke=50 md and the right plot diagram has effective permeability ke=500 md. It is noted that the plot diagrams of FIG. 12 (e.g., for a biaxial anisotropic formation) is similar to those of FIG. 4 (e.g., for a TI anisotropic formation).


Workflow for Prediction of Permeability Anisotropy



FIG. 13 is a flow chart illustrating a method 1300 comprising multi-level data processing based on multiple forward models with BA and TI anisotropies. The method 1300 provides for evaluation of formation permeability anisotropy using an integration of MCI log data with other sensor logs (e.g., NMR, sonic, or other logging tools).


At operation 1402, MCI measurement data captured by a triaxial MCI tool in a borehole extending through a subsurface geological formations is inputted after calibration, temperature correction and other preprocessing. This preprocessing does not include skin-effect correction. The MCI data may be multi-frequency data, and may be taken at multiple spacings. In some embodiments, the MCI measurement data can be single-frequency measurements for the respective arrays of the tool. Further, other log data captured by, for example, sonic or NMR logging tools are also inputted.


At operation 1304, inversion processing is applied to MCI measurement data to produce inverted parameters to be used for permeability evaluation. For example, a TI-based inversion can be applied to produce the inverted TI parameters of Rh and Rv. Further, a BA-based inversion can be applied to produce the inverted BA parameters of Rx, Ry, and Rz. Each of these TI- and BA-based inversion processings can be based on various formation models, including, but not limited to: radially one-dimensional (R1D) and zero-dimensional (0D) models. The inverted parameters can be used for the determination of, for example, the resistivity anisotropy (Rvh). The permeability anisotropy ratio (khv) can be determined from the resistivity anisotropy using, for example, equations (12)-(15) as discussed above. Operation 1304 can further include determining the laminar shale volume Vlam.


At operation 1306, the effective permeability ke (or a permeability component) can be determined using conventional log-derived permeabilities (e.g., using NMR, sonic, or resistivity logs). At operation 1308, permeability components can be determined using resistivity anisotropy data obtained from MCI resistivity measurement data and the effective permeability ke from conventional permeability logs. In one embodiment, equations (17)-(18) as discussed above can be used for determining the permeability components, kh and kv, in TI formations. The permeability components can be used to solve equations (19)-(20) as discussed above for recovering the horizontal and vertical permeabilities kshh and kshv of a pure shale formation. With a known laminar shale volume Vlam, equations (21) and (22) as discussed above can be solved to obtain the sand permeability ksd in laminated formations.


In other embodiments, resistivity and conventional permeability log data can be used for determining the permeability components in BA formations. The tri-axial resistivity components (Rx, Ry, and Rz) and the effective permeability (or a permeability component), from the previous operations, can be used to determine the tri-axial permeability components (kx, ky, and kz) and permeability anisotropy






(


e
.
g
.

,


k
x


k
y


,

and







k
y


k
z




)





using, for example, equations (34)-(36) as discussed above.


Once the above data is determined, equation (29) as discussed above can be used to estimate the CRP coefficient. In some embodiments, The resulting data, the calculated kh and kv or kx, ky, and kz, and ksd (and ksh) can be evaluated for data quality at operation 1310 before output for use in formation evaluation at 1312.


Benefits of the described methods and systems for permeability anisotropy assessment using both MCI data and conventional permeability logs include more accurate reservoir assessment, fracture detection, and oil development/production.


Workflow Validation



FIG. 14 is a set of plots that provides a synthetic data set without errors for predicting the permeability anisotropy ratio (khv) and the horizontal and vertical permeability components (kh and kv) from an integrated interpretation of Vlam, Rv/Rh, and the effective permeability (ke).


Track 1, the plot illustrated furthest to the left of FIG. 14, illustrates simulated horizontal and vertical resistivity (Rh and Rv having units of ohm-m) over a range of depths. The resistivity anisotropy (Rvh) is also illustrated in Track 1. Rh and Rv are calculated using equations (26) and (27), as discussed above, for a given Rsd=10 ohm-m, Rshh=Rshv=2 ohm-m, and the known laminar shale volume log Vlam is plotted in the adjacent Track 2. The laminar shale volume log Vlam is simulated using the below equation:

Vlam(z,AA,BB)=AA−BB·Σi=1Nshai·gi(z,bi,ci),  (37)

where z represents the logging depth (in units of feet), AA=0.95, BB=0.85, and Nsh=7.


gi(z, bi, ci) represents multiple, modified Gaussian functions, which can be expressed using the below equation:












g
i



(

z
,

b
i

,

c
i


)


=

e


-

1
2





(


z
-

b
i



c
i


)

2




,

i
=
1

,
2
,





,

N

s

h






(
38
)








where ai=1, bi=−70, −50, −25, 0, 25, 50, 70, and ci=5, 5, 10, 5, 10, 5, 5.


Track 3 illustrates simulated horizontal/vertical permeability (kh/kv), effective permeability (ke), and the permeability anisotropy ratio (khv=kh/kv) as determined by equations (23), (24), and (18) with a given ksd=100 md and kshh=kshv=20 md, and further with the known Vlam from Track 2. Track 4 illustrates the predicted horizontal/vertical permeability (kh/kv), sand permeability (ksd), and the permeability anisotropy ratio (khv=kh/kv) using the workflow described in FIG. 13.


Due to the permeability macro-anisotropy caused by laminations and the limitation of tool vertical resolution, the effective permeability ke from conventional logs is not the same as kh and kv. Further, it can be seen that all of ke, kh and kv are significantly different from the ksd log for the same lamination reason. However, by comparing Tracks 3 and 4, it can be seen that the predicted kh & kv, ksd, and khv in Track 4 are all identical to their true values in Track 3.



FIG. 15 is a set of plots that provides a synthetic data set with a maximum random error of ±20% for predicting the permeability anisotropy ratio (khv) and the horizontal/vertical permeability components (kh and kv) from an integrated interpretation of Vlam, Rv/Rh, and the effective permeability (ke). The laminar shale volume log Vlam is simulated with random error using the below equation:

Vlamerr(z,AA,BB)=(1+0.2*[2*randomErr−1])*Vlam(z,AA,BB)  (38)

where randomErr is a random error generated by computer (e.g., 0.0≤randomErr≤1.0). From the above, it is known that the max error is 20%.


Track 1, the plot illustrated furthest to the left of FIG. 15, illustrates simulated horizontal and vertical resistivity (Rh and Rv having units of ohm-m) over a range of depths. The resistivity anisotropy (Rvh) is also illustrated in Track 1. Rh and Rv are calculated using equations (26) and (27), as discussed above, for a given Rsd=10 ohm-m, Rshh=Rshv=2 ohm-m, and the known laminar shale volume log Vlam is plotted in the adjacent Track 2. Track 3 illustrates simulated horizontal/vertical permeability (kh/kv), effective permeability (ke), and the permeability anisotropy ratio (khv=kh/kv) as determined by equations (23), (24), and (18) with a given ksd=100 md and kshh=kshv=20 md, and further with the known Vlam from Track 2. Track 4 illustrates the predicted horizontal/vertical permeability (kh/kv), sand permeability (ksd), and the permeability anisotropy ratio (khv=kh/kv) using the workflow described in FIG. 13. Even though the data set of FIG. 15 includes error (as opposed to FIG. 14), it can be seen that the predicted permeability data obtained in Track 4 still compares well relative to the true values shown in Track 3.


Referring now to FIG. 16, a cross plot is illustrated for kh vs kv vs Vlam for the data set shown in FIG. 15. From the cross-plot the shale permeability ksh can be determined using the highest Vlam. It can be observed that the shale is isotropic and has a permeability of ksh=20 md.


Example System



FIG. 17 is a schematic block diagram of an example system 1700 for estimating subsurface formation and invasion properties, according to an example embodiment. The example system 1700 of FIG. 17 may be configured to perform one or more of the methods described above with reference to FIG. 13. The system 1700 is described in terms of a number of modules for performing the respective operations previously described. As used herein a “module” can be an engine, logic component, or mechanism capable of performing described operations and/or configured or arranged in a certain manner. Modules may constitute either software modules, with code embodied on a non-transitory machine-readable medium (i.e., such as any conventional storage device, such as volatile or non-volatile memory, disk drives or solid state storage devices (SSDs), etc.), or hardware-implemented modules. In certain example embodiments, one or more computer systems (e.g., a standalone, client, or server computer system) or one or more components of a computer system (e.g., a processor or a group of processors) may be configured by software (e.g., an application or application portion) or firmware (note that software and firmware can generally be used interchangeably herein as is known by a skilled artisan) as a module that operates to perform the described operations.


In various embodiments, a hardware-implemented module may be implemented mechanically or electronically. A hardware-implemented module is a tangible unit capable of performing certain operations and may be configured or arranged in a certain manner. In example embodiments, one or more computer systems (e.g., a standalone, client, or server computer system) or one or more processors may be configured by software (e.g., an application or application portion) as a hardware-implemented module that operates to perform certain operations as described herein. For example, a hardware-implemented module may include dedicated circuitry or logic that is permanently configured (e.g., within a special-purpose processor, application specific integrated circuit (ASIC), or logic array) to perform the identified operations. A hardware-implemented module may also include programmable logic or circuitry (e.g., as encompassed within a general-purpose processor or other programmable processor) that is temporarily configured by software or firmware to perform some or all of such operations.


The term “hardware-implemented module” should be understood to encompass a tangible entity, be that an entity that is physically constructed, permanently configured (e.g., hardwired), or temporarily or transitorily configured (e.g., programmed) to operate in a certain manner and/or to perform certain operations described herein. Considering embodiments in which hardware-implemented modules are temporarily configured (e.g., programmed), each of the hardware-implemented modules need not be configured or instantiated at any one instance in time. For example, where the hardware-implemented modules comprise a general-purpose processor configured using software, the general-purpose processor may be configured as respective different hardware-implemented modules at different times. Software may accordingly configure a processor, for example, to constitute a particular hardware-implemented module at one instance of time and to constitute a different hardware-implemented module at a different instance of time.


Hardware-implemented modules can provide information to, and receive information from, other hardware-implemented modules. Accordingly, the described hardware-implemented modules may be regarded as being communicatively coupled. Where multiple of such hardware-implemented modules exist contemporaneously, communications may be achieved through signal transmission (e.g., over appropriate circuits and buses) that connect the hardware-implemented modules. In embodiments in which multiple hardware-implemented modules are configured or instantiated at different times, communications between such hardware-implemented modules may be achieved, for example, through the storage and retrieval of information in memory structures to which the multiple hardware-implemented modules have access. For example, one hardware-implemented module may perform an operation and store the output of that operation in a memory device to which it is communicatively coupled. A further hardware-implemented module may then, at a later time, access the memory device to retrieve and process the stored output. Hardware-implemented modules may also initiate communications with input or output devices, and can operate on a resource (e.g., a collection of information).


Accordingly, the term “module” should be understood to encompass a tangible entity, be that an entity that is physically constructed, permanently configured (e.g., hardwired), non-transitory, or temporarily configured (e.g., programmed) to operate in a certain manner or to perform certain operations described herein. In some embodiments, modules or components may temporarily configured (e.g., programmed); and each of the modules or components need not be configured or instantiated at any one instance in time. For example, where the modules or components include a general-purpose processor configured using software, the general-purpose processor may be configured as respective different modules at different times. Software may accordingly configure the processor to constitute a particular module at one instance of time and to constitute a different module at a different instance of time.


Modules can provide information to, and receive information from, other modules. Accordingly, the described modules may be regarded as being communicatively coupled. Where multiples of such modules exist contemporaneously, communications may be achieved through signal transmission (e.g., over appropriate circuits and buses) that connect the modules. In embodiments in which multiple modules are configured or instantiated at different times, communications between such modules may be achieved, for example, through the storage and retrieval of information in memory structures to which the multiple modules have access. For example, one module may perform an operation and store the output of that operation in a memory device to which it is communicatively coupled. A further module may then, at a later time, access the memory device to retrieve and process the stored output. Modules may also initiate communications with input or output devices and can operate on a resource (e.g., a collection of information).


For purposes of the present description, the modules of FIG. 17 will be described in terms of the algorithms executed in each module, as may be executed by one or more processors, general purpose computer or other mechanism based on instructions stored in hardware in accordance with the description above.


In this example embodiment, the system 1700 includes a data access module 1704 configured to access MCI and other logging data from an initialization module 1702. An inversion module 1706 is configured to perform inversion in accordance with one or more of the example embodiments discussed with reference to FIG. 13, while a permeability determination module 1708 is configured to determine one or more permeability anisotropy components based on results of the inversion, according to the permeability anisotropy models and/or formulas discussed above. In some embodiments, the inversion module 1706 can be configured to perform the inversion processing as described herein. In one embodiment, the inversion module 1706 can be configured to calculate inverted TI parameters by performing an iterative inversion operation on the MCI data using a TI formation model that represents simulated formation resistivity characteristics that account for transverse formation isotropy to resistivity in the manner described herein. In other embodiments, the can be configured to calculate inverted BA parameters by performing an iterative inversion operation on the MCI data using a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity in the manner described herein.


After inversion processing, the permeability determination module 1708 can perform permeability anisotropy calculations as described herein. The system 1700 further comprises an output module 1710 configured to deliver the calculated permeability components. The output module 1710 may in some embodiments deliver numerical tables with estimated values for formation resistivity, laminar shale volume, and/or various permeability components and invasion resistivity at multiple different points along a borehole. In other embodiments, a graphical plot that maps the estimated values to the borehole positions may be printed in hard copy, and/or may be displayed on a display screen (e.g., video display unit 1810 as further described below in reference to FIG. 18).


Example Machine Architecture and Machine-Readable Medium


The various operations of example methods described herein may be performed, at least partially, by one or more processors that are temporarily configured (e.g., by software) or permanently configured to perform the relevant operations. Whether temporarily or permanently configured, such processors may constitute processor-implemented modules that operate to perform one or more operations or functions. The modules referred to herein may, in some example embodiments, comprise processor-implemented modules.


Similarly, the methods described herein may be at least partially processor-implemented. For example, at least some of the operations of a method may be performed by one or more processors or processor-implemented modules. The performance of certain operations may be distributed among the one or more processors, not only residing within a single machine, but deployed across a number of machines. In some example embodiments, the processor or processors may be located in a single location (e.g., within a home environment, an office environment or as a server farm), while in other embodiments the processors may be distributed across a number of locations.


The one or more processors may also operate to support performance of the relevant operations in a “cloud computing” environment or as a “software as a service” (SaaS). For example, at least some of the operations may be performed by a group of computers (as examples of machines including processors), with these operations being accessible via a network (e.g., the Internet) and via one or more appropriate interfaces (e.g., Application Program Interfaces (APIs).)



FIG. 18 is a diagrammatic representation of a machine in the example form of a computer system 1800 within which a set of instructions 1824 may be executed for causing the machine to perform any one or more of the methodologies discussed herein. For example, the surface computer system 366 (FIG. 3) or any one or more of its components may be provided by the system 1800.


In alternative embodiments, the machine operates as a standalone device or may be connected (e.g., networked) to other machines. In a networked deployment, the machine may operate in the capacity of a server or a client machine in a server-client network environment, or as a peer machine in a peer-to-peer (or distributed) network environment. The machine may be a server computer, a client computer, a personal computer (PC), a tablet PC, a set-top box (STB), a Personal Digital Assistant (PDA), a cellular telephone, a web appliance, a network router, switch or bridge, or any machine capable of executing a set of instructions (sequential or otherwise) that specify actions to be taken by that machine. Further, while only a single machine is illustrated, the term “machine” shall also be taken to include any collection of machines that individually or jointly execute a set (or multiple sets) of instructions to perform any one or more of the methodologies discussed herein.


The example computer system 1800 includes a processor 1802 (e.g., a central processing unit (CPU) a graphics processing unit (GPU) or both), a main memory 1804 and a static memory 1806, which communicate with each other via a bus 1808. The computer system 1800 may further include a video display unit 1810 (e.g., a liquid crystal display (LCD) or a cathode ray tube (CRT)). The computer system 1800 also includes an alpha-numeric input device 1812 (e.g., a keyboard), a cursor control device 1814 (e.g., a mouse), a disk drive unit 1816, a signal generation device 1818 (e.g., a microphone/speaker) and a network interface device 1820.


The disk drive unit 1816 includes a machine-readable or computer-readable storage medium 1822 on which is stored one or more sets of instructions 1824 (e.g., software) embodying any one or more of the methodologies or functions described herein. The instructions 1824 may also reside, completely or at least partially, within the main memory 1804 and/or within the processor 1802 during execution thereof by the computer system 1800, the main memory 1804 and the processor 1802 also constituting non-transitory machine-readable media. The instructions 1824 may further be transmitted or received over a network 1826 via the network interface device 1820.


While the machine-readable storage medium 1822 is shown in an example embodiment to be a single medium, the term “machine-readable medium” should be taken to include a single medium or multiple media (e.g., a centralized or distributed database and/or associated caches and servers) that store the one or more sets of instructions 1824. The term “machine-readable medium” shall also be taken to include any medium that is capable of storing a set of instructions for execution by the machine and that cause the machine to perform any one or more of the methodologies of this disclosure. The term “machine-readable medium” shall accordingly be taken to include, but not be limited to, solid-state memory devices of all types, as well as optical and magnetic media.


Although this disclosure has been described with reference to specific example embodiments, it will be evident that various modifications and changes may be made to these embodiments without departing from the broader scope of the disclosure. Accordingly, the specification and drawings are to be regarded in an illustrative rather than a restrictive sense.


The following numbered examples are illustrative embodiments in accordance with various aspects of the present disclosure.


1. A method of formation evaluation that includes accessing multicomponent induction (MCI) and permeability data captured by measurement tools in a borehole extending through a subsurface formation; determining inverted resistivity parameters by performing an iterative inversion operation on the MCI data; determining a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy; determining one or more permeability components based at least in part on the permeability data and permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio; and generating a visualization of the one or more permeability components for display on a display device.


2. The method of example 1, further including determining inverted transverse isotropy (TI) parameters by performing an iterative inversion operation on the MCI data.


3. The method of either of examples 1 or 2, in which the iterative inversion operation on the MCI data uses a TI formation model that represents simulated formation resistivity characteristics that account for transverse formation isotropy to resistivity.


4. The method of any of examples 1-3, further including determining inverted biaxial isotropy (BA) parameters by performing an iterative inversion operation on the MCI data.


5. The method of any of examples 1-4, in which the iterative inversion operation on the MCI data uses a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity.


6. The method of any of examples 1-5, in which determining one or more permeability components further includes determining biaxial permeability components for a TI anisotropic formation.


7. The method of any of examples 1-6, in which determining one or more permeability components further includes determining triaxial permeability components for a BA anisotropic formation.


8. The method of any of examples 1-7, in which determining one or more permeability components further includes determining permeability anisotropy based at least in part on the one or more permeability components.


9. The method of any of examples 1-8, further including determining horizontal and vertical permeabilities of a shale formation.


10. The method of any of examples 1-9, further including determining true sand permeability in a laminated formation.


11. The method of any of examples 1-10, in which the display device displays one or more permeability components at multiple points along the borehole.


12. A system includes a data access module to access multicomponent induction (MCI) and permeability data captured by measurement tools in a borehole extending through a subsurface formation; an inversion module to calculate inverted resistivity parameters by performing an iterative inversion operation based on the MCI data; a permeability determination module to calculate one or more permeability components based at least in part on the permeability data and the inverted resistivity parameters; and a display device for generating a visualization of the one or more permeability components.


13. The system of example 12, in which the inversion module is further configured to calculate inverted transverse isotropy (TI) parameters by performing an iterative inversion operation on the MCI data using a TI formation model that represents simulated formation resistivity characteristics that account for transverse formation isotropy to resistivity.


14. The system of any of the preceding examples, in which the inversion module is further configured to calculate inverted biaxial anisotropy (BA) parameters by performing an iterative inversion operation on the MCI data using a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity.


15. The system of any of the preceding examples, in which the permeability determination module is further configured to determine a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy.


16. The system of any of the preceding examples, in which the permeability determination module determines the one or more permeability components based at least in part on a permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio.


17. The system of any of the preceding examples, in which the permeability determination module is further configured to determine sand permeability in a laminated formation.


18. A computer readable storage medium having stored thereon instructions for causing a machine, in response to execution of the instructions by the machine, to perform operations that include accessing multicomponent induction (MCI) and permeability data captured by measurement tools in a borehole extending through a subsurface formation; determining inverted resistivity parameters by performing an iterative inversion operation on the MCI data; determining a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy; determining one or more permeability components based at least in part on the permeability data and permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio; and generating a visualization of the one or more permeability components for display on a display device.


19. The computer readable storage medium of example 18, in which the instructions are executable to further perform operations including determining horizontal and vertical permeabilities of a shale formation using a multi-component permeability model including isotropic sands and transversely isotropic shale.


20. The computer readable storage medium of example 18, in which the instructions are executable to further perform operations including determining sand permeability in a laminated formation using a bimodal permeability model.


In this description, references to “one embodiment” or “an embodiment,” or to “one example” or “an example,” are not intended necessarily to refer to the same embodiment or example; however, neither are such embodiments mutually exclusive, unless so stated or as will be readily apparent to those of ordinary skill in the art having the benefit of this disclosure. Thus, a variety of combinations and/or integrations of the embodiments and examples described herein may be included, as well as further embodiments and examples as defined within the scope of all claims based on this disclosure, and all legal equivalents of such claims.


In the foregoing Detailed Description, it can be seen that various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the following claims reflect, disclosed subject matter lies in less than all features of a single disclosed embodiment. Thus, the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment.

Claims
  • 1. A method of formation evaluation, comprising: accessing multicomponent induction (MCI) and isotropic permeability data captured by measurement tools in a borehole extending through a subsurface formation;determining inverted resistivity parameters by performing an iterative inversion operation on the MCI data, wherein determining the inverted resistivity parameters comprises calculating inverted biaxial anisotropy (BA) parameters by performing an iterative inversion operation on the MCI data using a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity;determining a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy;determining one or more permeability components based at least in part on the isotropic permeability data and permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio;generating a visualization of the one or more permeability components for display on a display device; andoperating a controlled device based at least a on the one or more permeability components.
  • 2. The method of claim 1, further comprising: determining inverted transverse isotropy (TI) parameters by performing an iterative inversion operation on the MCI data.
  • 3. The method of claim 2, wherein the iterative inversion operation on the MCI data uses a TI formation model that represents simulated formation resistivity characteristics that account for transverse formation isotropy to resistivity.
  • 4. The method of claim 1, wherein the iterative inversion operation on the MCI data uses a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity.
  • 5. The method of claim 1, wherein determining one or more permeability components further comprises determining biaxial permeability components for a TI anisotropic formation.
  • 6. The method of claim 1, wherein determining one or more permeability components further comprises determining triaxial permeability components for a BA anisotropic formation.
  • 7. The method of claim 1, wherein determining one or more permeability components further comprises determining permeability anisotropy based at least in part on the one or more permeability components.
  • 8. The method of claim 1, further comprising: determining horizontal and vertical permeabilities of a shale formation.
  • 9. The method of claim 1, further comprising: determining true sand permeability in a laminated formation.
  • 10. The method of claim 1, wherein the display device displays one or more permeability components at multiple points along the borehole.
  • 11. A system comprising: a data access module to access multicomponent induction (MCI) and isotropic permeability data captured by measurement tools in a borehole extending through a subsurface formation;an inversion module to calculate inverted resistivity parameters by performing an iterative inversion operation based on the MCI data;a permeability determination module to calculate one or more permeability components based at least in part on the isotropic permeability data and the inverted resistivity parameters; anda display device for generating a visualization of the one or more permeability components,wherein the inversion module is further configured to calculate inverted biaxial anisotropy (BA) parameters by performing an iterative inversion operation on the MCI data using a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity; andwherein a controlled device is operated based at least in part on the one or more permeability components.
  • 12. The system of claim 11, wherein the inversion module is further configured to calculate inverted transverse isotropy (TI) parameters by performing an iterative inversion operation on the MCI data using a TI formation model that represents simulated formation resistivity characteristics that account for transverse formation isotropy to resistivity.
  • 13. The system of claim 11, wherein the permeability determination module is further configured to determine a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy.
  • 14. The system of claim 13, wherein the permeability determination module determines the one or more permeability components based at least in part on a permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio.
  • 15. The system of claim 11, wherein the permeability determination module is further configured to determine sand permeability in a laminated formation.
  • 16. A non-transitory computer readable storage medium having stored thereon instructions for causing a machine, in response to execution of the instructions by the machine, to perform operations comprising: accessing multicomponent induction (MCI) and isotropic permeability data captured by measurement tools in a borehole extending through a subsurface formation;determining inverted resistivity parameters by performing an iterative inversion operation on the MCI data, wherein determining the inverted resistivity parameters comprises calculating inverted biaxial anisotropy (BA) parameters by performing an iterative inversion operation on the MO data using a BA formation model that represents simulated formation resistivity characteristics that account for transverse biaxial formation anisotropy to resistivity;determining a resistivity anisotropy ratio based at least in part on the inverted resistivity parameters and a permeability anisotropy ratio based at least in part on the formation resistivity anisotropy;determining one or more permeability components based at least in part on the isotropic permeability data and permeability-resistivity relationship between the resistivity anisotropy ratio and the permeability anisotropy ratio;generating a visualization of the one or more permeability components for display on a display device; andoperating a controlled device based at east in part on the one or more permeability components.
  • 17. The computer readable storage medium of claim 16, wherein the instructions are executable to further perform operations comprising: determining horizontal and vertical permeabilities of a shale formation using a multi-component permeability model comprising isotropic sands and transversely isotropic shale.
  • 18. The computer readable storage medium of claim 16, wherein the instructions are executable to further perform operations comprising: determining sand permeability in a laminated formation using a bimodal permeability model.
PCT Information
Filing Document Filing Date Country Kind
PCT/US2016/014679 1/25/2016 WO 00
Publishing Document Publishing Date Country Kind
WO2017/131608 8/3/2017 WO A
US Referenced Citations (11)
Number Name Date Kind
5663499 Semmelbeck Sep 1997 A
7157915 Schoen Jan 2007 B2
20020101235 Schoen Aug 2002 A1
20040190589 Zazovsky Sep 2004 A1
20060241867 Kuchuk Oct 2006 A1
20090145600 Wu et al. Jun 2009 A1
20100004866 Rabinovich Jan 2010 A1
20100185424 Massonnat Jul 2010 A1
20100305927 Suarez-Rivera Dec 2010 A1
20140257706 Biterge Sep 2014 A1
20170363765 Hou Dec 2017 A1
Foreign Referenced Citations (2)
Number Date Country
WO 2015051300 Apr 2015 WO
WO 2015161282 Oct 2015 WO
Non-Patent Literature Citations (10)
Entry
Adams (MIT), 2013, Permeability anisotropy and resistivity anisotropy of mudrocks.
Aldred, R.D., 2015, 3D petrophysical modeling to resolve complex reservoir interpretation problems: A case study from South East Asia: SPE 176243.
Georgi, D. et al., 2002, On the relationship between resistivity and permeability anisotropy: SPE 77715.
Hagiwara, T., 2015, To estimate permeability anisotropy from resistivity anisotropy: Resistivity SIG meeting in Apr. 22, 2015, Houston, USA.
International Search Report and The Written Opinion of the International Search Authority, or the Declaration, dated Sep. 12, 2016, PCT/US2016/014679, 11 pages, ISA/KR.
Klein et al., 1995, The petrophysics of electrically anisotropic reservoirs: SPWLA 36th Annual Logging Symposium, Paper HH.
Moran et al., EE, 1962, Theoretical analysis of pressure phenomena associated with wireline formation tester: Journal of petroleum technology, v14, 899-908.
Srisutthiyakorn et al., 2015, An improved Kozeny-Carman for irregular pore geometries: SEG extended abstract.
Tittman, 1990, Formation anisotropy: Reckoning with its effects. Schlumberger Oilfield Review, vol. 2(1), 16-23.
Wyllie et al., MB, 1952, Application of electrical resistivity measurements to problem of fluid flow in porous media: Bulletin of the American association of petroleum geologists: v36, 359-403.
Related Publications (1)
Number Date Country
20200309985 A1 Oct 2020 US