# Simulation of Monopole and Multipole Seismoelectric Logging

Authors: Zhiwen Cui; Jinxia Liu; Yujun Zhang; Kexie Wang; Hengshan Hu
Journal: Advances in Acoustics and Vibration
 (2011)

---

## Abstract

In a fluid-saturated porous formation, acoustics and electromagnetic waves are coupled based on Pride seismoelectric theory. An exact treatment of the nonaxisymmetric seismoelectric field excited by acoustic multipole sources is presented. The frequency wavenumber domain representations of the acoustic field and associated seismoelectric field due to acoustic multipole sources are formulated. The full waveforms of acoustic waves and electric and magnetic fields in the time domain propagation in borehole are simulated by using discrete wave number integration, and frequency versus axial-wave number responses are presented and analyzed.


---

## Body

## 1. Introduction

The study of wave propagation in a fluid-saturated porous medium is of considerable interest in acoustics and geophysics due to its important applications in various technical and engineering processes. The investigation of wave propagation in fluid-saturated porous media was early developed by Biot [1, 2]. One of the major findings of Biot’s theory was that there is a compressional slow wave in a fluid-saturated porous medium. The first clear experimental observation of this slow wave was reported by Plona [3]. Biot predicted the slow waves should have an important bearing on electrokinetic effect [4]. This predication has been quantitatively confirmed by Pride [5] and Hu [6]. Elastic waves propagating in fluid-saturated porous media generate a movement of the ions in the pore fluid. Such movement induces an electromagnetic (EM) field. Thompson and Gist [7] have made field measurement clearly demonstrating that seismic waves can induce electromagnetic disturbances in saturated sediments. Pride [8] derived the governing equations for the coupled acoustics and electromagnetic waves in fluid-saturated porous media. Pride and Haartsen [9] analyzed the basic properties of seismoelectric waves.The electric field induced by elastic waves is weak and attenuates in propagation. In order to detect the seismoelectric signal effectively, Haartsen and Pride [10] suggested measuring vertical electroseismic profile. Mikhailov et al. [11] measured the electric field converted from low frequency Stoneley waves in a borehole and made theoretical analysis. Seismoelectric logging method has been proposed to detect deep target formation. The advantage of seismoelectric logging is the distance both from transmitter to target formation and from the transmitter to the receiver is small, and signals can be received with relative high amplitude. Zhu et al. [12–14] made laboratory experiments and observed the seismoelectric conversion in model wells, and their experimental results confirm that seismoelectric logging could be a new borehole logging technique. Seismoelectric logging methods based on the excitation and reception of axisymmetric wave phenomena in a fluid-filled borehole embedded in permeable formations have been investigated [6, 15]. In the present paper, we consider the nonaxisymmetric seismoelectric field excited by acoustic multipole sources based on the entire Pride seismoelectric theory [9]. We theoretically formulate and numerically simulate both axisymmetric and nonaxially symmetric waves in a fluid-filled borehole surrounded by a homogeneous fluid-saturated porous formation. The frequency versus axial-wave number responses are presented, and time-domain transient waveforms of acoustic waves and electromagnetic waves propagation in borehole are simulated and compared. The current model may be an important tool in understanding and predicting borehole seismoelectric phenomena. This paper may be of interest for researchers in the field referred to above as well as experimental wave propagation in fluid-saturated porous media.


## 2. Theoretical Formulations

In this section we formulate a frequency wavenumber integral representation for the acoustic, electric, and magnetic fields excited by a multipole acoustic source in the borehole fluid. Consider a fluid-filled borehole of radiusa embedded in linear isotropic fluid-saturated porous medium and acoustic multipole sources in the borehole. In the borehole, the acoustic and electromagnetic fields are not coupled, but when the acoustic source in the borehole is excited, according to Pride’s theory [8, 9], the acoustic waves in porous formation will generate accompanying electric and magnetic fields, and acoustic wave induces radiating EM wave at borehole wall (discontinuous interface). We can receive electric or magnetic fields signal inside the borehole due to the boundary condition at the borehole wall.


### 2.1. The Acoustic Field in the Borehole Fluid

