A method may comprise: inserting into a wellbore penetrating a subterranean formation an apparatus comprising: a transmitter sub comprising a transmitter coil; a first receiver sub comprising a first receiver coil; a second receiver sub comprising a second receiver coil; and a third receiver sub comprising a third receiver coil; generating an electromagnetic wave at the transmitter coil; propagating the electromagnetic wave through the subterranean formation; receiving the electromagnetic wave in the first receiver coil, the second receiver coil, and the third receiver coil; generating a plurality of response signals, wherein the plurality of response signals comprises a first response signal in the first receiver coil, a second response signal in the second receiver coil, and a third response signal the third receiver coil; solving a series of linear equations defined by an orientation of the transmitter sub in the wellbore, orientation in the wellbore for each of the first receiver coil, the second receiver coil, and the third receiver coil, and the plurality of response signals to determine a Green's response tensor.

Patent
   10928542
Priority
Jun 07 2018
Filed
Jun 07 2018
Issued
Feb 23 2021
Expiry
Jul 03 2038
Extension
26 days
Assg.orig
Entity
Large
0
18
currently ok
1. A method comprising:
inserting into a wellbore penetrating a subterranean formation an apparatus comprising:
a transmitter sub comprising a transmitter coil;
a first receiver sub comprising a first receiver coil;
a second receiver sub comprising a second receiver coil; and
a third receiver sub comprising a third receiver coil;
generating an electromagnetic wave at the transmitter coil;
propagating the electromagnetic wave through the subterranean formation;
receiving the electromagnetic wave in the first receiver coil, the second receiver coil, and the third receiver coil;
generating a plurality of response signals, wherein the plurality of response signals comprises a first response signal in the first receiver coil, a second response signal in the second receiver coil, and a third response signal the third receiver coil;
solving a series of linear equations defined by an orientation of the transmitter sub in the wellbore, orientation in the wellbore for each of the first receiver coil, the second receiver coil, and the third receiver coil, and the plurality of response signals to determine a Green's response tensor
wherein the series of linear equations have the following form for each of the response signals:
V ( θ T , β T , θ R , β R ) = [ sin θ T cos β T sin θ T sin β T cos θ T ] [ z xx z x y z xz z y x z y y z y z z zx z z y z zz ] [ sin θ R cos β R sin θ R sin β R cos θ R ]
where V(θT, βT, θR, βR) is the response signal for receiver coil i, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, βR is a tilt angle of the receiver coil i, βR is an azimuth angle of the receiver coil i, zxx is a component of the response signal for receiver coil i oriented in an x direction caused by the transmitter coil in the x direction, zxy is a component of the response signal for receiver coil i oriented in a y direction caused by the transmitter coil in the x direction, zzz is a component of the response signal for receiver coil i oriented in an z direction caused by the transmitter coil in the x direction, zxy is a component of the response signal for receiver coil i oriented in the x direction caused by the transmitter coil in the y direction, zyy is a component of the response signal for receiver coil i oriented in the y direction caused by the transmitter coil in the y direction, zyz is a component of the response signal for receiver coil i oriented in the z direction caused by the transmitter coil in the y direction, zzx is a component of the response signal for receiver coil i oriented in the x direction caused by the transmitter coil in the z direction, zzy i is a component of the response signal for receiver coil i oriented in y direction caused by the transmitter coil in the z direction, zzz is a component of the response signal for receiver coil i oriented in the z direction caused by the transmitter coil in the z direction, θR is a tilt angle of the receiver coil i, and βR is an azimuth angle of the receiver coil i, wherein receiver coil i is one of the first receiver coil, the second receiver coil, or the third receiver coil.
2. A method comprising:
inserting into a wellbore penetrating a subterranean formation an apparatus comprising:
a transmitter sub comprising a transmitter coil;
a first receiver sub comprising a first receiver coil;
a second receiver sub comprising a second receiver coil; and
a third receiver sub comprising a third receiver coil;
generating an electromagnetic wave at the transmitter coil;
propagating the electromagnetic wave through the subterranean formation;
receiving the electromagnetic wave in the first receiver coil, the second receiver coil, and the third receiver coil;
generating a plurality of response signals, wherein the plurality of response signals comprises a first response signal in the first receiver coil, a second response signal in the second receiver coil, and a third response signal the third receiver coil;
solving a series of linear equations defined by an orientation of the transmitter sub in the wellbore, orientation in the wellbore for each of the first receiver coil, the second receiver coil, and the third receiver coil, and the plurality of response signals to determine a Green's response tensor
wherein the series of linear equations has the following form for the response signal corresponding to each receiver coil;
for the first receiver coil:
V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where v1T, βT, θR1, βR1) is the response signal of the first receiver coil, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, zxxR1 is a component of the first response signal oriented in an x direction caused by the transmitter coil in the x direction, zxyR1 is a component of the first response signal oriented in a y direction caused by the transmitter coil in the x direction, zXZR1 is a component of the first response signal oriented in a z direction caused by the transmitter coil in the x direction, zyxR1 is a component of the first response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR1 is a component of the first response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR1 is a component of the first response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR1 is a component of the first response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR1 is a component of the first response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR1 is a component of the first response signal oriented in the z direction caused by the transmitter coil in the z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the first receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where v2 T, β7, θR2, βR2) is the response signal of the second receiver coil, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR2 is a component of the second response signal oriented in the x direction, zyxR2 is a component of the second response signal oriented in the y direction, zyyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the z direction, θR2 is a tilt angle of the second receiver coil, and βR2 is an azimuth angle of the second receiver coil;
for the third receiver coil:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where V3TTR3R3) is the third response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the x direction, zyxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR3 i is a component of the third response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the z direction, θR3 is a tilt angle of the receiver coil, and βR3 is an azimuth angle of the third receiver coil.
10. An apparatus comprising:
a wellbore drill bit
a transmitter sub comprising a transmitter coil;
a first receiver sub comprising a first receiver coil;
a second receiver sub comprising a second receiver coil; and
a third receiver sub comprising a third receiver coil;
wherein the transmitter sub is configured to generate an electromagnetic wave at the transmitter coil and propagate the electromagnetic wave through a subterranean formation, wherein the first receiver coil is configured to receive the electromagnetic wave and generate a first response signal, wherein the second receiver coil is configured to receive the electromagnetic wave and generate a second response signal, and wherein the third receiver coil is configured to receive the electromagnetic wave and generate a third response signal;
at least one processor; and
a memory coupled to the processor to provide software that configures the processor to solve a series of linear equations defined by an orientation of the transmitter sub in the wellbore, an orientation of the first receiver coil, an orientation of the second receiver coil, and an orientation the third receiver coil, and each of the first response signal, the second response signal, and the third response signal to determine a Green's response tensor
wherein the series of linear equations has the following form for a response signal corresponding to each receiver, the response signal being the first response signal, the second response signal, or the third response signal;
for the first receiver:
V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where v1T, βT, θR1, βR1) is the first response signal, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, zxxR1 a component of the first response signal oriented in an x firection caused by the transmitter coil in the x direction, zxyR1 a component of the first response signal oriented in a y direction caused by the transmitter coil in the x direction, zxzR1 a component of the first response signal oriented in a z direction caused by the transmitter coil in the x direction, zyxR1 a component of the first response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR1 a component of the first response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR1 a component of the first response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR1 a component of the first response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR1 a component of the first response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR1 a component of the first response signal oriented in the z direction caused by the transmitter coil in the z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where V2T, βT, θR2, β2) is the second response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the x direction, zzxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the y direction, zzzR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the z direction, θR2 is a tilt angle of the second receiver coil, and ΔR2 is an azimuth angle of the second receiver coil;
for the third receiver:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where v3 T, βT, θR3, βRO is the third response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxyR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the x direction, zyxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the y direction, zxxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the z direction, zxxR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the z direction, θR3 is a tilt angle of the receiver coil, and βR3 is an azimuth angle of the third receiver coil.
14. A system comprising:
a tool comprising:
a transmitter sub comprising a transmitter coil;
a first receiver sub comprising a first receiver coil;
a second receiver sub comprising a second receiver coil; and
a third receiver sub comprising a third receiver coil;
wherein the transmitter sub is configured to generate an electromagnetic wave at the transmitter coil and propagate the electromagnetic wave through a subterranean formation, wherein the first receiver coil is configured to receive the electromagnetic wave and generate a first response signal, wherein the second receiver coil is configured to receive the electromagnetic wave and generate a second response signal, and wherein the third receiver coil is configured to receive the electromagnetic wave and generate a third response signal;
a conveyance electrically coupled to the tool; and
a surface processing unit comprising:
at least one processor; and
a memory coupled to the processor to provide software that configures the processor to solve a series of linear equations defined by an orientation of the transmitter sub, first receiver sub, second receiver sub, and third receiver sub orientation in the wellbore, and the first response signal, the second response signal, and the first response signal to determine a Green's response tensor wherein the series of linear equations has the following form for a response signal corresponding to each receiver, the response signal being the first response signal, the second response signal, or the third response signal;
for the first receiver coil:
V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where V1t, Δt, θR1, βR1) is the first response signal, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, zxxR1 is a component of the first response signal oriented in an x direction caused by the transmitter coil in the x direction, zxyR1 is a component of the first response signal oriented in a y direction caused by the transmitter coil in the x direction, zxzR1 is a component of the first response signal oriented in a z direction caused by the transmitter coil in the x direction, zyxR1 is a component of the first response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR1 is a component of the first response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR1 is a component of the first response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR1 is a component of the first response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR1 is a component of the first response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR1 is a component of the first response signal oriented in the z direction caused by the transmitter coil in the z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where V2T, βT, θR2, βR2) is the second response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the x direction, zyxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR2 is a component of the second response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR2 is a component of the second response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR2 is a component of the second response signal oriented in the z direction caused by the transmitter coil in the z direction, θR2 is a tilt angle of the second receiver coil, and βR2 is an azimuth angle of the second receiver coil;
for the third receiver:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where v3T, βT, θR3, βR) is the third response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, zxxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the x direction, zxyR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the x direction, zxzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the x direction, zyxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the y direction, zyyR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the y direction, zyzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the y direction, zzxR3 is a component of the third response signal oriented in the x direction caused by the transmitter coil in the z direction, zzxR3 is a component of the third response signal oriented in the y direction caused by the transmitter coil in the z direction, zzzR3 is a component of the third response signal oriented in the z direction caused by the transmitter coil in the z direction, θR3 is a tilt angle of the receiver coil, and βR3 is an azimuth angle of the third receiver coil.
3. The method of claim 2 wherein each azimuth angle is found by the following equations:
azimuth angle of the transmitter coil:
β T = n m · 360
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:

βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver coil:

βR2Toff
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:

βR2ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil.
4. The method of claim 3 wherein an assumption of zR1≈ZR2≈ZR3 is applied to each response signal equation such that:
V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · z yy R 22 + sin θ T sin β T cos ( n m · 360 ° + β 0 ff ) · z yz R 2 + cos θ T sin θ R 2 cos β R 2 · z zx R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.
5. The method of claim 4 wherein each of the response signals with respect to n is:
[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ z xx R 2 z xy R 2 z x z R 2 z y x R 2 z y y R 2 z y z R 2 z z x R 2 z z y R 2 z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each of the response signals with respect to n for the first response signal, the second response signal, and the third response signal from n to m.
6. The method of claim 5 wherein the step solving the series of linear equations comprises solving for a vector defined as:
[ z xx R 2 z xy R 2 z x z R 2 z y x R 2 z y y R 2 z y z R 2 z z x R 2 z z y R 2 z zz R 2 ] .
7. The method of claim 6 wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.
8. The method of claim 7 wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor.
9. The method of claim 8 wherein λ is determined by an L-curve method.
11. The apparatus of claim 10 wherein each azimuth angle is found by the following equations:
azimuth angle of the transmitter coil:
β T = n m · 360 °
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:

βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver:

βR3ToffΔ2
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:

βR3ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil,
and wherein an assumption of zR1≈ZR2≈ZR3 is applied to each response signal equation such that:
V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ1 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · z yy R 2 2 + sin θ T sin β T cos ( n m · 360 ° + β off ) · z yz R 2 + cos θ T sin θ R 2 cos β R 2 · z z x R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.
12. The apparatus of claim 11 wherein each response signal with respect to n is:
[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each response signal with respect to n for the first receiver, the second receiver, and the third receiver from n to m, wherein the step solving the series of linear equations comprises solving for a vector defined as:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] .
and wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.
13. The apparatus of claim 11 wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor and wherein λ is determined by an L-curve method.
15. The apparatus of claim 14 wherein each azimuth angle is found by the following equations:
azimuth angle of the transmitter coil:
β T = n m · 360 °
where βT is the azimuth angle of the transmitter coil, m is a number of times the receivers sample per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the transmitter coil:
β T = n m · 360 °
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:

βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver:

βR2Toff
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:

βR2ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil,
and wherein an assumption of zR1≈ZR2≈ZR3 is applied to each response signal equation such that:
V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ1 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · z yy R 2 2 + sin θ T sin β T cos ( n m · 360 ° + β off ) · z yz R 2 + cos θ T sin θ R 2 cos β R 2 · z z x R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.
16. The apparatus of claim 15 wherein each response signal with respect to n is:
[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each of the first response signal, second response signal, and third response signal with respect to index n to m, wherein the step solving the series of linear equations comprises solving for a vector defined as:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ]
and wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.
17. The apparatus of claim 15 wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:
[ z xx R 2 z xy R 2 z xz R 2 z yx R 2 z yy R 2 z yz R 2 z zx R 2 z zy R 2 z zz R 2 ] = ( B _ _ m × 9 H = B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor and wherein λ is determined by an L-curve method.

Modern petroleum drilling and production operations may demand a great quantity of information relating to the parameters and conditions downhole. Such information typically includes the location and orientation of the borehole and drilling assembly, earth formation properties, and parameters of the downhole drilling environment. The collection of information relating to formation properties and downhole conditions is commonly referred to as “logging”, and can be performed during the drilling process itself (hence the term “logging while drilling” or “LWD,” frequently used interchangeably with the term “measurement while drilling” or “MWD”).

Various measurement tools exist for use in LWD. One such tool is the resistivity tool, which includes one or more antennas for transmitting an electromagnetic signal into the formation and one or more antennas for receiving a formation response. Physical phenomena that determine how the measurement is made may vary with frequency, for example. In some examples, the amplitude and/or the phase of the receive signals are compared to the amplitude and/or phase of the transmit signals to measure the formation resistivity. In other cases, the amplitude and/or phase of multiple receive signals are compared to each other to measure the formation resistivity.

When plotted as a function of depth or tool position in the borehole, the logging tool measurements are termed “logs.” Resistivity logging may be used in well logging to determine geological correlation of formation strata and detect and quantify potentially productive formation zones. Such logs may provide indications of hydrocarbon concentrations and other information useful to drillers and completion engineers. In particular, azimuthally-sensitive logs may provide information useful for steering the drilling assembly because they can inform the driller when a target formation bed has been entered or exited, thereby enabling modifications to the drilling program that will provide much more value and higher success than would be the case using only seismic data.

Generally, in real formations, a resistivity of the formation varies in different directions, for example, a formation resistivity may vary in the x, y, and z coordinates. Previous methods of determining a formation resistivity may have assumed a uniaxial vertical 1-dimensional (V-1D) structure where resistivity is assumed to be constant on an x-y plane, for example, with resistivity only varying in the z direction. In some uncomplicated formations, an assumption of V1-D may be appropriate as the assumption of V1-D may closely match the true formation resistivity. However, in complicated formations with real resistivity varying in three directions, the assumption of V-1D may not produce an accurate enough picture of formation resistivity. Using an assumption of V-1D may result in an incomplete picture of formation resistivity and thereby a sub-optimal drilling path may be followed. As a consequence, a full Green's Tensor for resistivity in complicated formations may not be able to be calculated using the assumption of V-1D.

A better understanding of the various disclosed embodiments may be obtained when the following detailed description is considered in conjunction with the attached drawings, in which:

FIG. 1 shows an illustrative logging while drilling (LWD) environment.

FIG. 2 shows an illustrative drillstring with two logging tool modules.

FIG. 3 shows an illustrative logging tool and a surface system.

FIG. 4 shows an illustrative logging tool.

It should be understood that the drawings and corresponding detailed description do not limit the disclosure, but on the contrary, they provide the foundation for understanding all is modifications, equivalents, and alternatives falling within the scope of the appended claims.

This disclosure may generally relate to apparatuses, systems and methods for producing deep formation evaluations using logging while drilling (“LWD”) tools with multiple subs. In other examples, this disclosure may relate to apparatuses, systems and methods for producing deep formation evaluations using wireline logging having multiple subs.

The resistivity of a given formation may be isotropic (equal in all directions) or anisotropic (unequal in different directions). In electrically anisotropic formations, the anisotropy is generally attributable to extremely fine layering during the sedimentary buildup of the formation. As a result, in a formation Cartesian coordinate system oriented such that the x-y plane is parallel to the formation layers and the z axis is perpendicular to the formation layers, resistivities Rx and Ry in the x and y directions, respectively, may tend to be similar, but resistivity Rz in the z direction may tend to be different. The resistivity in a direction parallel to the formation plane (i.e., the x-y plane) may be referred to as the horizontal resistivity, Rh, and the resistivity in the direction perpendicular to the plane of the formation (i.e., the z direction) may be referred to the vertical resistivity, Rv. The index of anisotropy, η, is defined as r|=[Rv/Rh]1/2.

As previously discussed, a complex formation may have resistivity that is directionally dependent in three dimension. Such a formation may be referred to as anisotropic. A resistivity of a formation may also be described as isotropic or uniaxial. For a formation that is isotropic, Equation 1 may describe the resistivity. An isotropic formation may have no variation in resistivity with orientation in the formation. Equation 2 may describe a uniaxial formation where the resistivity does not vary in the x-y plane but varies with z. Although illustrated in the x-y plane in Equation 2, the plane may also be x-z or y-z. As previously discussed, an assumption that a formation behaves in a uniaxial manner may be referred to as V1-D. Equation 3 illustrates the case of an anisotropic formation resistivity where the resistivity is directionally dependent in three dimensions.
F(x,y,z)=[f(x,y,z)]  (1)
F(x,y,z)=[f_xy(x,y,z),f_z(x,y,z,)]  (2)
F(x,y,z)=[f_x(x,y,z),f_y(x,y,z),f_z(x,y,z)]  (3)

The disclosed apparatuses, systems and methods may be best understood in the context of the larger systems in which they operate. FIG. 1 illustrates a diagrammatic view of an exemplary logging while drilling (LWD) and/or measurement while drilling (MWD) wellbore operating environment 100 in which the present disclosure may be implemented. As depicted in FIG. 1, 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, the drill bit 114 creates a wellbore 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 the annulus around drill string 108, and into a retention pit 124. The drilling fluid transports cuttings from the wellbore 116 into the pit 124 and aids in maintaining the integrity of the wellbore 116. Various materials can be used for drilling fluid, including oil-based fluids and water-based fluids.

As depicted in FIG. 1, logging tools 126 are integrated into the bottom-hole assembly 125 near the drill bit 114. As the drill bit 114 extends the wellbore 116 through the formations 118, logging tools 126 collect measurements relating to various formation properties as well as the orientation of the tool and various other drilling conditions. The bottom-hole assembly 125 may also include a telemetry sub 128 to transfer measurement data to a surface receiver 130 and to receive commands from the surface. In some embodiments, the telemetry sub 128 communicates with a surface receiver 130 using mud pulse telemetry. In other cases, 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. Notably, one or more of the bottom-hole assembly 125, the logging tools 126, and the telemetry sub 128 may also operate using a non-conductive cable (e.g. slickline, etc.) with a local power supply, such as batteries and the like. When employing non-conductive cable, communication may be supported using, for example, wireless protocols (e.g. EM, acoustic, etc.) and/or measurements and logging data may be stored in local memory for subsequent retrieval at the surface, as is appreciated by those skilled in the art.

Each of the logging tools 126 may include a plurality of tool components, spaced apart from each other, and communicatively coupled with one or more wires. Logging tools 126 may include tools such as the one shown in FIG. 4 in order to perform resistivity, or conductivity logging. The telemetry sub 128 may include wireless telemetry or logging capabilities, or both, such as to transmit or later provide information indicative of received energy/waveforms to operators on the surface or for later access and data processing for the evaluation of formation 118 properties.

The logging tools 126, including the resistivity logging tool, may also include one or more computing device 150 communicatively coupled with one or more of the plurality of tool components. Computing device 150 may comprise components such as processors, memory, input/output devices, interfaces, and the like. The computing device 150 may be configured to control or monitor the performance of the tools 126, process logging data, and/or carry out the methods of the present disclosure.

In some embodiments, one or more of the logging tools 126 may communicate with a surface receiver 130, such as a wired drillpipe. In other cases, the one or more of the logging tools 126 may communicate with a surface receiver 130 by wireless signal transmission. In at least some cases, one or more of the logging tools 126 may receive electrical power from a wire that extends to the surface, including wires extending through a wired drillpipe. In at least some instances the methods and techniques of the present disclosure may be performed by a computing device (not shown) on the surface. In some embodiments, the computing device may be included in surface receiver 130. For example, surface receiver 130 of wellbore operating environment 100 at the surface may include one or more of wireless telemetry, processor circuitry, or memory facilities, such as to support substantially real-time processing of data received from one or more of the logging tools 126. In some embodiments, data is processed at some time subsequent to its collection, wherein the data may be stored on the surface at surface receiver 130, stored downhole in telemetry sub 128, or both, until it is retrieved for processing.

FIG. 2 shows an illustrative example of a deep formation evaluation logging tool that includes two LWD tool modules 202 and 206 at different locations and orientations along a drillstring. In the example shown, a resistivity logging tool receive antenna 212 and a corresponding receive antenna position measurement device 222a may be housed within LWD tool module 202, while a resistivity logging tool transmit antenna 216 and a corresponding transmit antenna position measurement device 222b (components of an “at bit” instrument) are housed within LWD tool module 206. The position measurement devices may locate the position of each corresponding antenna, which may be expressed, for example, in terms of each antenna's tilt angle (θr and θT relative to the zr and zt axes respectively; generally fixed and known), each antenna's azimuthal angle (αr and αt relative to the x axis), each LWD tool module's inclination angle (φr and φt) and the distance d′ between the antennas.

As a further complication to measuring formation resistivity, boreholes are generally perpendicular to formation beds. The angle between the axis of the well bore and the orientation of the formation beds (as represented by the normal vector) has two components. These components are the dip angle and the azimuth angle. The dip angle is the angle between the borehole axis and the normal vector for the formation bed. The azimuth angle is the direction in which the borehole's axis “leans away from” the normal vector. Electromagnetic resistivity logging measurements are a complex function of formation resistivity, formation anisotropy, and the formation dip and azimuth angles, which may all be unknown. A triaxial induction well logging tool may be used to detect formation properties such as resistivity anisotropy, which is one of the important parameters in evaluation subterranean formations such as sand-shale reservoirs or fractured reservoirs. However, the resistivity anisotropy parameter cannot be obtained without performing a numerical inversion process. Specifically, numerical inversion may be required to obtain accurate formation resistivity anisotropy parameters. The log inversion required for anisotropy determination may involve a large number of inversion parameters to be determined by an algorithm referred to as the ID vertical inversion. Generally, this algorithm may require large amounts of processing time and be sensitive to noise from logging, the logging environment characteristics and borehole correction, which could result in errors in the inverted vertical resistivity.

FIG. 3 illustrates a diagrammatic view of a conveyance logging wellbore operating environment 300 in which the present disclosure may be implemented. As depicted in FIG. 3, a hoist 306 may be included as a portion of a platform 302, such as that coupled to derrick 304, and used with a conveyance 342 to raise or lower equipment such as resistivity logging tool 310 into or out of a borehole. Resistivity logging tool 310 may include, for example, tools such as the one shown in FIG. 4. A conveyance 342 may provide a communicative coupling between the resistivity logging tool 310 and a logging facility 344 at the surface. The conveyance 342 may include wires (one or more wires), slicklines, cables, or the like, as well as tubular conveyances such as coiled tubing, joint tubing, or other tubulars, and may include a downhole tractor. Additionally, power can be supplied via the conveyance 342 to meet power requirements of the tool. The resistivity logging tool 310 may have a local power supply, such as batteries, downhole generator and the like. When employing non-conductive cable, coiled tubing, pipe string, or downhole tractor, communication may be supported using, for example, wireless protocols (e.g. EM, acoustic, etc.), and/or measurements and logging data may be stored in local memory for subsequent retrieval. The logging facility 344 may include a computing device 350 capable of carrying out the methods and techniques of the present disclosure. In this manner, information about the formation 318 may be obtained by resistivity logging tool 310 and processed by a computing device, such as computing device 350. In some embodiments, computing device 350 is equipped to process the received information in substantially real-time, while in some embodiments, computing device 350 can be equipped to store the received information for processing at some subsequent time.

FIG. 4 illustrates an example wellbore tool 400 that may be used in the systems and methods described herein. Wellbore tool 400 may comprise transmitter sub 405 and one or more receiver subs 410, 411, and 412. In some examples, transmitter sub 405 may be referred to as TX and receiver subs 410, 411, and 412 may be referred to as RX1, RX2, and RX3 respectively. Transmitter sub 405 may comprise a transmitter coil 406 which may be an electromagnetic wave source such as a monopole, dipole, quadrupole, or other higher order wave source. Each of the receiver subs 410, 411, and 412 may comprise three or more receiver coils per sub configured to receive an electromagnetic wave from transmitter sub 405. Receiver subs 410, 411, and 412 may be disposed on wellbore tool 400 a distance 115, 116, 117 from transmitter sub 105. Distance 115, 116, 117 may also be referred to as S1, S2, and S3 respectively.

As illustrated in FIG. 4, first receiver coils 420 are not co-axial with receiver sub 410. There may be an axial offset between first receiver coils 420 and a centerline of receiver sub 410 which may be notated as θR1. Similarly, for transmitter coil 406, second receiver coils 421, and third receiver coils 422, there may be an axial offset of coils from a centerline of the respective subs notated as θT, θR2, and θR3 respectively. In addition to axial offset, each of the first receiver coils 420, second receiver coils 421, and third receiver coils 422 may have an azimuthal offset relative to transmitter coil 406.

Table 1 illustrates the notations for wellbore tool 400 of FIG. 1. The tilt angle of the transmitter coil is notated as θT and the tilt angle, or azimuthal offset, of each of the receiver coils is notated as θR1, θR2, and θR3 for RX1, RX2, and RX3 respectively. The azimuth angle is dependent on the tool's rotated position in the wellbore. βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil which can be measured before the tool is inserted into the wellbore. βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil and βΔ2 the difference in azimuthal angle between the third receiver coil and second receiver coil. βΔ1 and βΔ2 can take any value, but must meeting the following requirements: (1) βΔ1≠βΔ2, (2) βΔ1≠0°, and (3) βΔ2≠0°.

As will be further illustrated below, a system of linear equations may be generated from the wellbore tool's measurements in the wellbore. In Table 1, m represents the total number of measurements per turn of the wellbore tool and n represents the current measurement number. For example, a wellbore tool may take m=32 measurements per turn and n would count 1 through 32. If m is 32, a system of 32 linear equations may be generated.

TABLE 1
Tilt Distance Measured
Antenna # Angle Azimuth Angle From TX Response
TX θT β T = n m · 360 ° ( n = 1 , 2 , , m ) 0 N.A.
RX1 θR1 βR1 = βT + βoff + βΔ1 S1 V1T, βT, θR1, βR1)
RX2 θR2 βR2 = βT + βoff S2 V2T, βT, θR2, βR2)
RX3 θR3 βR3 = βT + βoff + βΔ2 S3 V3T, βT, θR3, βR3)

