An analysis of wave dissipation at the Hendijan mud coast, the Persian Gulf

Effects of non-rigid muddy bed on the wave climate at the Hendijan coast along the northwestern part of the Persian Gulf have been examined through field measurements and numerical wave transformation modeling. The field survey included measurements of wave characteristics at an offshore and a nearshore station, and mud sampling to obtain the thickness of the fluid mud layer and its rheological properties. Comparisons of wave spectra at the two stations show energy dissipation along the wave trajectory with higher dissipation in the wave period band around 6?s, because depending on the site a given frequency band tends to be more effective in wave–mud interaction. Dissipation induced by the non-rigid bed is introduced into the REF/DIF wave transformation model through the application of viscoelastic constitutive equations for fluid mud. Numerical outputs of the nearshore wave height, for which the viscoelastic parameters included in the model were obtained independently from oscillatory frequency-sweep tests, are found to be comparable with measured values at the nearshore station. This implies that the model is useful for estimating the design wave conditions in the study area.     

A comparative study on the rheology and wave dissipation of kaolinite and natural Hendijan Coast mud, the Persian Gulf

The objective of this paper is to investigate the rheological behavior of kaolinite and Hendijan mud, located at the northwest part of the Persian Gulf, and the dissipative role of this muddy bed on surface water waves. A series of laboratory rheological tests was conducted to investigate the rheological response of mud to rotary and cyclic shear rates. While a viscoplastic Bingham model can successfully be applied for continuous controlled shear-stress tests, the rheology of fluid mud displays complex viscoelastic behavior in time-periodic motion. The comparisons of the behavior of natural Hendijan mud with commercial kaolinite show rheological similarities. A large number of laboratory wave-flume experiments were carried out with a focus on the dissipative role of the fluid mud. Assuming four rheological models of viscous, Kelvin-Voigt viscoelastic, Bingham viscoplastic, and viscoelastic-plastic for fluid mud layer, a numerical multi-layered model was applied to analyze the effects of different parameters of surface wave and muddy bed on wave attenuation. The predicted results based on different rheological models generally agree with the obtained wave-flume data implying that the adopted rheological model does not play an important role in the accuracy of prediction.     

A two-dimensional horizontal wave propagation and mud mass transport model

The present study offers a two-dimensional horizontal wave propagation and morphodynamic model for muddy coasts. The model can be applied on a general three-dimensional bathymetry of a soft muddy coast to calculate wave damping, fluid mud mass transport and resulting bathymetry change under wave actions. The wave propagation model is based on time-dependent mild slope equations including the wave energy dissipation due to the wave-mud interaction of bottom mud layers as well as the combined effects of the wave refraction, diffraction and breaking. The constitutive equations of the visco-elastic–plastic model are adopted for the rheological behavior of fluid mud. The mass transport velocity within the fluid mud layer is calculated combining the Stokes’ drift, the mean Eulerian velocity and the gravity-driven mud flow. The results of the numerical model are compared against a series of conducted wave basin experiments, wave flume experiments and field observations. Comparisons between the computed results with both the field and laboratory data reveal the capability of the proposed model to predict the wave transformation and mud mass transport.     

EXPERIMENTAL STUDY ON WAVE ENERGY DISSIPATION AND COHESIVE SEDIMENT TRANSPORT IN SILT COAST

EmmmL STonY ON w^rs EareGY DlsSunnox AunCOIIES1VE SEDonT TRANSPORT 1N SirT COASTShaong HUl and Onyx WAI2Abstract:The interahon betwee the wave and fluld mud laye Plays an twrtan fOle in the develOPmen ofsilt coast. Sediment is essenhally transported in the form of rheological flow of mud laye unds thewave achon, and on the other hand, the fluid Inud layer dams the wave consfory. ms papestUdies the laws of wave energy dissiPation and mud bed deformhon, and the moveInen of…     

Mud erosion by waves: a laboratory study