We consider acoustic multipole sources inside the borehole. The Dipole source is constructed from two point sources of opposite sign placed close together in the same horizontal plane. Equivalently, the dipole can be viewed as a point force oriented horizontally and pointed at an azimuthal angle ofθ0. Similarly, a quadrupole source can be viewed as two closely spaced dipole sources pointing in opposite directions. Higher-order multipoles have similar interpretations. Following Kurkjian and Chang [16], with emphasizing the radius of a multipole source being very small and |ηfr0|≪1, the representation in the frequency-wavenumber domain of the displacement potential associated with n acoustic multipole sources (r>r0) is given by(1)Φ0(r,θ,kz,ω)=γεnKn(ηfr)cos⁡n(θ-θ0),
where γ=V0(ω)(ηfr0/2)n/n!,  εn is Neumann’s factor (εn=1 for n=0; εn=2 for n=otherwise), V0(ω) is the source spectrum, and ηf2=kz2-ω2/cf2 is the radial wavenumber with respect to the acoustic velocity of bore fluid cf.Kn(x) is the nth order modified Bessel function of the second kind.The total acoustic wave field in the borehole is the sum of the source contribution and the reflected acoustic field. The radiation conditions prescribe that the reflected field should be finite at the axis of the borehole, thus the reflected waves displacement potential associated withn acoustic multipole sources is(2)Φ1(r,θ,kz,ω)=γAnIn(ηfr)cos⁡n(θ-θ0),
where In(x) is the nth-order modified Bessel function of the second kind. Thus, the total field inside the borehole in the frequnecy-wavenumber domain is given by(3)Φ(r,θ,kz,ω)=γ[AnIn(ηfr)+εnKn(ηfr)]cos⁡n(θ-θ0),
the reflection coefficients An can be derived by the boundary conditions at the borehole wall. The radial component of the displacement field in the borehole fluid is given by ur=∂Φ/∂r and the radial component of normal stress in the borehole fluid is equal to the negative of borehole fluid pressure, that is,
(4)τrr=-P(r,θ,kz,ω)=-ρfω2Φ.


### 2.2.  The Electromagnetic Field in the Borehole Fluid

In cylindrical coordinate system(r,θ,z), it is convenient for us to take the z-component of vector E and H, and they satisfy the Homheltz-type wave equations without electric current sources in the borehole,(5)∇2Ez+ke2Ez=0,(6)∇2Hz+ke2Hz=0,
where ke2=ω2μεb,εb=εf(1+iσf/ωεf),εf is the borehole fluid’s electrical permittivity, σf  is the borehole fluid’s electrical conductivity, andμ is the magnetic permeability. For really fluid and all porous media of interest it can be assumed that μ=μ0. In borehole fluid and fluid-saturated porous media we take μ=μ0 throughout the paper.Considering the radiation conditions, that prescribe the reflected electromagnetic field should be finite at the axis of the borehole, we can take the solutions forEz and Hz as (7)Ez(r,θ,kz,ω)=γAnEIn(ηfr)cos⁡n(θ-θ0),(8)Hz(r,θ,kz,ω)=γAnHIn(ηfr)sinn(θ-θ0),
where ηe2=kz2-ke2. As we will see, in order to satisfy the boundary conditions of continuity of tangential electromagnetic fields at borehole wall r=a, the solutions in general require a combination of both TE and TM fields. The modes with Ez and Hz components present are hybrid modes. Transverse field components are obtained from Maxwell’s equations, (9)Eθ=iηe2(ωμ∂Hz∂r-kzrEz∂θ),(10)Hθ=-iηe2(ωεb∂Ez∂r+kzr∂Hz∂θ),(11)Er=-iηe2(ωμr∂Hz∂θ+kz∂Ez∂r),(12)Hr=iηe2(ωεbr∂Ez∂θ-kz∂Hz∂r).Reflected electromagnetic field inside the borehole in the time-space domain is given by(13)Π(r,θ,z,t)=12π∬∞∞Π(r,θ,kz,ω)exp⁡[i(kzz-ωt)]dkzdω;
here we use Π to denote any component of electric or magnetic field given by (7)–(12). The reflection coefficients AnE, AnH can be derived by the boundary conditions at the borehole wall.


### 2.3. Seismoelectric Waves in the Porous Formation