FIG. 3 illustrates a multicomponent tool which can provide measurements V1, V2, and V3. The measurements may be used to derive Green's response tensor illustrated in Equation 4. The elements of the Green's response tensor represent the impedance response signal received from each RX on the tool. Each element has two subscripts. The first subscript denotes the component of TX and the second subscript denotes the component of RX. For example, Zxy means the response signal of an RX (receiver coil) oriented in Y coordinate direction caused by signal from TX oriented in X direction. The X, Y, and Z directions are defined by the tool coordinate system wherein the coordinate system origin is at the center of TX and the Z direction aligns with the direction from the center of RX2 (the middle RX) to the center of TX. The X direction is perpendicular to Z direction and parallel to layer boundary and the Y direction is perpendicular to both X and Z direction.

[ Z x x Z x y Z x z Z y x Z y y Z y z Z z x Z z y Z zz ] ( 4 )

As previously discussed, an assumption of uniaxial vertical-1D formation structure has been used to simplify resistivity measurements. Equation 5 illustrates that a V1-D assumption will reduce the Greens response tensor to 5 non-zero components. In real formations with anisotropic resistivity each component of the Green's response tensor is needed to accurately report resistivity in each direction.

[ Z x x 0 Z x z 0 Z y y 0 Z z x 0 Z zz ] ( 5 )