The role of mud erosion under waves in governing cohesive sediment transport in estuarial and coastal waters is well known. A laboratory study was conducted in order to elucidate the mechanism by which soft muds erode under progressive waves in a flume. Two types of cohesive sediment were used, a commercial kaolinite and an estuarial mud. Beds were formed by pouring in a pre-prepared sediment-water slurry and allowing the deposit to consolidate for a period ranging from 2 to 14 days. A multi-layered hydrodynamic model, which considers the mud to be viscoelastic, has been developed and used to evaluate the bed shear stress at the oscillating mud-water interface. The viscoelastic property of the mud has been confirmed by rheological measurements, and model results on velocity, pressure and wave attenuation verified against flume data. Concentration profiles indicate a distinct evolutionary pattern resulting in a highly stratified suspension. Just above the bed, a thin layer of fluid mud is generated. Above this layer, the suspension concentration is significantly lower. This two-layered feature of the concentration profile is related to the oscillatory response of the mud and water layers, and the associated momentum exchange and mass diffusion characteristics. An expression relating the rate of erosion to the bed shear stress in excess of bed shear resistance has been developed. Generation of fluid mud during erosion is a significant feature of the role of waves over mud.     

A study of dissipation of wind-waves by mud at Cassino Beach,Brazil: Prediction and inversion

The impact of a non-rigid seafloor on the wave climate at Cassino Beach, Brazil, May–June 2005 is studied using field measurements and a numerical wave model. The measurements consist of wave data at four locations; rheology and mud thickness from grab samples; and an estimate of the horizontal distribution of mud based on echo-soundings. The dissipation of waves by a non-rigid bottom is represented in the wave model by treating the mud layer as a viscous fluid. Applied for 431 time periods, the model without this type of dissipation has a strong tendency to overpredict nearshore wave energy, except during a period of large storm waves. Two model variations which include this dissipation have a modest tendency to underpredict the nearshore wave energy. An inversion methodology is developed and applied to infer an alternate mud distribution which, when used with the wave model, yields the observed waveheights.     

Implementation of viscoelastic mud-induced energy attenuation in the third-generation wave model,SWAN

The interaction of waves with fluid mud can dissipate the wave energy significantly over few wavelengths. In this study, the third-generation wave model, SWAN, was advanced to include attenuation of wave energy due to interaction with a viscoelastic fluid mud layer. The performances of implemented viscoelastic models were verified against an analytical solution and viscous formulations for simple one-dimensional propagation cases. Stationary and non-stationary test cases in the Surinam coast and the Atchafalaya Shelf showed that the inclusion of the mud-wave interaction term in the third-generation wave model enhances the model performance in real applications. A high value of mud viscosity (of the order of 0.1 m²/s) was required in both field cases to remedy model overestimation at high frequency ranges of the wave spectrum. The use of frequency-dependent mud viscosity value improved the performance of model, especially in the frequency range of 0.2–0.35 Hz in the wave spectrum. In addition, the mud-wave interaction might affect the high frequency part of the spectrum, and this part of the wave spectrum is also affected by energy transfer from wind to waves, even for the fetch lengths of the order of 10 km. It is shown that exclusion of the wind input term in such cases might result in different values for parameters of mud layer when inverse modeling procedure was employed. Unlike viscous models for wave-mud interaction, the inverse modeling results to a set of mud parameters with the same performance when the viscoelastic model is used. It provides an opportunity to select realistic mud parameters which are in more agreement with in situ measurements.     

Analysis of attenuation and dispersion of propagating wave due to the coexistence of three fluid phases in the pore volume

In exploration geophysics, the efforts to extract subsurface information from wave characteristics exceedingly depend on the construction of suitable rock physics model. Analysis of different rock physics models reveals that the strength and magnitude of attenuation and dispersion of propagating wave exceedingly depend on wave-induced fluid flow at multiple scales. In current work, a comprehensive analysis of wave attenuation and velocity dispersion is carried out at broad frequency range. Our methodology is based on Biot's poroelastic relations, by which variations in wave characteristics associated with wave-induced fluid flow due to the coexistence of three fluid phases in the pore volume is estimated. In contrast to the results of previous research, our results indicate the occurrence of two-time pore pressure relaxation phenomenon at the interface between fluids of disparate nature, that is, different bulk modulus, viscosity and density. Also, the obtained results are compatible with numerical results for the same 1D model which are accounted using Biot's poroelastic and quasi-static equation in frequency domain. Moreover, the effects of change in saturation of three-phase fluids were also computed which is the key task for geophysicist. The outcomes of our research reveal that pore pressure relaxation phenomenon significantly depends on the saturation of distinct fluids and the order of saturating fluids. It is also concluded that the change in the saturation of three-phase fluid significantly influences the characteristics of the seismic wave. The analysis of obtained results indicates that our proposed approach is a useful tool for quantification, identification and discrimination of different fluid phases. Moreover, our proposed approach improves the accuracy to predict dispersive behaviour of propagating wave at sub-seismic and seismic frequencies.     