For the problem of modeling the propagation of coupled electromagnetic and mechanical disturbances in an isotropic-porous material, Pride [8] has derived equations that control such “seismoelectric” phenomena. According to Pride’s theory, without the applied force and electric current sources existing, and assuming an time dependence exp(-iwt), we may write the equations for the coupled electromagnetic and acoustics in macroscopically homogenous, isotropic, fluid saturated porous media as follows(14)∇·τ=-ω2(ρu+ρfw),(15)τ=(H-2Gb)(∇·u)I+C(∇·w)I+Gb(∇u+∇uT),(16)-pf=C(∇·u)+M(∇·w),(17)-iωw=(-∇pf+ρfω2u)k(ω)η+L(ω)E,(18)J=(-∇p+pfω2u)L(ω)+σ(ω)E,(19)∇×E=iωμH,(20)∇×H=-iωεE+J,Those governing equations are the Biot equations for porous media acoustic along with the Maxwell equations for the electric and magnetic fieldsE and H. Here τ is the bulk stress in the porous medium, pf is the pressure in the pore fluid, u is the displacement in the solid, and w is the relative fliud-solid motion. The symbol ρ denotes the bulk density of the porous medium, ρ=(1-ϕ)ρs+ϕρf,ϕis the porosity of the medium, ρs is the solid density, and ρf denotes the fliud density. Equations (17) and (18) are in the form of Darcy law and Ohm law, through which acoustic and electromagnetic fields are coupled, where J is the electric-current density and -iωw is the Darcy filtration velocity. Where H, C, M, and Gb are four moduli of isotropic porous media. As the most important coefficient is set to zero, Pride’s equations will be separated into Biot’s equations for elastic field and Maxwell equations for electromagnetic field. Here σ(ω) is the frequency-dependent electrical conductivity of the medium, k(ω) is the dynamic permeability, η is the fluid dynamic viscosity, and L(ω) is the frequency-dependent electrokinetic coupling coefficient.The expressions for σ(ω), k(ω), and L(ω) are given in Pride [8].The solution to the coupled equations of seismoelectric waves motions can be obtained by separating the basic fieldu, w, E into its compressional, and vertically and horizontally polarized shear components, and written in terms of scalar potential functions as(21)u=∇(Φpf+Φps)+∇×(Ψsh+Ψem)ez+∇×∇×(Γsh+Γem)ez,(22)w=∇(αpfΦpf+αpsΦps)+∇×(αshΨsh+αemΨem)ez+∇×∇×(αshΓsh+αemΓem)ez,(23)E=∇(βpfΦpf+βpsΦps)+∇×(βshΨsh+βemΨem)ez+∇×∇×(βshΓsh+βemΓem)ez,
where Φj is the compressional wave potential, ez is the unit vector in the z direction, Ψj is the horizontal polarized shear potential, and Γj is vertically polarized shear wave potential. Each of these potentials satisfy Helmholtz-type wave equation. The factors αpf, αps, αsh, αem, βpf, βps, βsh, and βem as well as the formulas of the body wave velocities can be found in [9, 15, 17].Outside the borehole, because the radiation conditions prescribe no incoming wave existing at the position wherer tends to infinity, the solutions in the frequency-axial wavenumber with respect to an nth-order multipole source are, similar to only acoustic fields case [18, 19],
(24)Φj(r,θ,kz,ω)=γBnjKn(ηjr)cos⁡n(θ-θ0),(25)Ψj(r,θ,kz,ω)=γCnjKn(ηjr)cos⁡n(θ-θ0),(26)Γj(r,θ,kz,ω)=γDnjKn(ηjr)cos⁡n(θ-θ0),
whereηj2=kz2-kj2(kj2=ω2sj2) are the radical wavenumbers of the seismoelectric waves whose complex slowness sj(j=pf,ps,sh,em) expressions can be found in [9, 15, 17]. Here, Bnj, Cnj, and Dnj which characterize the outgoing seismoelectric waves, are functions of the axial wavemumber kz and the angular frequency ω. Once u, w, E are known, all other quantities such as stress tensor, pore fluid pressure, and magnetic vector can be deduced by (15)-(16) and (19). The components of seisoelectric filed in the formation can be written in terms of the potentials using (23) and (26). The detailed expressions of those formulations can be found in [17].


### 2.4. Boundary Conditions and Resolution

The unknown coefficientsAn, AnE, AnH, Bnpf, Bnps, Cnsh, Cnem, Dnsh, andDnem which appear in the expressions for the field potentials are determined by the boundary conditions. The boundary conditions of seismoelectric waves for a cylindrical interface (at r=a) between a fluid and a fluid-saturated porous media for the case of permeable interface are(27)urI=urII+wrII,-PI=τrrII,0=τrzII;0=τrθII,PI=pfII;EθI=EθII,EzI=EzII,HθI=HθII,HzI=HzII,
where the superscript I and II represent the media inside and outside the borehole, respectively. In the present paper, we do not consider the possibility of a sealed or partially sealed boundary condition. With the expressions of wr, ur, τrr, τrz, τrθ, Eθ,Ez, Hθ, Hz deduced from (14)–(26), (7) –(10), and urI=∂(Φ0+Φ1)/∂r, PI=ρfω2(Φ0+Φ1) boundaries (27) may be written in a matrix form as(28)M̅A=B,
where M̅ are 9×9 complex matrices whose elements are given in the [17] and (29)A={An,Bnpf,Bnps,Cnsh,Dnsh,Cnem,Dnem,AnE,AnH}T,(30)B={b1,b2,0,0,b5,0,0,0,0}T,
where b1=εn[ηfKn+1(ηfa)-(n/a)Kn(ηfa)], b2=b5=-ρfω2εnKn(ηfa).


## 2.1. The Acoustic Field in the Borehole Fluid

We consider acoustic multipole sources inside the borehole. The Dipole source is constructed from two point sources of opposite sign placed close together in the same horizontal plane. Equivalently, the dipole can be viewed as a point force oriented horizontally and pointed at an azimuthal angle ofθ0. Similarly, a quadrupole source can be viewed as two closely spaced dipole sources pointing in opposite directions. Higher-order multipoles have similar interpretations. Following Kurkjian and Chang [16], with emphasizing the radius of a multipole source being very small and |ηfr0|≪1, the representation in the frequency-wavenumber domain of the displacement potential associated with n acoustic multipole sources (r>r0) is given by(1)Φ0(r,θ,kz,ω)=γεnKn(ηfr)cos⁡n(θ-θ0),
where γ=V0(ω)(ηfr0/2)n/n!,  εn is Neumann’s factor (εn=1 for n=0; εn=2 for n=otherwise), V0(ω) is the source spectrum, and ηf2=kz2-ω2/cf2 is the radial wavenumber with respect to the acoustic velocity of bore fluid cf.Kn(x) is the nth order modified Bessel function of the second kind.The total acoustic wave field in the borehole is the sum of the source contribution and the reflected acoustic field. The radiation conditions prescribe that the reflected field should be finite at the axis of the borehole, thus the reflected waves displacement potential associated withn acoustic multipole sources is(2)Φ1(r,θ,kz,ω)=γAnIn(ηfr)cos⁡n(θ-θ0),
where In(x) is the nth-order modified Bessel function of the second kind. Thus, the total field inside the borehole in the frequnecy-wavenumber domain is given by(3)Φ(r,θ,kz,ω)=γ[AnIn(ηfr)+εnKn(ηfr)]cos⁡n(θ-θ0),
the reflection coefficients An can be derived by the boundary conditions at the borehole wall. The radial component of the displacement field in the borehole fluid is given by ur=∂Φ/∂r and the radial component of normal stress in the borehole fluid is equal to the negative of borehole fluid pressure, that is,
(4)τrr=-P(r,θ,kz,ω)=-ρfω2Φ.


