共查询到20条相似文献，搜索用时 765 毫秒
1.
Jing Ba José M. Carcione Hong Cao Fengchang Yao Qizhen Du 《Geophysical Prospecting》2013,61(3):599-612
We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data. 相似文献
2.
Philip Tillotson Angus Ian Best Jeremy Sothcott Clive McCann Wang Shangxu Xiang‐Yang Li 《Geophysical Prospecting》2011,59(1):111-119
P‐ and S‐wave velocity and attenuation coefficients (accurate to ±0.3% and ±0.2 dB/cm, respectively) were measured in synthetic porous rocks with aligned, penny‐shaped fractures using the laboratory ultrasonic pulse‐echo method. Shear‐wave splitting was observed by rotating the S‐wave transducer and noting the maximum and minimum velocities relative to the fracture direction. A block of synthetic porous rock of fracture density 0.0201 ± 0.0068 and fracture size 3.6 ± 0.38 mm (measured from image analysis of X‐ray CT scans) was sub‐sampled into three 20–30 mm long, 50 mm diameter core plugs oriented at 0°, 45° and 90° to the fracture normal (transversely isotropic symmetry axis). Full waveform data were collected over the frequency range 500–1000 kHz for both water and glycerin saturated cores to observe the effect of pore fluid viscosity at 1 cP and 100 cP, respectively. The shear‐wave splitting observed in the 90° core was 2.15 ± 0.02% for water saturated and 2.39 ± 0.02% for glycerin saturated, in agreement with the theory that suggests that the percentage splitting should be 100 times the fracture density and independent of the saturating fluid. In the 45° core, by contrast, splitting was 0.00 ± 0.02% for water saturation and ?0.77 ± 0.02% for glycerin saturation. This dependence on fracture orientation and pore fluid viscosity is consistent with the poro‐visco‐elastic theory for aligned, meso‐scale fractures in porous rocks. The results suggest the possible use of shear‐ or converted‐wave data to discriminate between fluids on the basis of viscosity variations. 相似文献
3.
Hua‐You Chai Tian‐Bin Li Kok‐Kwang Phoon Elton J. Chen Dian‐Ji Zhang 《Geophysical Prospecting》2017,65(4):992-1003
In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms. 相似文献
4.
Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients. 相似文献
5.
Lars O. Løseth 《Geophysical Prospecting》2011,59(1):145-160
In marine controlled‐source electromagnetic (CSEM) surveys the subsurface is explored by emitting low‐frequency signals from an electric dipole source close to the sea‐bed. The main goal is often to detect and describe possible thin resistive layers beneath the sea‐bed. To gain insight into how CSEM signals propagate, it is informative to study a stratified model. The electromagnetic field is then given in terms of integrals over TE‐ and TM‐polarized plane‐wave constituents. An asymptotic evaluation of the field integrals for large propagation distances results in explicit spatial expressions for the field components and the derived expressions can be used to analyse how the CSEM signals propagate. There are two major signal pathways in a standard CSEM model. One of these pathways is via the thin resistive layer and the resulting response is accounted for by a pole in the reflection response for the TM mode. The signal is propagating nearly vertically down to the resistor from the source, then guided while attenuated along the reservoir, before propagating nearly vertically up to the receiver. The response is slightly altered by the sea‐bed interface and further modified in shallow water due to multiple reflections between the sea‐surface and sea‐bed at both the source and receiver sides. The other major signal pathway is via the resistive air half‐space, the so‐called airwave. The airwave is generated by the TE mode and interacts with the subsurface via vertically propagating signals reflected between the sea‐surface and subsurface at both the source and receiver sides. 相似文献
6.
Dongdong Wang Yongxin Gao Cheng Yao Baozhen Wang Mengqiang Wang 《Geophysical Prospecting》2020,68(6):1958-1979
We investigate the seismoelectric/electroseismic wavefields excited by a point source in an air/seawater/three-layered porous medium configuration containing a hydrocarbon layer. The results show that if an explosive source for excitation is used, receivers at seafloor can record the coseismic electromagnetic fields accompanying the P, S, fluid acoustic waves and the interface responses converted from the acoustic waves at seafloor interface and from the seismic waves at the interfaces beneath the seafloor. Employing a vertical electric dipole source shows that, with the exception of the interface responses converted from electromagnetic waves at seafloor, the interface responses converted from transmitted electromagnetic waves at the interfaces beneath the seafloor can also be identified. Given that the strength of the explosive source is within excitation capability of industry air guns, the generated interface responses from the hydrocarbon layer can be detected by current electromagnetic sensors considering the low ambient noise at the seafloor. Our results demonstrate the feasibility of the seismoelectric method applied to marine hydrocarbon exploration. Electroseismic modelling results suggest that it is not practical to employ this method to prospect marine hydrocarbon layer due to the weak interface response signal, unless a much larger current is injected into seafloor. 相似文献
7.
Ho Seuk Bae Changsoo Shin Young Ho Cha Yunseok Choi Dong‐Joo Min 《Geophysical Prospecting》2010,58(6):997-1010
Although waveform inversion has been intensively studied in an effort to properly delineate the Earth's structures since the early 1980s, most of the time‐ and frequency‐domain waveform inversion algorithms still have critical limitations in their applications to field data. This may be attributed to the highly non‐linear objective function and the unreliable low‐frequency components. To overcome the weaknesses of conventional waveform inversion algorithms, the acoustic Laplace‐domain waveform inversion has been proposed. The Laplace‐domain waveform inversion has been known to provide a long‐wavelength velocity model even for field data, which may be because it employs the zero‐frequency component of the damped wavefield and a well‐behaved logarithmic objective function. However, its applications have been confined to 2D acoustic media. We extend the Laplace‐domain waveform inversion algorithm to a 2D acoustic‐elastic coupled medium, which is encountered in marine exploration environments. In 2D acoustic‐elastic coupled media, the Laplace‐domain pressures behave differently from those of 2D acoustic media, although the overall features are similar to each other. The main differences are that the pressure wavefields for acoustic‐elastic coupled media show negative values even for simple geological structures unlike in acoustic media, when the Laplace damping constant is small and the water depth is shallow. The negative values may result from more complicated wave propagation in elastic media and at fluid‐solid interfaces. Our Laplace‐domain waveform inversion algorithm is also based on the finite‐element method and logarithmic wavefields. To compute gradient direction, we apply the back‐propagation technique. Under the assumption that density is fixed, P‐ and S‐wave velocity models are inverted from the pressure data. We applied our inversion algorithm to the SEG/EAGE salt model and the numerical results showed that the Laplace‐domain waveform inversion successfully recovers the long‐wavelength structures of the P‐ and S‐wave velocity models from the noise‐free data. The models inverted by the Laplace‐domain waveform inversion were able to be successfully used as initial models in the subsequent frequency‐domain waveform inversion, which is performed to describe the short‐wavelength structures of the true models. 相似文献
8.
Xiaoyu Zhang Qunshan Zhang George A. McMechan Gladys Gonzalez 《Geophysical Prospecting》2017,65(3):669-686
Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region. 相似文献
9.
This study investigates the control of jacket‐type offshore platforms. The deck displacement of jacket‐type offshore platforms can be controlled using both passive and active control mechanisms. Among the passive control mechanisms, a tuned mass damper concept is studied in this paper. Active control mechanisms considered here include the active mass damper, the active tendon mechanism and the propeller thruster. An optimal frequency domain approach to active control of wave‐excited platforms is used in which the H2 norm of the transfer function from the external disturbance to the regulated output is minimized. In this study, the hydrodynamic drag force is evaluated using the JONSWAP wave spectrum. Unlike conventional linearization approaches, the influence of non‐linearity in the drag force is retained in this scheme by expressing the non‐linear force components in terms of higher‐order convolutions of the water‐particle velocities. To demonstrate the effectiveness of this scheme, the platform performance with and without control devices under different sea states is evaluated. It is demonstrated that the control devices are useful in reducing the displacement response of jacket‐type offshore platforms, especially when the wave forces are concentrated at frequencies close to the natural frequencies of the platform. This becomes especially significant in deep waters because the natural frequencies of jacket‐type platforms fall closer to the dominant wave frequencies in deep waters. Adding control devices to deep water platforms will ensure a reduction both in the global response of the platform and in localized effects, such as the fatigue of welded joints. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
10.
Hysteresis represents a loop in a rating curve and is a phenomenon which closely resembles that occurring in stress–strain curves used for studying the elastic properties of solid substances in engineering mechanics. Earlier hysteresis‐based studies used for defining floodwave propagation in open channels have qualitatively shown that hysteresis is an index of energy loss during floodwave propagation. Using the concept of elasticity, this paper introduces a new term called flow strain (defined as the ratio of change in discharge to the initial discharge) for investigating hysteresis. The usefulness of this new term is evaluated with use of four dam‐break studies. The study reveals that:
- 1 flow strain is a function of three wave speeds, Seddon speed, Lagrange speed, and elastic speed;
- 2 a single linear reservoir concept frequently used in flood routing is a specific variant of the Seddon speed formula;
- 3 the non‐linear storage–discharge relationship, widely used in overland flow modelling, is a variant of the kinematic wave representation;
- 4 the discharge ordinates on the recession part of a hydrograph follow a simple first‐order autoregressive form;
- 5 the hysteresis, phase difference and logarithmic decrement all define attenuation and are indices of energy loss during floodwave propagation.
11.
The phenomenon of acoustic waves inducing electric fields in porous media is called the seismoelectric effect. Earlier investigators proposed the usage of seismoelectric effect for well logging. Soil texture has a strong influence on the coupled wave fields during shallow surface explorations. In this article, we study the borehole pure shear‐horizontal wave and the coupling transverse‐electric field (acoustic–electrical coupling wave fields) in the partially saturated soil. Combined with related theories, we expand the formation parameters to partially saturated forms and discuss the influence of soil texture conditions on the seismoelectric wave fields. The results show that the elastic and electrical properties of porous media are sensitive to water saturation. The compositions of the acoustic and electric fields for different soil textures do not change, but the waveforms differ. We also use the secant integral method to simulate the interface‐converted electromagnetic waves. The results show that interface response strength is greatly influenced by soil texture. In addition, considering the sensitivity of the inducing electric field to fluid salinity, we also simulate the time‐domain waveforms of electric field for different pore fluid salinity levels. The results show that as the salinity increases, the electric field amplitude decreases monotonically. The above conclusions have certain significance for the application of borehole shear wave and its coupled electric fields for resource exploration, saturation assessment and groundwater pollution monitoring. 相似文献
12.
Influence of water pressure dynamics and fluid flow on the streaming‐potential response for unsaturated conditions 
下载免费PDF全文
下载免费PDF全文 V. Allègre L. Jouniaux F. Lehmann P. Sailhac R. Toussaint 《Geophysical Prospecting》2015,63(3):694-712
Streaming‐potentials are produced by electrokinetic effects in relation to fluid flow and are used for geophysical prospecting. The aim of this study is to model streaming potential measurements for unsaturated conditions using an empirical approach. A conceptual model is applied to streaming potential measurements obtained from two drainage experiments in sand. The streaming potential data presented here show a non‐monotonous behaviour with increasing water saturation, following a pattern that cannot be predicted by existing models. A model involving quasi‐static and dynamic components is proposed to reproduce the streaming potential measurements. The dynamic component is based on the first time derivative of the driving pore pressure. The influence of this component is investigated with respect to fluid velocity, which is very different between the two experiments. The results demonstrate that the dynamic component is predominant at the onset of drainage in experiments with the slowest water flow. On the other hand, its influence appears to vanish with increasing drainage velocity. Our results suggest that fluid flow and water distribution at the pore scale have an important influence on the streaming potential response for unsaturated conditions. We propose to explain this specific streaming potential response in terms of the behaviour of both rock/water interface and water/air interfaces created during desaturation processes. The water/air interfaces are negatively charged, as also observed in the case of water/rock interfaces. Both the surface area and the flow velocity across these interfaces are thought to contribute to the non‐monotonous behaviour of the streaming potential coefficient as well as the variations in its amplitude. The non‐monotonous behaviour of air/water interfaces created during the flow was highlighted as it was measured and modelled by studies published in the literature. The streaming potential coefficient can increase to about 10 to 40 when water saturation decreases. Such an increase is possible if the amount of water/air interfaces is increased in sufficient amount, which can be the case. 相似文献
13.
Recently, mode converted shear waves (C‐waves) have been shown to enable overpressure prediction in media where primary wave acquisition is inhibited by gas and fluid effects – C‐wave moveout is analysed and a long standing relationship between differential stress and primary‐wave (P‐wave) velocity is modified and employed. Though pore‐pressure prediction based on C‐waves is supported by empirical evidence from laboratory and field experiments, a theoretical justification has yet to be developed. In this research note, we provide a supporting algebra for the original relationship between pore pressure and C‐wave velocity. 相似文献
14.
The nature of intra‐annual variability in the non‐advective heat fluxes affecting streams and rivers in Devon, UK was investigated through detailed monitoring of study reaches in an upland moorland catchment, below a regulating reservoir, and flowing through deciduous woodland and coniferous forest during the period May 1995 to April 1996. A clear pattern of seasonal variation was evident, whereby net radiation provided a heat source during the summer but a heat sink in the winter, as incoming short‐wave radiation declined and outgoing long‐wave radiation increased. Sensible transfer added heat to the study reaches in the summer but removed it during the winter, and bed conduction acted as a heat sink in the summer period but as a heat source in the winter months. Friction and evaporation added and removed heat, respectively, from the study reaches throughout the year, but the magnitude of these fluxes reflected seasonal variations in discharge and in wind speed. Water temperature generally followed the net non‐advective heat energy budget, which was positive in summer but negative in winter. Although a general pattern of seasonal variability in the non‐advective heat energy budget was evident, detailed differences in the nature and extent of intra‐annual variability were apparent between the study reaches and particularly between forested and non‐forested sites. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
15.
Heike Sütterlin Rainer Trittler Sebastian Bojanowski Ernst A. Stadlbauer Klaus Kümmerer 《洁净——土壤、空气、水》2007,35(1):81-87
Thermocatalytic low temperature conversion (LTC) is a new method for the disposal of sewage sludge. Using this method, sludge is converted into a residual solid (coal) along with reaction water, oil, and non‐condensable gases. The oil can be used as an energy source and the coal as a substitute for charcoal. To this end, it is important to determine whether there are any easily available contaminants present in the coal generated by the process. Contaminants that can be strongly sorbed by sewage sludge solids are, e. g., pharmaceuticals and disinfectants. As an example the fate of the persistent and strong adsorbing disinfectant benzalkonium chloride (BAC) has been investigated within the LTC process. The sewage sludge was spiked with BAC and then subjected to the LTC process. The resulting coal was extracted and analyzed using LC‐LC/ESI‐MS/MS (ion trap). BAC could not be detected in the LTC coal, although it could be extracted from the spiked sludge before the LTC treatment. It can thus be concluded that the investigated compound is not easily available in the coal, and hence that its use does not present a risk. 相似文献
16.
Guangquan Li 《Geophysical Prospecting》2020,68(3):910-917
Biot theory was based on two ideas: the coupling factor to quantify the kinetic energy of fluid and Darcy permeability to quantify the dissipation function. As Biot theory did not well predict attenuation of ultrasonic S wave, we modify the theory to better characterize the S wave attenuation. The range of the coupling factor is at first estimated in view of fluid mechanics. Application of the original theory to water-saturated Boise sandstone and brine-saturated Berea sandstone shows that the model prediction significantly underestimates the S wave attenuation ultrasonically measured. For this reason, we replace Darcy permeability with variable permeability to improve the fluid momentum equation. The new model yields predictions of phase velocity and the quality factor both close to the ultrasonic measurements. The reason why the improved model is superior to Biot theory is that variable permeability is based on the Stokes boundary layer at the fluid–solid interface, thus accurately quantifying the viscous stress between the two phases. Finally, the length scale of the viscous stress is calculated in the mesoscopic sense. 相似文献
17.
Analysis of time‐lapse travel‐time and amplitude changes to assess reservoir compartmentalization 下载免费PDF全文
下载免费PDF全文 Fluid depletion within a compacting reservoir can lead to significant stress and strain changes and potentially severe geomechanical issues, both inside and outside the reservoir. We extend previous research of time‐lapse seismic interpretation by incorporating synthetic near‐offset and full‐offset common‐midpoint reflection data using anisotropic ray tracing to investigate uncertainties in time‐lapse seismic observations. The time‐lapse seismic simulations use dynamic elasticity models built from hydro‐geomechanical simulation output and a stress‐dependent rock physics model. The reservoir model is a conceptual two‐fault graben reservoir, where we allow the fault fluid‐flow transmissibility to vary from high to low to simulate non‐compartmentalized and compartmentalized reservoirs, respectively. The results indicate time‐lapse seismic amplitude changes and travel‐time shifts can be used to qualitatively identify reservoir compartmentalization. Due to the high repeatability and good quality of the time‐lapse synthetic dataset, the estimated travel‐time shifts and amplitude changes for near‐offset data match the true model subsurface changes with minimal errors. A 1D velocity–strain relation was used to estimate the vertical velocity change for the reservoir bottom interface by applying zero‐offset time shifts from both the near‐offset and full‐offset measurements. For near‐offset data, the estimated P‐wave velocity changes were within 10% of the true value. However, for full‐offset data, time‐lapse attributes are quantitatively reliable using standard time‐lapse seismic methods when an updated velocity model is used rather than the baseline model. 相似文献
18.
Colin M. Sayers 《Geophysical Prospecting》2002,50(4):393-401
Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. A shear wave propagating in an azimuthally anisotropic medium splits into two components with different polarizations if the source polarization is not aligned with the principal axes of the medium. For vertical propagation of shear waves in a horizontally layered medium containing vertical fractures, the shear‐wave splitting depends on the shear compliance of the fractures, but is independent of their normal compliance. If the fractures are not perfectly vertical, the shear‐wave splitting also depends on the normal compliance of the fractures. The normal compliance of a fluid‐filled fracture decreases with increasing fluid bulk modulus. For dipping fractures, this results in a decrease in shear‐wave splitting and an increase in shear‐wave velocity with increasing fluid bulk modulus. The sensitivity of the shear‐wave splitting to fluid bulk modulus depends on the interconnectivity of the fracture network, the permeability of the background medium and on whether the fracture is fully or partially saturated. 相似文献
19.
Ho Seuk Bae Sukjoon Pyun Wookeen Chung Seung‐Goo Kang Changsoo Shin 《Geophysical Prospecting》2012,60(3):413-432
We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open. 相似文献
20.
The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly. 相似文献