Fluidization and representative wave transformation on muddy beds

Irregular wave-induced mud fluidization and wave spectrum transformation on muddy profiles are studied through representative wave technique. The constitutive equations of visco-elastic model are adopted for the rheological behavior of fluid mud, while the behavior of stationary mud is assumed to be elastic. A set of representative waves are employed to investigate wave–mud interaction. The results are verified using real field data. Comparing the performance of common representative waves, it is concluded that the phenomena can be better predicted by root mean square wave.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明：本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.

     

岩石的粘弹性谐振Q模型

实验表明,在应变小于10^-6范围内,砂岩对地震波的吸收主要由孔隙流体的局部运动引起,而且Q值随频率的变化出现谐振现象。据此,作者认为,地壳上部地震波能量的耗损可以表示为以下两种成分的线性组合:1.由滑动摩擦、热弹性驰豫、位错运动等引起的能量消散,它们主要与岩石的固相成分有关,可近似地用恒Q模型描述;2.由孔隙流体运动引起的能量消散,可以用谐振公式描述。根据这一认识,利用描述因果关系的Kramers-Krnig关系式可推导出表示地壳上部岩石粘弹性（复弹性模量,相速度频散和衰减函数）的公式,它们综合地描述了由各种机制引起的波的频散和吸收,并在谐振Q值等于参考常数Q值时退化为目前常用的Futterman模型。作为这种谐振Q模型的应用,介绍了它用于Q值测量结果外推和频散一吸收研究以及粘弹性介质中反射地震道合成的结果。     

岩石的粘弹性谐振Q模型

实验表明,在应变小于10-6范围内,砂岩对地震波的吸收主要由孔隙流体的局部运动引起,而且Q值随频率的变化出现谐振现象。据此,作者认为,地壳上部地震波能量的耗损可以表示为以下两种成分的线性组合:1.由滑动摩擦、热弹性驰豫、位错运动等引起的能量消散,它们主要与岩石的固相成分有关,可近似地用恒Q模型描述;2.由孔隙流体运动引起的能量消散,可以用谐振公式描述。根据这一认识,利用描述因果关系的Kramers-Krnig关系式可推导出表示地壳上部岩石粘弹性（复弹性模量,相速度频散和衰减函数）的公式,它们综合地描述了由各种机制引起的波的频散和吸收,并在谐振Q值等于参考常数Q值时退化为目前常用的Futterman模型。作为这种谐振Q模型的应用,介绍了它用于Q值测量结果外推和频散一吸收研究以及粘弹性介质中反射地震道合成的结果。     

Viscoelastic representation of surface waves in patchy saturated poroelastic media

Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1–1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot’s differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.     

The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.     

Wave propagation analysis of porous rocks with the thermal activated relaxation mechanism

Wave attenuation and phase velocity dispersion in the temperature domain are more complicated than those in the frequency domain. To describe wave propagation properties in the temperature domain, a so-called thermal activation mechanism model is built on the experimental result that increasing the temperature or decreasing frequency could obtain similar results on the attenuation. A rheological model (the Zener model) is employed to describe viscoelastic attenuation in saturated porous rocks. The Arrhenius relation is introduced to describe the thermal activation mechanism. The wave propagation model with thermal effects in porous media is then obtained, and 1-D P-wave and S-wave propagation characteristics are analyzed in numeric process, respectively.Two attenuation mechanisms are found in this model, the Biot loss and the thermal activation relaxation. The thermal relaxation attenuation peak and the Biot attenuation peak are observed in both frequency and temperature spectra. These two peaks move towards each other when the temperature increases on frequency spectra. The thermal relaxation peak shifts towards higher frequencies while the Biot peak shifts towards lower frequencies. At some temperature, these two peaks will superpose. The combination of the thermal relaxation and the Biot loss leads to the complexity of wave velocity curves. Similar phenomena could be observed on temperature spectra. The thermal relaxation features may relate to a so-called “local heat transfer” mechanism. These two peaks in the temperature domain have been observed in the experiments by other investigators. The characteristics of velocity and attenuation are more remarkable for high porosity rock samples. The model is helpful for the understanding of wave propagation in the temperature domain.     

Angular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous Media