## 2.2.  The Electromagnetic Field in the Borehole Fluid

In cylindrical coordinate system(r,θ,z), it is convenient for us to take the z-component of vector E and H, and they satisfy the Homheltz-type wave equations without electric current sources in the borehole,(5)∇2Ez+ke2Ez=0,(6)∇2Hz+ke2Hz=0,
where ke2=ω2μεb,εb=εf(1+iσf/ωεf),εf is the borehole fluid’s electrical permittivity, σf  is the borehole fluid’s electrical conductivity, andμ is the magnetic permeability. For really fluid and all porous media of interest it can be assumed that μ=μ0. In borehole fluid and fluid-saturated porous media we take μ=μ0 throughout the paper.Considering the radiation conditions, that prescribe the reflected electromagnetic field should be finite at the axis of the borehole, we can take the solutions forEz and Hz as (7)Ez(r,θ,kz,ω)=γAnEIn(ηfr)cos⁡n(θ-θ0),(8)Hz(r,θ,kz,ω)=γAnHIn(ηfr)sinn(θ-θ0),
where ηe2=kz2-ke2. As we will see, in order to satisfy the boundary conditions of continuity of tangential electromagnetic fields at borehole wall r=a, the solutions in general require a combination of both TE and TM fields. The modes with Ez and Hz components present are hybrid modes. Transverse field components are obtained from Maxwell’s equations, (9)Eθ=iηe2(ωμ∂Hz∂r-kzrEz∂θ),(10)Hθ=-iηe2(ωεb∂Ez∂r+kzr∂Hz∂θ),(11)Er=-iηe2(ωμr∂Hz∂θ+kz∂Ez∂r),(12)Hr=iηe2(ωεbr∂Ez∂θ-kz∂Hz∂r).Reflected electromagnetic field inside the borehole in the time-space domain is given by(13)Π(r,θ,z,t)=12π∬∞∞Π(r,θ,kz,ω)exp⁡[i(kzz-ωt)]dkzdω;
here we use Π to denote any component of electric or magnetic field given by (7)–(12). The reflection coefficients AnE, AnH can be derived by the boundary conditions at the borehole wall.


## 2.3. Seismoelectric Waves in the Porous Formation