A general formula for a response signal of an RX is illustrated in Equation 6 and the expanded formula is illustrated in equation 7.

V = ( θ T , β T , θ R , β R ) = [ sin θ T cos β T sin θ T sin β T cos θ T ] [ Z x x Z x y Z x z Z y x Z y y Z y z Z z x Z z y Z zz ] [ sin θ R cos β R sin θ R sin β R cos θ R ] ( 6 ) = sin θ T cos β T sin θ R cos β R · Z xx + sin θ T cos β T sin θ R sin β R · Z xy + sin θ T cos β T cos β R · Z xz + sin θ T sin β T sin θ R cos β R · Z yx + sin θ T sin β T sin θ R sin β R · Z yy + sin θ T sin β T cos β R · Z yz + cos θ T sin θ R cos β R · Z zx + cos θ T sin θ R sin β R · Z zy + cos θ T cos β R · Z zz ( 7 )

In Equation 6, the first vector [sin θT cos βT sin θT sin βT cos θT] corresponds to the orientation of the transmitter and the second vector

[ sin θ R cos β R sin θ R sin β R cos θ R ]
corresponds to the orientation of the receiver. Equation 6 may be applied to each receiver using the formulas from Table 1. The general response signal for RX1, RX2 and RX3 are displayed below in equations 8, 9, and 10 respectively.
V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1  (8)
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2  (9)
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3  (10)

