Biot theory research has been extended to the multi-scale heterogeneity in actual rocks. Focused on laboratory frequency bandwidth studies, we discuss the relationships between double-porosity and BISQ wave equations,...Biot theory research has been extended to the multi-scale heterogeneity in actual rocks. Focused on laboratory frequency bandwidth studies, we discuss the relationships between double-porosity and BISQ wave equations, analytically derive the degeneration method for double-porosity＇s return to BISQ, and give three necessary conditions which the degeneration must satisfy. By introducing dynamic permeability and tortuosity theory, a full set of dynamic double-porosity wave equations are derived. A narrow band approximation is made to simplify the numerical simulation for dynamic double-porosity wavefields. Finally, the pseudo-spectral method is used for wave simulation within the laboratory frequency band （50 kHz）. Numerical results have proved the feasibility for dynamic double-porosity＇s description of squirt flow and the validity of the quasi-static approximation method.展开更多
Phenomena of reflection and refraction of plane harmonic waves at a plane interface between an elastic solid and doubleporosity dual-permeability material are investigated. The elastic solid behaves non-dissipatively,...Phenomena of reflection and refraction of plane harmonic waves at a plane interface between an elastic solid and doubleporosity dual-permeability material are investigated. The elastic solid behaves non-dissipatively, while double-porosity dual-permeability materials behave dissipatively to wave propagation due to the presence of viscosity in pore fluids. All the waves(i.e., incident and reflected) in an elastic medium are considered as homogeneous(i.e., having the same directions of propagation and attenuation), while all the refracted waves in double-porosity dual-permeability materials are inhomogeneous(i.e., having different directions of propagation and attenuation). The coefficients of reflection and refraction for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected and refracted waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy among various reflected and refracted waves. The effect of incident direction on the partition of the incident energy is analyzed with a change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure.It has been confirmed from numerical interpretation that during the reflection/refraction process, conservation of incident energy is obtained at each angle of incidence.展开更多
The present work deals with the reflection of plane seismic waves at the stress-free plane surface of double-porosity dualpermeability material. The incidence of two main waves（i.e., P1 and SV） is considered. As a r...The present work deals with the reflection of plane seismic waves at the stress-free plane surface of double-porosity dualpermeability material. The incidence of two main waves（i.e., P1 and SV） is considered. As a result of the incident waves,four reflected（three longitudinal and one shear） waves are found in the medium. The expressions of reflection coefficients for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy for fully closed as well as perfectly open pores. Effect of incident direction on the partition of the incident energy is analyzed with the change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure. It has been confirmed from the numerical interpretation that during the reflection process, conservation of incident energy is obtained at each angle of incidence.展开更多
Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in ...Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in the double-porosity medium,two porous skeletons are usually assumed,namely,host and inclusions.Of them,the volume ratio of inclusion skeletons is low.All previous studies have ignored the consideration of local fluid flow velocity field in inclusions,and therefore they can not completely describe the physical process of local flow oscillation and should not be applied to the situation where the fluid kinetic energy in inclusions cannot be neglected.In this paper,we analyze the local fluid flow velocity fields inside and outside the inclusion,rewrite the kinetic energy function and dissipation function based on the double-porosity medium model containing spherical inclusions,and derive the reformulated Biot-Rayleigh(BR)equations of elastic wave propagation based on Hamilton’s principle.We present simulation examples with different rock and fluid types.Comparisons between BR equations and reformulated BR equations show that there are significant differences in wave response characteristics.Finally,we compare the reformulated BR equations with the previous theories and experimental data,and the results show that the theoretical results of this paper are correct and effective.展开更多
Reservoir boundary shape has a great influence on the transient pressure response of oil wells located in arbitrarily shaped reservoirs. Conventional analytical methods can only be used to calculate transient pressure...Reservoir boundary shape has a great influence on the transient pressure response of oil wells located in arbitrarily shaped reservoirs. Conventional analytical methods can only be used to calculate transient pressure response in regularly shaped reservoirs. Under the assumption that permeability varies exponentially with pressure drop, a mathematical model for well test interpretation of arbitrarily shaped deformable reservoirs was established. By using the regular perturbation method and the boundary element method, the model could be solved. The pressure behavior of wells with wellbore storage and skin effects was obtained by using the Duhamel principle. The type curves were plotted and analyzed by considering the effects of permeability modulus, arbitrary shape and impermeable region.展开更多
基金supported by the National Basic Research Program of China (973 Program, 2007CB209505)the International Cooperative Project of the Ministry of Science and Technology of China (2006DFB62030)
文摘Biot theory research has been extended to the multi-scale heterogeneity in actual rocks. Focused on laboratory frequency bandwidth studies, we discuss the relationships between double-porosity and BISQ wave equations, analytically derive the degeneration method for double-porosity＇s return to BISQ, and give three necessary conditions which the degeneration must satisfy. By introducing dynamic permeability and tortuosity theory, a full set of dynamic double-porosity wave equations are derived. A narrow band approximation is made to simplify the numerical simulation for dynamic double-porosity wavefields. Finally, the pseudo-spectral method is used for wave simulation within the laboratory frequency band （50 kHz）. Numerical results have proved the feasibility for dynamic double-porosity＇s description of squirt flow and the validity of the quasi-static approximation method.
文摘Phenomena of reflection and refraction of plane harmonic waves at a plane interface between an elastic solid and doubleporosity dual-permeability material are investigated. The elastic solid behaves non-dissipatively, while double-porosity dual-permeability materials behave dissipatively to wave propagation due to the presence of viscosity in pore fluids. All the waves(i.e., incident and reflected) in an elastic medium are considered as homogeneous(i.e., having the same directions of propagation and attenuation), while all the refracted waves in double-porosity dual-permeability materials are inhomogeneous(i.e., having different directions of propagation and attenuation). The coefficients of reflection and refraction for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected and refracted waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy among various reflected and refracted waves. The effect of incident direction on the partition of the incident energy is analyzed with a change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure.It has been confirmed from numerical interpretation that during the reflection/refraction process, conservation of incident energy is obtained at each angle of incidence.
文摘The present work deals with the reflection of plane seismic waves at the stress-free plane surface of double-porosity dualpermeability material. The incidence of two main waves（i.e., P1 and SV） is considered. As a result of the incident waves,four reflected（three longitudinal and one shear） waves are found in the medium. The expressions of reflection coefficients for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy for fully closed as well as perfectly open pores. Effect of incident direction on the partition of the incident energy is analyzed with the change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure. It has been confirmed from the numerical interpretation that during the reflection process, conservation of incident energy is obtained at each angle of incidence.
基金supported by the National Natural Science Foundation of China(Grant No.41104066)RIPED Youth Innovation Foundation(Grant No.2010-A-26-01)+1 种基金the National Basic Research Program of China(Grant No.2014CB239006)the Open fund of SINOPEC Key Laboratory of Geophysics(Grant No.WTYJY-WX2013-04-18)
文摘Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in the double-porosity medium,two porous skeletons are usually assumed,namely,host and inclusions.Of them,the volume ratio of inclusion skeletons is low.All previous studies have ignored the consideration of local fluid flow velocity field in inclusions,and therefore they can not completely describe the physical process of local flow oscillation and should not be applied to the situation where the fluid kinetic energy in inclusions cannot be neglected.In this paper,we analyze the local fluid flow velocity fields inside and outside the inclusion,rewrite the kinetic energy function and dissipation function based on the double-porosity medium model containing spherical inclusions,and derive the reformulated Biot-Rayleigh(BR)equations of elastic wave propagation based on Hamilton’s principle.We present simulation examples with different rock and fluid types.Comparisons between BR equations and reformulated BR equations show that there are significant differences in wave response characteristics.Finally,we compare the reformulated BR equations with the previous theories and experimental data,and the results show that the theoretical results of this paper are correct and effective.
文摘Reservoir boundary shape has a great influence on the transient pressure response of oil wells located in arbitrarily shaped reservoirs. Conventional analytical methods can only be used to calculate transient pressure response in regularly shaped reservoirs. Under the assumption that permeability varies exponentially with pressure drop, a mathematical model for well test interpretation of arbitrarily shaped deformable reservoirs was established. By using the regular perturbation method and the boundary element method, the model could be solved. The pressure behavior of wells with wellbore storage and skin effects was obtained by using the Duhamel principle. The type curves were plotted and analyzed by considering the effects of permeability modulus, arbitrary shape and impermeable region.