For the problem of modeling the propagation of coupled electromagnetic and mechanical disturbances in an isotropic-porous material, Pride [8] has derived equations that control such “seismoelectric” phenomena. According to Pride’s theory, without the applied force and electric current sources existing, and assuming an time dependence exp(-iwt), we may write the equations for the coupled electromagnetic and acoustics in macroscopically homogenous, isotropic, fluid saturated porous media as follows(14)∇·τ=-ω2(ρu+ρfw),(15)τ=(H-2Gb)(∇·u)I+C(∇·w)I+Gb(∇u+∇uT),(16)-pf=C(∇·u)+M(∇·w),(17)-iωw=(-∇pf+ρfω2u)k(ω)η+L(ω)E,(18)J=(-∇p+pfω2u)L(ω)+σ(ω)E,(19)∇×E=iωμH,(20)∇×H=-iωεE+J,Those governing equations are the Biot equations for porous media acoustic along with the Maxwell equations for the electric and magnetic fieldsE and H. Here τ is the bulk stress in the porous medium, pf is the pressure in the pore fluid, u is the displacement in the solid, and w is the relative fliud-solid motion. The symbol ρ denotes the bulk density of the porous medium, ρ=(1-ϕ)ρs+ϕρf,ϕis the porosity of the medium, ρs is the solid density, and ρf denotes the fliud density. Equations (17) and (18) are in the form of Darcy law and Ohm law, through which acoustic and electromagnetic fields are coupled, where J is the electric-current density and -iωw is the Darcy filtration velocity. Where H, C, M, and Gb are four moduli of isotropic porous media. As the most important coefficient is set to zero, Pride’s equations will be separated into Biot’s equations for elastic field and Maxwell equations for electromagnetic field. Here σ(ω) is the frequency-dependent electrical conductivity of the medium, k(ω) is the dynamic permeability, η is the fluid dynamic viscosity, and L(ω) is the frequency-dependent electrokinetic coupling coefficient.The expressions for σ(ω), k(ω), and L(ω) are given in Pride [8].The solution to the coupled equations of seismoelectric waves motions can be obtained by separating the basic fieldu, w, E into its compressional, and vertically and horizontally polarized shear components, and written in terms of scalar potential functions as(21)u=∇(Φpf+Φps)+∇×(Ψsh+Ψem)ez+∇×∇×(Γsh+Γem)ez,(22)w=∇(αpfΦpf+αpsΦps)+∇×(αshΨsh+αemΨem)ez+∇×∇×(αshΓsh+αemΓem)ez,(23)E=∇(βpfΦpf+βpsΦps)+∇×(βshΨsh+βemΨem)ez+∇×∇×(βshΓsh+βemΓem)ez,
where Φj is the compressional wave potential, ez is the unit vector in the z direction, Ψj is the horizontal polarized shear potential, and Γj is vertically polarized shear wave potential. Each of these potentials satisfy Helmholtz-type wave equation. The factors αpf, αps, αsh, αem, βpf, βps, βsh, and βem as well as the formulas of the body wave velocities can be found in [9, 15, 17].Outside the borehole, because the radiation conditions prescribe no incoming wave existing at the position wherer tends to infinity, the solutions in the frequency-axial wavenumber with respect to an nth-order multipole source are, similar to only acoustic fields case [18, 19],
(24)Φj(r,θ,kz,ω)=γBnjKn(ηjr)cos⁡n(θ-θ0),(25)Ψj(r,θ,kz,ω)=γCnjKn(ηjr)cos⁡n(θ-θ0),(26)Γj(r,θ,kz,ω)=γDnjKn(ηjr)cos⁡n(θ-θ0),
whereηj2=kz2-kj2(kj2=ω2sj2) are the radical wavenumbers of the seismoelectric waves whose complex slowness sj(j=pf,ps,sh,em) expressions can be found in [9, 15, 17]. Here, Bnj, Cnj, and Dnj which characterize the outgoing seismoelectric waves, are functions of the axial wavemumber kz and the angular frequency ω. Once u, w, E are known, all other quantities such as stress tensor, pore fluid pressure, and magnetic vector can be deduced by (15)-(16) and (19). The components of seisoelectric filed in the formation can be written in terms of the potentials using (23) and (26). The detailed expressions of those formulations can be found in [17].


## 2.4. Boundary Conditions and Resolution

The unknown coefficientsAn, AnE, AnH, Bnpf, Bnps, Cnsh, Cnem, Dnsh, andDnem which appear in the expressions for the field potentials are determined by the boundary conditions. The boundary conditions of seismoelectric waves for a cylindrical interface (at r=a) between a fluid and a fluid-saturated porous media for the case of permeable interface are(27)urI=urII+wrII,-PI=τrrII,0=τrzII;0=τrθII,PI=pfII;EθI=EθII,EzI=EzII,HθI=HθII,HzI=HzII,
where the superscript I and II represent the media inside and outside the borehole, respectively. In the present paper, we do not consider the possibility of a sealed or partially sealed boundary condition. With the expressions of wr, ur, τrr, τrz, τrθ, Eθ,Ez, Hθ, Hz deduced from (14)–(26), (7) –(10), and urI=∂(Φ0+Φ1)/∂r, PI=ρfω2(Φ0+Φ1) boundaries (27) may be written in a matrix form as(28)M̅A=B,
where M̅ are 9×9 complex matrices whose elements are given in the [17] and (29)A={An,Bnpf,Bnps,Cnsh,Dnsh,Cnem,Dnem,AnE,AnH}T,(30)B={b1,b2,0,0,b5,0,0,0,0}T,
where b1=εn[ηfKn+1(ηfa)-(n/a)Kn(ηfa)], b2=b5=-ρfω2εnKn(ηfa).


## 3. Simulations

In this section, we present synthetic acoustic waves and electromagnetic field in the borehole excited by acoustic monopole, dipole, and quadrupole sources. The formation in our numerical example is permeable fluid saturated porous media and the medium input parameters describing the porous formation are given in Table1. The borehole radius a is 0.1 m. The same medium parameters and source function as in [20]. The source is a cosine pulse [19]. If set acoustic transducer, electrode, and magnetometer in the borehole fluid, one could measure acoustic, electric and magnetic fields, respectively. For acoustic logging, only the axisymmetric modes contribute to the pressure on the axis [21]. However, for electroacoustic logging, beside axisymmetric modes, dipole modes of  nonaxisymmetric modes also contribute to the electric and magnetic field on the axis [20, 22]. For example, we only select some position (0.05, 0, z) off the axis to record the full waveforms. The spectral response is calculated similar to [23].Table 1
Medium properties used at input for numerical calculations.
PropertyFastformationSlowformationBoreholefluidPorosityϕ0.20.2—Permeabilityκ0, μm21.01.0—Tortuosityα∞3.03.0—Bulk modulus solidKs, Pa35.70×10933.40×109—Bulk modulus fluidKf, Pa2.25×1092.25×1092.25×109Frame bulk modulusKb, Pa14.39×1092.61×109—Frame shear modulusGb, Pa13.99×1092.83×109—Density solidρs, Kg/m326502650—Density fluidρf, Kg/m3100010001000Fluid viscosityη, Pa.s0.0010.0010.001SalinityC0, mol/L0.010.010.01Permittivity of fluidεf80ε080ε080ε0Permittivity of solidεs4ε04ε0—