Wave-induced fluid flow generates a dominant attenuation mechanism in porous media. It consists of energy loss due to P-wave conversion to Biot (diffusive) modes at mesoscopic-scale inhomogeneities. Fractured poroelastic media show significant attenuation and velocity dispersion due to this mechanism. The theory has first been developed for the symmetry axis of the equivalent transversely isotropic (TI) medium corresponding to a poroelastic medium containing planar fractures. In this work, we consider the theory for all propagation angles by obtaining the five complex and frequency-dependent stiffnesses of the equivalent TI medium as a function of frequency. We assume that the flow direction is perpendicular to the layering plane and is independent of the loading direction. As a consequence, the behaviour of the medium can be described by a single relaxation function. We first consider the limiting case of an open (highly permeable) fracture of negligible thickness. We then compute the associated wave velocities and quality factors as a function of the propagation direction (phase and ray angles) and frequency. The location of the relaxation peak depends on the distance between fractures (the mesoscopic distance), viscosity, permeability and fractures compliances. The flow induced by wave propagation affects the quasi-shear (qS) wave with levels of attenuation similar to those of the quasi-compressional (qP) wave. On the other hand, a general fracture can be modeled as a sequence of poroelastic layers, where one of the layers is very thin. Modeling fractures of different thickness filled with CO₂ embedded in a background medium saturated with a stiffer fluid also shows considerable attenuation and velocity dispersion. If the fracture and background frames are the same, the equivalent medium is isotropic, but strong wave anisotropy occurs in the case of a frameless and highly permeable fracture material, for instance a suspension of solid particles in the fluid.     

Interactions among waves,current, and mud: Numerical and laboratory studies

Interactions between waves, current, mud and turbulence are very complicated in the coastal and estuarine turbid waters. It is still necessary to improve our understanding of the fundamental physical processes governing the cohesive sediment transport in the coastal and estuarine waters. A numerical model is developed to study the interactions among waves, current, and mud. An eddy viscosity model for wave and current is proposed in order to close the equations of wave motion or of current motion in a combined flow, respectively. The equations of mud transport are derived based on the visco-elastic properties of mud. Coupling the equations of wave motion or of current motion for water layer with those of mud layer can give (1) wave height; (2) distributions of current velocities in the water layer; (3) distributions of transport velocities at the water–mud interface; and (4) distributions of mass transport velocities within the mud layer. These modeled results are in a reasonable agreement with experimental results. Results suggest that (1) the rate of wave attenuation increases in the opposing currents (currents against in the direction in which the waves propagate) and decreases in the following currents (currents in the same direction as that in which the waves propagate); (2) the opposing currents would have more significant effects on the rate of wave height attenuation than the following currents; (3) the effect of current on the rate of wave attenuation on the muddy bottom is larger than that on the rigid bottom; (4) mud transport rate increased in the following currents but decreased in the opposing currents; and (5) the rate of wave height attenuation on the mud bottom is one order of magnitude larger than that on the rigid bottom.     

A field study of coastal dynamics on a muddy coast offshore of Cassino beach,Brazil

Mud deposits near sandy beaches, found throughout the world, are of scientific and societal interest as they form important natural sea defenses by efficiently damping storm waves. A multi-national field experiment to study these phenomena was performed offshore Cassino beach in southern Brazil starting in 2004. This experiment aimed to investigate the formation of an offshore mud deposit, to characterize wave attenuation over potentially mobile muddy bottoms, and to evaluate the performance of models for wave transformation over heterogeneous beds through the measurement of water waves, near-bottom currents, bathymetry, and changes in bottom sediment characteristics. The main instrumentation was a set of wave sensors deployed in a transect from the shoreline across sandy and muddy deposits offshore to a depth of 25 m. Additional sensors, including current meters and optical backscatter sensors, were concentrated at stations in the middle of the mud deposit and in the surf zone to document aspects of the wave boundary layer and lutocline dynamics. This fieldwork also involved the geological and geotechnical characterization of the mud deposit using seismic equipment, echo-sounders, cores, surficial sampling and an in-situ density meter. These sediment samples were subsequently analyzed for density, grain size distribution, mineralogy, rheology and sedimentary structures. In addition, video and radar monitoring equipment were installed to measure the long-term aspects of surf zone damping by fluid mud and any associated morphodynamic responses. This paper provides a summary of environmental conditions monitored during the experiment and describes the major findings of the various investigations. Although data collection was more difficult than anticipated and dramatic wave attenuation involving the onshore transport of fluid mud into the surf zone region was not observed during the instrumented interval, the new methodologies developed and comprehensive observations obtained during this effort are being used to improve our understanding of shoaling wave dynamics and sediment transport in the coastal zone in regions with significant cohesive sediment deposits.