By solving equations 8-9, the Green's response tensor may be obtained. In some examples, certain parameters may be already known based on preprocessing of survey information. For example, θT, θR1R2R3, βTR1R2R3 may be known based on tool characterizations. Furthermore, an assumption may be made regarding the relative spacing of each of the receiver subs relative to the transmitter sub. Since the distance between each RX is small compared to each TX-RX spacing, each spacing may be assumed to be ZR1≈ZR2≈ZR3. Additionally, βoff off may be measured on the wellbore tool. The assumptions are applied to Equations 8-10 to transform the equations to Equations 11-13.

V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · Z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · Z zz R 2 ( 11 ) V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · Z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · Z yy R 22 + sin θ T sin β T cos ( n m · 360 ° + β 0 ff ) · Z yz R 2 + cos θ T sin θ R 2 cos β R 2 · Z zx R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · Z zz R 2 ( 12 ) V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · Z zz R 2 ( 13 )

For each of Equations 11-13, n ranges from 1 to m with m corresponding to the number of measurements per turn. As previously discussed, the number of measurements can be any arbitrary number. Equations 11-13 can be combined into a linear system where each component of the m-by-9 matrix Bm×9 may be determined by the coefficients of Equations 11-13. The linear system is illustrated in Equation 14.

[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] ( 14 )

The linear system of Equation 14 is an over-determined linear equation. However, the rank matrix B is 9 so the system may be solved by least square matrix inversion illustrated in Equation 14 where H denotes Hermitian operator.

[ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] ( 15 )

In some examples where the measurements are corrupted by noise, a Tikhonov regularization may be applied to the method. Equation 16 illustrates a Tikhonov regularization applied to Equation 15. The regularization factor λ may be determined by the L-curve method.

[ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]

The methods presented herein were validated by numerical analysis with a signal-noise-ratio of 20 dB. The results of the test are displayed in Table 2.