### 3.1. Fast Formation

Figure1 shows spectral response and synthetic waveforms of acoustic waves and electric field for fast formation generated by an acoustic monopole source. In our calculation, the source is set to be pressure source, and the center frequency is 6 kHz. The pressure amplitude is 10 MPa measured at 1 cm away from the source center. The numerical result of borehole pressure waveforms is on the order of those by Kurkjian [16]. From Figures 1(b) and 1(d), we can find there are four wave groups in electric field Ez  that is different from three in acoustic waveforms. A clear seismoelectric response A-A group for the electric field is relatively weak but visible in full waveforms, which were recorded at all receivers simultaneously. And the following groups for the electric field having apparent velocities are comparable to those of borehole acoustic waves. Because of the nonequilibrium of free charges in the pore fluid at the borehole wall, that is, at the interface, the acoustic wave generates a propagating electromagnetic wave (EM wave) [9]. The almost same arrival time of small wave components in full waveforms indicates that they propagate with higher speed than that of any acoustic waves.Spectral representation and time response of borehole pressure and electric fields generated by monopole source. (a) and (c) are spectral of pressure and electric fields, respectively. (b) and (d) are time responses of pressure and electric fields, respectively.
(a)(b)(c)(d)In the spectral representation, the response allows the identification of the pesudo-Rayleigh wave (pR.1), which has an about cut-off frequency of 7 kHz. The frequency of the excited source pulse is enough to reveal this wave mode. The excitation is stronger at low frequency for Stoneley wave without cut-off frequency. Spectral representation of electric field (Figure1(c)) Ez shows that the borehole seismoelectric conversion dominates at low frequency. Figures 1(c) and 1(d) shows that P wave response may be identified easily in both spectral representation and waveforms of electric field. In this case the seismoelectric conversion efficiency of P wave is biggest, Stoneley wave is second.Figure2 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for fast formation generated by an acoustic dipole source. In the spectral representations of borehole pressure (Figure 2(a)), the two flexual modes obviously appear in the frequency range 0–20 kHz; however, in the spectral representations of magnetic field, the lowest mode is dominating. From Figures 2(b) and 2(d), the waveforms of borehole pressure with 3 kHz center frequency dominates by flexural guided waves; however, that of magnetic field has two groups with different velocity. Again, from Figure 2(d) we see a relatively weak but visible A-A group; it arrives simultaneously on all receivers. And the following groups having apparent velocities of acoustic flexural waves velocities of borehole acoustic waves.Spectral representation and time response of borehole pressure and magnetic fields for fast formation generated by dipole source. (a) and (c) are spectral of pressure and magnetic field respectively. (b) and (d) are time responses of pressure and magnetic field respectively.
(a)(b)(c)(d)Figure3 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr  for fast formation generated by an acoustic quadrupole source. In the spectral representations of borehole pressure (Figure 3(a)), the two screw modes obviously appear in the range of 0–20 kHz; however, in the spectral representations of magnetic field, the lowest mode is dominating. From Figures 3(b) and 3(d), the waveforms of acoustic waves quadrupole sources with the center frequency 5 kHz mainly dominated by screw waves and that of magnetic field consists of groups having apparent velocities of acoustic screw waves. And we can find almost no EM waves in magnetic field full waveforms.Spectral representation and time response of borehole pressure and magnetic fields for fast formation generated by quadrupole source. (a) and (c) are spectral of pressure and magnetic fields respectively. (b) and (d) are time responses of pressure and magnetic fields respectively.
(a)(b)(c)(d)


### 3.2. Slow Formation

Figure4 shows spectral response and synthetic waveforms for acoustic waves and electric field for slow formation generated by an acoustic monopole source. From Figures 4(b) and 4(d), we can find there are three wave groups in electric field Ez different from two in acoustic waveforms. A clear seismoelectric response A-A group for the electric field Ez is relatively weak but visible in full waveforms, which were recorded at all receivers simultaneously. And the following groups for the electric field having apparent velocities are comparable to those of borehole acoustic waves.Spectral representation and time response of borehole pressure and electric fields for slow formation generated by monopole source. (a) and (c) are spectral of pressure and electric fields, respectively. (b) and (d) are time responses of pressure and electric fields, respectively.
(a)(b)(c)(d)In the spectral representation, the pesudo-Rayleigh wave disappears and the response allows the identification of P wave and Stoneley wave. The excitation is stronger at low frequency for Stoneley wave for slow formation same as fast formation. Spectral representation of electric field (Figure4(c)) Ez shows that the borehole seismoelectric conversion dominates at low frequency. Figures 4(c) and 4(d) shows that P wave response may be identified easily in both spectral representation and waveforms of electric field. Different from fast formation case, the seismoelectric conversion efficiency of Stoneley wave is biggest, P wave is second.Figure5 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for slow formation generated by an acoustic dipole source.  In the spectral representations of acoustic (Figure 5 (a)), one flexural mode appears in the frequency range 0–20 kHz, and in the spectral representations of magnetic field (Figure 5 (b)), the flexural mode is dominating. From Figures 5(b) and 5(d), the waveforms of borehole pressure with 1 kHz center frequency dominates by flexural guided waves; however, that of magnetic field has two groups with different velocity. Again, from Figure 5(d) we see a relatively weak but visible A-A group; it arrives simultaneously on all receivers. And the following group has apparent velocities of acoustic flexural waves velocities of borehole acoustic waves. Again, it shows the borehole seismoelectric conversion efficiency is strong at low frequency range.Spectral representation and time response of borehole pressure and magnetic fields for slow formation generated by dipole source. (a) and (c) are spectral of pressure and magnetic fields, respectively. (b) and (d) are time responses of pressure and magnetic fields, respectively.
(a)(b)(c)(d)Figure6 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for slow formation generated by an acoustic quadrupole source. In the spectral representations of acoustic (Figure 6(a)), a screw mode appears in the range of 0–10 kHz; however, in the spectral representations of magnetic field, this mode is more evident. From Figures 6(b) and 6(d), the waveforms of acoustic waves excited by quadrupole source with the center frequency 3 kHz mainly dominated by screw waves and that of magnetic field consists of groups having apparent velocities of acoustic screw waves. We almost cannot find the EM wave in the full waveforms of magnetic field Hr.Spectral representation and time response of borehole pressure and magnetic fields for slow formation generated by quadrupole source. (a) and (c) are spectral of pressure and magnetic fields, respectively. (b) and (d) are time responses of pressure and magnetic fields, respectively.
(a)(b)(c)(d)


## 3.1. Fast Formation

Figure1 shows spectral response and synthetic waveforms of acoustic waves and electric field for fast formation generated by an acoustic monopole source. In our calculation, the source is set to be pressure source, and the center frequency is 6 kHz. The pressure amplitude is 10 MPa measured at 1 cm away from the source center. The numerical result of borehole pressure waveforms is on the order of those by Kurkjian [16]. From Figures 1(b) and 1(d), we can find there are four wave groups in electric field Ez  that is different from three in acoustic waveforms. A clear seismoelectric response A-A group for the electric field is relatively weak but visible in full waveforms, which were recorded at all receivers simultaneously. And the following groups for the electric field having apparent velocities are comparable to those of borehole acoustic waves. Because of the nonequilibrium of free charges in the pore fluid at the borehole wall, that is, at the interface, the acoustic wave generates a propagating electromagnetic wave (EM wave) [9]. The almost same arrival time of small wave components in full waveforms indicates that they propagate with higher speed than that of any acoustic waves.Spectral representation and time response of borehole pressure and electric fields generated by monopole source. (a) and (c) are spectral of pressure and electric fields, respectively. (b) and (d) are time responses of pressure and electric fields, respectively.
(a)(b)(c)(d)In the spectral representation, the response allows the identification of the pesudo-Rayleigh wave (pR.1), which has an about cut-off frequency of 7 kHz. The frequency of the excited source pulse is enough to reveal this wave mode. The excitation is stronger at low frequency for Stoneley wave without cut-off frequency. Spectral representation of electric field (Figure1(c)) Ez shows that the borehole seismoelectric conversion dominates at low frequency. Figures 1(c) and 1(d) shows that P wave response may be identified easily in both spectral representation and waveforms of electric field. In this case the seismoelectric conversion efficiency of P wave is biggest, Stoneley wave is second.Figure2 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for fast formation generated by an acoustic dipole source. In the spectral representations of borehole pressure (Figure 2(a)), the two flexual modes obviously appear in the frequency range 0–20 kHz; however, in the spectral representations of magnetic field, the lowest mode is dominating. From Figures 2(b) and 2(d), the waveforms of borehole pressure with 3 kHz center frequency dominates by flexural guided waves; however, that of magnetic field has two groups with different velocity. Again, from Figure 2(d) we see a relatively weak but visible A-A group; it arrives simultaneously on all receivers. And the following groups having apparent velocities of acoustic flexural waves velocities of borehole acoustic waves.Spectral representation and time response of borehole pressure and magnetic fields for fast formation generated by dipole source. (a) and (c) are spectral of pressure and magnetic field respectively. (b) and (d) are time responses of pressure and magnetic field respectively.
(a)(b)(c)(d)Figure3 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr  for fast formation generated by an acoustic quadrupole source. In the spectral representations of borehole pressure (Figure 3(a)), the two screw modes obviously appear in the range of 0–20 kHz; however, in the spectral representations of magnetic field, the lowest mode is dominating. From Figures 3(b) and 3(d), the waveforms of acoustic waves quadrupole sources with the center frequency 5 kHz mainly dominated by screw waves and that of magnetic field consists of groups having apparent velocities of acoustic screw waves. And we can find almost no EM waves in magnetic field full waveforms.Spectral representation and time response of borehole pressure and magnetic fields for fast formation generated by quadrupole source. (a) and (c) are spectral of pressure and magnetic fields respectively. (b) and (d) are time responses of pressure and magnetic fields respectively.
(a)(b)(c)(d)