TABLE 2
Components True Values Test Results
Zxx   0.3118 + 0.1039i   0.3139 + 0.1046i
Zxy   0.0054 − 0.0011i   0.0054 − 0.0011i
Zxz −0.2009 + 0.0502i −0.2025 + 0.0506i
Zyx −0.1385 − 0.1172i −0.1390 − 0.1176i
Zyy   0.0003 + 0.0015i   0.0003 + 0.0015i
Zyz   0.2019 − 0.0505i   0.2035 − 0.0509i
Zzx   0.1322 − 0.1119i   0.1335 − 0.1129i
Zzy −0.0051 − 0.0010i −0.0052 − 0.0010i
Zzz −0.1299 − 0.0162i −0.1301 − 0.0163i

Discussed above are methods and systems for improving current technology. Specifically, improving the quality of electrical images that may be used to determine formation properties. Current problems faced by imaging tools reduce and degrade electrical images. Poor quality electrical images make drilling operations and/or production operations difficult. Operators may make improper decisions based off low quality images. This may extend drilling time, increase costs, reduce production, and/or negatively affect a well. There is a need to produce high quality electrical images that are reliable.

In wellbore imaging operations, certain assumptions may be made about the formation's resistivity trends in 3-D space. These assumptions, referred to as V1-D, may produce inaccurate representations of formation resistivity. As disclosed above, a method which does not make the assumption of V1-D will improve wellbore resistivity images. Specifically, the methods and system may eliminate effects that cause inaccuracies in the imaging data, thereby increasing accuracy of the electrical images.

Accordingly, this disclosure describes apparatus and methods that may relate to subterranean operations. The apparatus, methods, and compositions may further be characterized by one or more of the following statements:

Statement 1. A method comprising:

Statement 2. The method of statement 1 wherein the series of linear equations have the following form for each of the response signals:

V ( θ T , β T , θ R , β R ) = [ sin θ T cos β T sin θ T sin β T cos θ T ] [ Z xx Z xy Z xz Z yx Z yy Z yz Z zx Z zy Z zz ] [ sin θ R cos β R sin θ R sin β R cos θ R ]
where V(θTTRR) is the response signal for receiver coil i, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, θR is a tilt angle of the receiver coil i, βR is an azimuth angle of the receiver coil i, Zxx is a component of the response signal for receiver coil i oriented in an X direction caused by the transmitter coil in the X direction, Zxy is a component of the response signal for receiver coil i oriented in a Y direction caused by the transmitter coil in the X direction, Zxz is a component of the response signal for receiver coil i oriented in an Z direction caused by the transmitter coil in the X direction, Zyx is a component of the response signal for receiver coil i oriented in the X direction caused by the transmitter coil in the Y direction, Zyy is a component of the response signal for receiver coil i oriented in the Y direction caused by the transmitter coil in the Y direction, Zyz is a component of the response signal for receiver coil i oriented in the Z direction caused by the transmitter coil in the Y direction, Zzx is a component of the response signal for receiver coil i oriented in the X direction caused by the transmitter coil in the Z direction, Zzy i is a component of the response signal for receiver coil i oriented in Y direction caused by the transmitter coil in the Z direction, Zzz is a component of the response signal for receiver coil i oriented in the Z direction caused by the transmitter coil in the Z direction, θR is a tilt angle of the receiver coil i, and βR is an azimuth angle of the receiver coil i, wherein receiver coil i is one of the first receiver coil, the second receiver coil, or the third receiver coil.

Statement 3. The method of any preceding statement wherein the series of linear equations has the following form for the response signal corresponding to each receiver coil;

for the first receiver coil:

V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where V1TTR1R1) is the response signal of the first receiver coil, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, ZxxR1 is a component of the first response signal oriented in an X direction caused by the transmitter coil in the X direction, ZxyR1 is a component of the first response signal oriented in a Y direction caused by the transmitter coil in the X direction, ZxzR1 is a component of the first response signal oriented in a Z direction caused by the transmitter coil in the X direction, ZyxR1 is a component of the first response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR1 is a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR1 is a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR1 is a component of the first response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR1 is a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR1 is a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the first receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where V2TTR2R2) is the response signal of the second receiver coil, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR2 is a component of the second response signal oriented in the X direction, ZyxR2 is a component of the second response signal oriented in the Y direction, ZyyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR2 is a tilt angle of the second receiver coil, and βR2 is an azimuth angle of the second receiver coil;
for the third receiver coil:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where V3TTR3R3) is the response signal of the third receiver coil, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the X direction, ZyxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR3 i is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR3 is a tilt angle of the third receiver coil, and βR is an azimuth angle of the third receiver coil.

Statement 4. The method of any preceding statement wherein each azimuth angle is found by the following equations:

azimuth angle of the transmitter coil:

β T = n m · 360
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:
βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver coil:
βR2Toff
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:
βR3ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil.

Statement 5. The method of any preceding statement wherein an assumption of ZR1≈ZR2≈ZR3 is applied to each response signal equation such that:

V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · Z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · Z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · Z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · Z yy R 22 + sin θ T sin β T cos ( n m · 360 ° + β 0 ff ) · Z yz R 2 + cos θ T sin θ R 2 cos β R 2 · Z zx R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · Z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · Z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.

Statement 6. The method of any preceding statement wherein each of the response signals with respect to n is:

[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each of the response signals with respect to n for the first response signal, the second response signal, and the third response signal from n to m.

Statement 7. The method of any previous statement wherein the step solving the series of linear equations comprises solving for a vector defined as:

[ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] .

Statement 8. The method of any previous statement wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:

[ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.

Statement 9. The method of any previous statement wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:

[ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor.

Statement 10. The method of any previous statement wherein λ is determined by an L-curve method.

Statement 11. An apparatus comprising:

a memory coupled to the processor to provide software that configures the processor to solve a series of linear equations defined by an orientation of the transmitter sub in the wellbore, an orientation of the first receiver coil, an orientation of the second receiver coil, and an orientation the third receiver coil, and each of the first response signal, the second response signal, and the third response signal to determine a Green's response tensor.

Statement 12. The apparatus of statement 11 wherein the series of linear equations has the following form for a response signal corresponding to each receiver, the response signal being the first response signal, the second response signal, or the third response signal;

for the first receiver:

V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where V1TTR1R1) is the first response signal, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, ZxxR1 a component of the first response signal oriented in an X direction caused by the transmitter coil in the X direction, ZxyR1 a component of the first response signal oriented in a Y direction caused by the transmitter coil in the X direction, ZxzR1 a component of the first response signal oriented in a Z direction caused by the transmitter coil in the X direction, ZyxR1 a component of the first response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR1 a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR1 i a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR1 a component of the first response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR1 a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR1 i a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where V2TTR2R2) is the second response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the X direction, ZyxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR2 is a tilt angle of the second receiver coil, and βR2 is an azimuth angle of the second receiver coil;
for the third receiver:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where V3TTR3R3) is the third response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the X direction, ZyxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR3 i is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR3 is a tilt angle of the receiver coil, and βR3 is an azimuth angle of the third receiver coil.