## 3.2. Slow Formation

Figure4 shows spectral response and synthetic waveforms for acoustic waves and electric field for slow formation generated by an acoustic monopole source. From Figures 4(b) and 4(d), we can find there are three wave groups in electric field Ez different from two in acoustic waveforms. A clear seismoelectric response A-A group for the electric field Ez is relatively weak but visible in full waveforms, which were recorded at all receivers simultaneously. And the following groups for the electric field having apparent velocities are comparable to those of borehole acoustic waves.Spectral representation and time response of borehole pressure and electric fields for slow formation generated by monopole source. (a) and (c) are spectral of pressure and electric fields, respectively. (b) and (d) are time responses of pressure and electric fields, respectively.
(a)(b)(c)(d)In the spectral representation, the pesudo-Rayleigh wave disappears and the response allows the identification of P wave and Stoneley wave. The excitation is stronger at low frequency for Stoneley wave for slow formation same as fast formation. Spectral representation of electric field (Figure4(c)) Ez shows that the borehole seismoelectric conversion dominates at low frequency. Figures 4(c) and 4(d) shows that P wave response may be identified easily in both spectral representation and waveforms of electric field. Different from fast formation case, the seismoelectric conversion efficiency of Stoneley wave is biggest, P wave is second.Figure5 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for slow formation generated by an acoustic dipole source.  In the spectral representations of acoustic (Figure 5 (a)), one flexural mode appears in the frequency range 0–20 kHz, and in the spectral representations of magnetic field (Figure 5 (b)), the flexural mode is dominating. From Figures 5(b) and 5(d), the waveforms of borehole pressure with 1 kHz center frequency dominates by flexural guided waves; however, that of magnetic field has two groups with different velocity. Again, from Figure 5(d) we see a relatively weak but visible A-A group; it arrives simultaneously on all receivers. And the following group has apparent velocities of acoustic flexural waves velocities of borehole acoustic waves. Again, it shows the borehole seismoelectric conversion efficiency is strong at low frequency range.Spectral representation and time response of borehole pressure and magnetic fields for slow formation generated by dipole source. (a) and (c) are spectral of pressure and magnetic fields, respectively. (b) and (d) are time responses of pressure and magnetic fields, respectively.
(a)(b)(c)(d)Figure6 shows spectral response and synthetic waveforms of acoustic waves and magnetic field Hr for slow formation generated by an acoustic quadrupole source. In the spectral representations of acoustic (Figure 6(a)), a screw mode appears in the range of 0–10 kHz; however, in the spectral representations of magnetic field, this mode is more evident. From Figures 6(b) and 6(d), the waveforms of acoustic waves excited by quadrupole source with the center frequency 3 kHz mainly dominated by screw waves and that of magnetic field consists of groups having apparent velocities of acoustic screw waves. We almost cannot find the EM wave in the full waveforms of magnetic field Hr.Spectral representation and time response of borehole pressure and magnetic fields for slow formation generated by quadrupole source. (a) and (c) are spectral of pressure and magnetic fields, respectively. (b) and (d) are time responses of pressure and magnetic fields, respectively.
(a)(b)(c)(d)


## 4. Conclusions and Discussions

In this paper, multipole seismoelectric logging waveforms are simulated in a fluid-filled borehole embedded in permeable formations based on the entire Pride seismoelectric theory. We consider the nonaxisymmetric seismoelectric field excited by acoustic multipole sources. Expressions for the multipole seismoelectric field quantities of interest were derived from the potential formulations. The full waveforms of acoustic waves and electric and magnetic fields in the time domain propagation in borehole are simulated by using a discrete wave number integration. The results show that electric and magnetic field detectors in a borehole can detect both the EM waves and the seismoelectric field. One can measure three components of electric field and magnetic field when higher order multipole sources (n≥1) are excited. The results shows that one can measure tens of millivolt electric fields as long as one selects an acoustic source with pressure amplitude being 1 MPa measuring 0.1 m from the source center. However, it is the weakness of electrical and magnetic signals in the presence of the big background noise that is one main obstacle to application of seismoelectric tools. Therefore, enhancing the signal-to-noise ratio is the key for seismoelectric logging measurement. The amplitudes of electromagnetic field should be related to formation parameters through coupling coefficients L(ω); the seismoelectric conversion efficiencies of the Stoneley wave, the flexural guided wave and the screw guided waves are slightly different due to borehole effect in nature. The borehole seismoelectric conversion mainly depends on the relative motion between a solid and a fluid, which is generated by different acoustic modes. The relationship between borehole seismoelectric conversion efficiencies and the formation parameters needs further research in detail.



---

*Source: 107827-2011-03-27.xml*