Statement 13. The apparatus of any of statements 11-12 wherein each azimuth angle is found by the following equations:

azimuth angle of the transmitter coil:

β T = n m · 360
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:
βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver:
βR2Toff
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:
βR3ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil,
and wherein an assumption of ZR1≈ZR2≈ZR3 is applied to each response signal equation such that:

V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · Z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · Z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · Z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · Z yy R 22 + sin θ T sin β T cos ( n m · 360 ° + β 0 ff ) · Z yz R 2 + cos θ T sin θ R 2 cos β R 2 · Z zx R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · Z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · Z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.

Statement 14. The apparatus of any of statements 11-13 wherein each response signal with respect to n is:

[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ Z xx R 2 Z x y R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each response signal with respect to n for the first receiver, the second receiver, and the third receiver from n to m, wherein the step solving the series of linear equations comprises solving for a vector defined as:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ]
and wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.

Statement 15. The apparatus of any of statements 11-14 wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor and wherein λ is determined by an L-curve method.

Statement 16. A system comprising:

Statement 17. The system of statement 16 wherein the series of linear equations has the following form for a response signal corresponding to each receiver, the response signal being the first response signal, the second response signal, or the third response signal;

for the first receiver coil:

V1TTR1R1)=sin θT cos βT sin θR1 cos βR1·ZxxR1+sin θT cos βT sin θR1 sin βR1·ZxyR1+sin θT cos βT cos βR1·ZxzR1+sin θT sin βT sin θR1 cos βR1·ZyxR1+sin θT sin βT sin θR1 sin βR1·ZyyR1+sin θT sin βT cos βR1·ZyzR1+cos θT sin θR1 cos βR1·ZzxR1+cos θT sin θR1 sin βR1·ZzyR1+cos θT cos βR1·ZzzR1
where V1TTR1R1) is the first response signal, θT is a tilt angle of the transmitter coil, βT is an azimuth angle of the transmitter coil, ZxxR1 is a component of the first response signal oriented in an X direction caused by the transmitter coil in the X direction, ZxyR1 is a component of the first response signal oriented in a Y direction caused by the transmitter coil in the X direction, ZxzR1 is a component of the first response signal oriented in a Z direction caused by the transmitter coil in the X direction, ZyxR1 is a component of the first response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR1 is a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyxR1 is a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR1 is a component of the first response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR1 is a component of the first response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR1 i is a component of the first response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR1 is a tilt angle of the first receiver coil, and βR1 is an azimuth angle of the receiver coil;
for the second receiver coil:
V2TTR2R2)=sin θT cos βT sin θR2 cos βR2·ZxxR2+sin θT cos βT sin θR2 sin βR2·ZxyR2+sin θT cos βT cos βR2·ZxzR2+sin θT sin βT sin θR2 cos βR2·ZyxR2+sin θT sin βT sin θR1 sin βR2·ZyyR22+sin θT sin βT cos βR2·ZyzR2+cos θT sin θR2 cos βR2·ZzxR2+cos θT sin θR2 sin βR2·ZzyR2+cos θT cos βR2·ZzzR2
where V2TTR2R2) is the second response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the X direction, ZyxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR2 is a component of the second response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR2 is a component of the second response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR2 is a component of the second response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR2 is a tilt angle of the second receiver coil, and βR2 is an azimuth angle of the second receiver coil;
for the third receiver:
V3TTR3R1)=sin θT cos βT sin θR3 cos βR3·ZxxR3+sin θT cos βT sin θR3 sin βR3·ZxyR3+sin θT cos βT cos βR3·ZxzR3+sin θT sin βT sin θR3 cos βR3·ZyxR3+sin θT sin βT sin θR3 sin βR3·ZyyR3+sin θT sin βT cos βR3·ZyzR3+cos θT sin θR3 cos βR3·ZzxR3+cos θT sin θR3 sin βR3·ZzyR3+cos θT cos βR3·ZzzR3
where V3TTR3R3) is the third response signal, θT is the tilt angle of the transmitter coil, βT is the azimuth angle of the transmitter coil, ZxxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the X direction, ZxyR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the X direction, ZxzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the X direction, ZyxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Y direction, ZyyR3 i is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Y direction, ZyzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Y direction, ZzxR3 is a component of the third response signal oriented in the X direction caused by the transmitter coil in the Z direction, ZzxR3 is a component of the third response signal oriented in the Y direction caused by the transmitter coil in the Z direction, ZzzR3 is a component of the third response signal oriented in the Z direction caused by the transmitter coil in the Z direction, θR3 is a tilt angle of the receiver coil, and βR3 is an azimuth angle of the third receiver coil.

Statement 18. The apparatus of any of statements 16-17 wherein each azimuth angle is found by the following equations:

azimuth angle of the transmitter coil:

β T = n m · 360
where βT is the azimuth angle of the transmitter coil, m is a number of times the receivers sample per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the transmitter coil:

β T = n m · 360
where βT is the azimuth angle of the transmitter coil, m is a number of times the first receiver coil samples per turn of the apparatus in the wellbore, and n is an index from 1 to m;
azimuth angle of the first receiver coil:
βR1ToffΔ1
where βR1 is the azimuth angle of the first receiver coil, βoff is the difference in azimuthal angle between the transmitter coil and the second receiver coil, βΔ1 is the difference in azimuthal angle between the first receiver coil and the second receiver coil,
azimuth angle of the second receiver:
βR2Toff
where βR2 is the azimuth angle of the second receiver coil,
azimuth angle of the third receiver coil:
βR3ToffΔ2
where βR2 is the azimuth angle of the second receiver coil and βΔ2 the difference in azimuthal angle between the second receiver coil and third receiver coil, and wherein an assumption of ZR1≈ZR2≈ZR3 is applied to each response signal equation such that:

V 1 ( n ) = sin θ T cos β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z xx R 2 + sin θ T cos β T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z xz R 2 + sin θ T sin β T sin θ R 1 cos ( n m · 360 ° + β off + β Δ 1 ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin β R 1 · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 1 ) · Z yz R 2 + cos θ T sin θ R 1 cos ( n m · 360 ° + β off + β Δ1 ) · Z zx R 2 + cos θ T sin θ R 1 sin ( n m · 360 ° + β off + β Δ 1 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 1 ) · Z zz R 2 V 2 ( n ) = sin θ T cos β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z xx R 2 + sin θ T cos β T sin θ R 2 sin ( n m · 360 ° + β off ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off ) · Z xz R 2 + sin θ T sin β T sin θ R 2 cos ( n m · 360 ° + β off ) · Z yx R 2 + sin θ T sin β T sin θ R 1 sin ( n m · 360 ° + β off ) · Z yy R 22 + sin θ T sin β T cos ( n m · 360 ° + β 0 ff ) · Z yz R 2 + cos θ T sin θ R 2 cos β R 2 · Z zx R 2 + cos θ T sin θ R 2 sin ( n m · 360 ° + β off ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off ) · Z zz R 2 V 3 ( n ) = sin θ T cos β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z xx R 2 + sin θ T cos β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z xy R 2 + sin θ T cos β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z xz R 2 + sin θ T sin β T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z yx R 2 + sin θ T sin β T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z yy R 2 + sin θ T sin β T cos ( n m · 360 ° + β off + β Δ 2 ) · Z yz R 2 + cos θ T sin θ R 3 cos ( n m · 360 ° + β off + β Δ 2 ) · Z zx R 2 + cos θ T sin θ R 3 sin ( n m · 360 ° + β off + β Δ 2 ) · Z zy R 2 + cos θ T cos ( n m · 360 ° + β off + β Δ 2 ) · Z zz R 2
where V1(n) is the first response signal with respect to index n, V2(n) is the second response signal with respect to index n, and V3(n) is the third response signal with respect to index n.

Statement 19. The apparatus of any of statements 16-18 wherein each response signal with respect to n is:

[ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ] = B _ _ m × 9 · [ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ]
where Bm×9 is an m by 9 matrix comprising coefficients from each of the first response signal, second response signal, and third response signal with respect to index n to m, wherein the step solving the series of linear equations comprises solving for a vector defined as:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ]
and wherein the vector is solved for using least square matrix inversion using a Hermitian operator such that:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where H is the Hermitian operator.

Statement 20. The apparatus of any of statements 16-18 wherein the vector is solved for using least square matrix inversion using a Hermitian operator and Tikhonov regularization is applied such that:

[ Z xx R 2 Z xy R 2 Z x z R 2 Z y x R 2 Z y y R 2 Z y z R 2 Z z x R 2 Z z y R 2 Z zz R 2 ] = ( B _ _ m × 9 H · B _ _ m × 9 + λ I _ _ ) - 1 · B _ _ m × 9 H · [ V 1 ( n = 1 ) V 2 ( n = 1 ) V 3 ( n = 1 ) V 1 ( n = m ) V 2 ( n = m ) V 3 ( n = m ) ]
where λ is a regularization factor and wherein λ is determined by an L-curve method.

The preceding description provides various examples of the wellbore tools and methods of use disclosed herein which may contain different method steps and alternative combinations of components. It should be understood that, although individual examples may be discussed herein, the present disclosure covers all combinations of the disclosed examples, including, without limitation, the different component combinations, method step combinations, and properties of the system. It should be understood that the compositions and methods are described in terms of “comprising,” “containing,” or “including” various components or steps, the compositions and methods can also “consist essentially of” or “consist of” the various components and steps. Moreover, the indefinite articles “a” or “an,” as used in the claims, are defined herein to mean one or more than one of the element that it introduces.

For the sake of brevity, only certain ranges are explicitly disclosed herein. However, ranges from any lower limit may be combined with any upper limit to recite a range not explicitly recited, as well as, ranges from any lower limit may be combined with any other lower limit to recite a range not explicitly recited, in the same way, ranges from any upper limit may be combined with any other upper limit to recite a range not explicitly recited. Additionally, whenever a numerical range with a lower limit and an upper limit is disclosed, any number and any included range falling within the range are specifically disclosed. In particular, every range of values (of the form, “from about a to about b,” or, equivalently, “from approximately a to b,” or, equivalently, “from approximately a-b”) disclosed herein is to be understood to set forth every number and range encompassed within the broader range of values even if not explicitly recited. Thus, every point or individual value may serve as its own lower or upper limit combined with any other point or individual value or any other lower or upper limit, to recite a range not explicitly recited.

Fan, Yijing, Pan, Li, Wu, Hsu-Hsiang

Patent Priority Assignee Title
Patent Priority Assignee Title
7545145, May 03 2006 Baker Hughes Incorporated Method and apparatus for tensorial micro-resistivity imaging in oil-based muds
7756641, Dec 20 2004 Schlumberger Technology Corporation Computation of sensitivity for resistivity measurements
7987074, Mar 08 2006 ExxonMobil Upstream Research Company Efficient computation method for electromagnetic modeling
8274289, Dec 15 2006 Halliburton Energy Services, Inc Antenna coupling component measurement tool having rotating antenna configuration
20080215243,
20100188089,
20100277177,
20110140701,
20120026314,
20150185344,
20160002977,
20160116616,
20160187525,
20170096887,
20170261637,
20200018149,
GB2409723,
WO2013086270,
////
Executed onAssignorAssigneeConveyanceFrameReelDoc
May 14 2018PAN, LIHalliburton Energy Services, IncASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0498020950 pdf
May 14 2018FAN, YIJINGHalliburton Energy Services, IncASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0498020950 pdf
Jun 07 2018Halliburton Energy Services, Inc.(assignment on the face of the patent)
Jun 08 2018WU, HSU-HSIANGHalliburton Energy Services, IncASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0498020950 pdf
Date Maintenance Fee Events
Jul 19 2019BIG: Entity status set to Undiscounted (note the period is included in the code).
Jun 20 2024M1551: Payment of Maintenance Fee, 4th Year, Large Entity.


Date Maintenance Schedule
Feb 23 20244 years fee payment window open
Aug 23 20246 months grace period start (w surcharge)
Feb 23 2025patent expiry (for year 4)
Feb 23 20272 years to revive unintentionally abandoned end. (for year 4)
Feb 23 20288 years fee payment window open
Aug 23 20286 months grace period start (w surcharge)
Feb 23 2029patent expiry (for year 8)
Feb 23 20312 years to revive unintentionally abandoned end. (for year 8)
Feb 23 203212 years fee payment window open
Aug 23 20326 months grace period start (w surcharge)
Feb 23 2033patent expiry (for year 12)
Feb 23 20352 years to revive unintentionally abandoned end. (for year 12)