Inversion of river-bottom sediment parameters using mechanically sampled specimens and subbottom profiling data

The study of river dynamics requires knowledge of physical parameters, such as porosity, permeability, and wave propagation velocity, of river-bottom sediments. To do so, sediment properties are determined on mechanically sampled specimens and from subbottom profiling. However, mechanical sampling introduces disturbances that affect test results, with the exception of grain-size distribution. In this study, we perform inversion of acoustic data using the grain-size distribution of mechanically sampled specimens and the relation between porosity and permeability from the Kozeny–Carman equation as prior information. The wave reflection coefficient of the water–silt interface is extracted from the raw subbottom profile. Based on the effective density fluid model, we combine the Kozeny–Carman equation and the wave reflection coefficient. We use experimental data from two Yellow River reservoirs to obtain the wave velocity and density of multiple sections and their spatial variations, and find that the inversion and testing results are in good agreement.     

Explicit use of the Biot coefficient in predicting shear-wave velocity of water-saturated sediments

Predicting the shear‐wave (S‐wave) velocity is important in seismic modelling, amplitude analysis with offset, and other exploration and engineering applications. Under the low‐frequency approximation, the classical Biot–Gassmann theory relates the Biot coefficient to the bulk modulus of water‐saturated sediments. If the Biot coefficient under in situ conditions can be estimated, the shear modulus or the S‐wave velocity can be calculated. The Biot coefficient derived from the compressional‐wave (P‐wave) velocity of water‐saturated sediments often differs from and is less than that estimated from the S‐wave velocity, owing to the interactions between the pore fluid and the grain contacts. By correcting the Biot coefficients derived from P‐wave velocities of water‐saturated sediments measured at various differential pressures, an accurate method of predicting S‐wave velocities is proposed. Numerical results indicate that the predicted S‐wave velocities for consolidated and unconsolidated sediments agree well with measured velocities.     

筑坝拦截对黑河河道沉积物粒度空间分布的影响

河流沉积物对流域环境变化具有敏感响应,其粒度参数能反映沉积环境中物质来源和水动力环境.本文以黑河流域上中游为研究区域,探究河流沉积物粒度对流域环境变化的响应.从黑河上中游干流22个主要控制断面采集河床沉积物样品,采用筛分法和吸管法对沉积物样品粒度参数进行测定,并分析其空间分布规律对筑坝拦截为主的环境改变响应.研究结果表明：受梯级水库建设影响,黑河上中游泥沙粒径大小差异显著.干流库区泥沙粒径较自然河段明显减小,分选很好,呈正偏或极正偏尖锐分布,而且在库区不同沉积高度上表现出分层沉积特征;坝下游河段因遭受强烈冲刷,较自然河段泥沙粒径粗化显著,分选变差,偏态趋向极正偏,峰态尖锐化;水库回水区受水库壅水及下泄清水的双重制约,泥沙粒度参数介于自然河段和坝下游河段之间,整体分选中等,呈极正偏尖锐分布;沉积环境分析表明,上游支流河段沉积物粒度特征受泥沙供给和物源特征的影响较水动力条件显著,干流河段沉积物粒度特征主要受水动力条件控制.研究结果既符合河流上中游沉积物粒度分布规律,也反映了河流环境变化对沉积物粒度组成的影响.     

Reflection and Transmission Coefficients for an Incident Plane Shear Wave at an Interface Separating Two Dissimilar Poroelastic Solids

Using Biot’s poroelasticity theory, we derive expressions for the reflection and transmission coefficients for a plane shear wave incident on an interface separating two different poroelastic solids. The coefficients are formulated as a function of the wave incidence angle, frequency and rock properties. Specific cases calculated include the boundary between water-saturated sand and water-saturated sandstone and the gas–water interface in sand. The results show a very different interface response to that of an incident P wave. Plane SV wave incidence does not significantly excite the Biot slow P wave if the frequency of the wave is below the transition frequency. Above this frequency, an incident plane SV wave can generate a mode-converted slow Biot P wave which is actually a normal propagating wave and not highly attenuating as in the usual (diffusive) case. For an incident SV wave onto a gas–water interface, even at very high frequency, there is no significant Biot second P wave produced. For small incident angles, the gas–water interface is essentially transparent. With increasing angles, there can arise an unusual "definitive angle" in the reflection/transmission coefficient curves which is related to the change of fluid viscosity on both sides of the interface and provides a possible new means for underground fluid assessment.     

模拟天然气水合物的岩石物理特性模型实验

针对水合物沉积的悬浮、颗粒接触和胶结三种微观模式,制作一组微弱胶结非固结高孔隙度人造样品和颗粒填充渐变的三维物理模型. 通过物理模型实验研究天然气水合物物性参数的敏感性. 实验结果表明：在弱颗粒间胶结物和低有效应力的固结差的沉积物中,声波对孔隙流体性质较敏感. 随着温度的降低颗粒胶结,改变原有沉积物的岩石物理特性,速度、弹性模量和频率升高,声波衰减和Vp/VS减小,沉积层内的反射波消隐.     

A new laboratory technique for determining the compressional wave properties of marine sediments at sonic frequencies and in situ pressures

We describe a new laboratory technique for measuring the compressional wave velocity and attenuation of jacketed samples of unconsolidated marine sediments within the acoustic (sonic) frequency range 1–10 kHz and at elevated differential (confining – pore) pressures up to 2.413 MPa (350 psi). The method is particularly well suited to attenuation studies because the large sample length (up to 0.6 m long, diameter 0.069 m) is equivalent to about one wavelength, thus giving representative bulk values for heterogeneous samples. Placing a sediment sample in a water‐filled, thick‐walled, stainless steel Pulse Tube causes the spectrum of a broadband acoustic pulse to be modified into a decaying series of maxima and minima, from which the Stoneley and compressional wave, velocity and attenuation of the sample can be determined. Experiments show that PVC and copper jackets have a negligible effect on the measured values of sediment velocity and attenuation, which are accurate to better than ± 1.5% for velocity and up to ± 5% for attenuation. Pulse Tube velocity and attenuation values for sand and silty‐clay samples agree well with published data for similar sediments, adjusted for pressure, temperature, salinity and frequency using standard equations. Attenuation in sand decreases with pressure to small values below Q^?1 = 0.01 (Q greater than 100) for differential pressures over 1.5 MPa, equivalent to sub‐seafloor depths of about 150 m. By contrast, attenuation in silty clay shows little pressure dependence and intermediate Q^?1 values between 0.0206–0.0235 (Q = 49–43). The attenuation results fill a notable gap in the grain size range of published data sets. Overall, we show that the Pulse Tube method gives reliable acoustic velocity and attenuation results for typical marine sediments.     

IDENTIFICATION OF SEA-FLOOR SEDIMENTS USING UNDERWAY ACOUSTICS*

Laboratory measurements, and a few in situ observations, show that saturated marine sediments have interdependent mechanical and acoustical properties, Acoustically, of particular importance are the acoustic impedance, velocity of sound and the sound attenuation coefficient of the sediment. The first two properties can be measured relatively easily from a surface ship; the measurement of attenuation however, is more problematical. It is suggested that this can be achieved by a quantitative treatment of the acoustic data collected during routine sub-bottom profiling over a variable thickness of superficial sediments. In the assessment of four different sediment locations in the Irish Sea it was found that quantitative treatment of the acoustic signals yielded both a value of the attenuation coefficient as well as a measure of the frequency dependence of the attenuation. In addition a statistical analysis of the signal intensities seems to provide an indication of the relative roughness of the bottom and sub-bottom interfaces. From the wide range of information provided the mechanical properties of the sea-floor sediment may be estimated.     

基于非饱和Biot-Stoll模型的海底沉积物介质声频散特性研究

利用含气非饱和Biot-Stoll模型研究了声波在海底表层沉积物介质中的传播,讨论了骨架耗散、含气饱和度对快纵波、慢纵波和横波速度和衰减的影响,并与Biot模型的结果进行了对比.研究结果表明:孔隙流体黏滞耗散与骨架耗散共同影响声波传播速度和衰减,低频情况下骨架耗散引起的衰减占主要地位,高频情况下骨架耗散引起的衰减较小;少量气体(<1%)的引入显著改变了快纵波速度,气体含量的变化对快纵波衰减影响很大,低频情况下气体对慢纵波速度的影响不大,而对横波速度的影响较大,气体含量的变化对慢纵波和横波衰减影响较小.利用超声波测量系统测量了一例杭州湾海底沉积物样品的纵波速度和衰减,当含气量趋近0%时,Biot-Stoll模型预测的纵波速度和实验测量结果较为一致.     

Acoustic wave propagation in saturated porous media: reformulation of the Biot/Squirt flow theory

A model of wave propagation in fluid-saturated porous media is developed where the principal fluid/solid interaction mode affecting the propagation of the acoustic wave results from the conjunction of the Biot and the Squirt flow mechanism. The difference between the original Biot/Squirt (BISQ) flow theory and the new theory, which we call the reformulated BISQ, is that the average fluid pressure term appearing in the dynamic equation for a two component solid/fluid continuum is independent of squirt flow length. P-velocity and attenuation relate to measurable rock physical parameters: the Biot's poroelastic constants, porosity, permeability, pore fluid compressibility and viscosity. Modelling shows that velocity and attenuation dispersion obtained using the reformulated BISQ theory are of the same order of magnitude as those obtained using the original BISQ theory. Investigation on permeability effect on velocity and attenuation dispersion indicate that the transition zone in velocity and attenuation peak, occurring both at the relaxation frequency, shifts toward high frequency when permeability decreases. This behaviour agrees with Biot's theory prediction.     

套管井声波全波列中地层波的幅度特性

套管井声全渡中的地层波可以检测固井第二界面的水泥胶结质量状况.地层波不仅受地层衰减、套管波干扰等因素影响,更主要受第二界面两侧的声阻抗差值的影响.为此,基于固体介质的声势函数和界面的边界条件,推导了固-固界面纵横波的反射系数和透射系数,以及第一临界角入射情况下反射系数和透射系数的近似公式,并且分析了透射系数与声阻抗之间的关系.数值计算结果发现,地层波幅度随第二界面两侧声阻抗差增大而减小,且这种关系是非线性的.基于临界入射的透射系数,对经过衰减校正和干扰校正的地层波进行声阻抗校正,消除地层声阻抗对地层波幅度影响.使用地层波幅度评价第二界面的水泥胶结状况,并在实际资料处理中见到了好的效果.     

深海海底沉积物声学特性原位测量试验研究

由于实验室测试环境条件与深海海底的原位温压环境存在较大的差异,取样测量的沉积物声学参数通常偏离海底原位状态的真实值.该文利用深水型压载式沉积声学原位测量系统,在水深超过5000 m的西太平洋海域开展了沉积物声学原位测量试验,准确获取了深海底原位状态下沉积物的声速和声衰减系数,并同步采集了沉积物柱状样品.结果显示,实验室测量的沉积物声速、声衰减系数均高于原位测量结果.通过温压校正和数据—模型对比表明,实验室声速比与原位声速比存在较大的差异,难以通过温压校正模型和频散模型完全改正,这可能与深海沉积物的结构扰动有关.较之浅海沉积物,软弱的深海沉积物结构更易受到采样过程的扰动,因而深海沉积物的声学特性测量更为困难,这表明在深海海底声学研究采用原位测量技术具有必要性.该研究标志着我国海底声学研究由浅海走向深海,对推动海底声学的深入研究具有重要意义.     

Numerical simulation of reservoir sediment and effects on hydro‐dynamic response of arch dams

Based on the dynamic theory for saturated porous media by Biot (Journal of the Acoustical Society of America 1956; 28 : 168–178), a numerical model is presented to analyse the reflection behaviours of reservoir sediment and compared with those from the visco‐elastic model. It is concluded that the two models give very similar results of reflection coefficient α within the frequency range of interest. Then, using the two models, the change of the reflection coefficients α with various sedimentation parameters and excitation frequencies are studied in detail. The results are further used in the analysis of response functions of hydro‐dynamic pressures on, and structural displacements of the Xiang Hong Dian arch dam, for which some results from a field vibration test are available. It appears that effects of water compressibility with sediment reflection on hydro‐dynamic pressures and structural response are not significant for this specific case. Copyright © 2001 John Wiley & Sons, Ltd.     

The paleoenvironmental significance of spatial distributions of grain size in groundwater‐recharged lakes: A case study in the hinterland of the Badain Jaran Desert,northwest China

For lakes in desert hinterlands that are not recharged by river runoff, sediment input solely comes from wind transport. While the processes of sediment transport and deposition in these lakes differ significantly from those with river discharge, the spatial distribution of sediment grain size in these groundwater‐recharged lakes remains largely unknown. Moreover, whether the grain size distribution in these lake sediments can be used as a proxy in the study of past climatic change and environmental evolution studies is unclear. In this study, five lakes with a range of surface areas that had no runoff recharge were selected from the hinterland of the Badain Jaran Desert of north‐western China, and a total of 108 samples of lake surface sediments were collected to examine the spatial distribution of grain size. Moreover, an end‐member‐modeling algorithm was used to calculate end members from all grain size measurements. Our results showed that both the median and mean grain sizes in the lake sediments decreased from the nearshore to the offshore, deep‐water zone. However, the lowest median and mean grain sizes were not found in the center of the lakes, in contrast to lakes recharged by surface runoff. The median grain size of sediment in the lake center was negatively correlated with lake level, and thus could help reveal lake evolution at low resolutions. Moreover, EM1 and EM2 were interpreted as wind transported sediment, and sediment perturbed by lake waves after wind transport, respectively. The modal grain size of EM1 varied slightly between lakes, while changes in the modal grain size of EM2 were related to lake area. Given the positive relationship found between EM2 content and lake level, changes in the EM2 content (%) can serve as a rough indicator of lake level fluctuations at low temporal resolutions. Copyright © 2017 John Wiley & Sons, Ltd.     

Extension of White's layered model to the full frequency range

The low‐frequency theory of the White model to predict the dispersion and intrinsic attenuation in a single porous skeleton saturated with periodic layers of two immiscible fluids is extended to the full frequency range using the Biot theory. The extension is similar to the Dutta–Odé model for spherical inhomogeneities. Below the layer resonance frequency, the acoustic bulk properties for several gas–water fractions are in good agreement with the original White model. Deviations start to occur at higher frequencies due to the growing importance of resonance phenomena that were neglected in the original White model. The full model predicts significantly higher damping at sonic frequencies than the original White model. We also show that attenuation is significantly dependent on porosity variations. With realistic rock and fluid properties, a maximum attenuation of about 0.3 is found at seismic frequencies.     

一种新型海底沉积物声学原位测量系统的研制及应用

为了能够精确地测量海底表层沉积物的声学参数,自主研制了一种新型海底沉积物声学原位测量系统,与国内外传统的声学原位测量系统相比,该系统能够实时显示声波波形,调整测量参数,其工作方式除了站位式测量之外,还实现了拖行式连续测量,极大地提高了工作效率.根据前期海试情况,对海底仪器结构进行了重新设计,使之可以同时测量海底沉积物及海底海水的声学参数,同时建立了双向数字信道,解决了测量过程中系统信号的干扰问题.该系统的结构分为两部分:甲板控制单元和水下测量单元,整套系统通过主机控制程序进行控制,采用GPS定位系统测定仪器的大地坐标.为了检验系统的稳定性及准确性,分别进行了实验室水槽实验和海上试验.利用水声测量设备对测量系统进行实验室水槽标定分析,实验结果表明系统测量值相对误差仅为0.04%,测量结果具有较高的精度.海上试验在青岛胶州湾和东海海域进行,获得了试验区域海底沉积物声速和声衰减系数的测量数据,将测量数据与他人的研究结果进行对比分析,结果表明测量数据与前人研究结果一致,较为准确.该原位测量系统在站位式测量和拖行式测量中都能够快速准确地测量出沉积物声速和声衰减系数,可以作为海底底质声学测量的调查设备.     

Variations in texture of beach sediments in the vicinity of the Tirumalairajanar river mouth of India

The distribution of grain size parameters along 11 km stretch of the beach sediments between Karikal and Nagore,reveals that the mean grain size exhibits a marked decreasing trend on either side of the mouth of the Tirumalairajanar River which flow from west to east.The sediments are mainly of medium to coarse grained,moderately sorted,near-symmetrical skewed to fine skewed and leptokurtic to mesokurtic in nature.Interrelationship of various parameters shows bimodal nature of sediments having dominance of medium to coarse sand.The major part of the sediment fall in a coarse to fine grained category(sand and silt).Based on the CM(Coarser one percentile value in micron) pattern,the sediment fall in rolling and suspension field.These factors includes the sediments discharged from the river mixes with offshore sediments and with the sediments eroded from a source rock.The effect of wave sorting, and the northward drifting of sediments by littoral current are understandable.Results indicate that the Tirumalairajanar River is the most important source for modern sediments in the study area.The agitation by waves is an important sorting mechanism in the study area,and the net sediment transport in the study area is northward.The findings are based on the grain sizes and also corroborated by shortterm observations of the beach sediment dynamics and transport during the monsoon and summer seasons between Karaikal and Nagore region.     

Particle size characteristics of suspended and bed sediments in the rhône river

Ten large volume water samples were taken from the Rhône River (Switzerland-France) in November, 1989 for recovery of total suspended sediment by continuous flow centrifugation. the samples were freeze-dried and analysed for particle size, organic carbon, total nitrogen, and carbonate. for comparative purposes, four bed sediments collected in July, 1989 are also described. the Rhône can be subdivided into three sections on the basis of the origins of the water. the first section is the Upper Rhône River draining into Lake Geneva. Waters are derived from glaciers, with low temperature and conductivity and high turbidity. Suspended sediment is coarse, has a bimodal distribution, and is low in both organic matter and carbonate. the second reach is from Lake Geneva to the confluence with the Saône at Lyon and has warmer water with higher conductivity and very low turbidity. Suspended sediment is higher in organic matter, with high carbonate originating from the lake. the final section is from Lyon to Arles, with warmer water and higher conductivity and turbidity due to modification by the Saône. Sediment is rich in organic matter, which May, account for an observed decline in oxygen in the river waters downstream from Lyon. Carbonate in these sediments also decreases due to increased turbidity from the Saône. Suspended sediments other than from the Upper Rhône show a remarkable consistency in grain size, predominantly in the fine silts (mode 9-11 μm). This consistency indicates a high degree of suitability for geochemical analysis. Bed sediments were bimodal throughout, with a dominant coarse population in two out of the four samples. Grain size statistical parameters could be easily explained by application of the theory of mixing of two major populations in the sand size (bed traction load) and the fine silt/clay size (suspended sediment load).     

Determining hydraulic resistance in gravel‐bed rivers from the dynamics of their water surfaces

Traditionally, approaches to account for the effect of the boundary roughness of a gravel‐bed river have used a grain‐size index of the bed surface as a surrogate for hydraulic resistance. The use of a single grain‐size does not take into account the spatial heterogeneity in the bed surface and how this heterogeneity imparts resistance on the flow, nor the way in which this relationship changes with variables such as flow stage. A new technique to remotely quantify hydraulic resistance is proposed. It is based on measuring the dynamics of a river's water surface and relating this to the actual hydraulic resistance created by a rough sediment boundary. The water surface dynamics are measured using a new acoustic technique, grazing angle sound propagation (GRASP). This proposed method to measure hydraulic resistance is based on a greater degree of physical reasoning, and this is discussed in the letter. By measuring acoustically the temporal dynamics of turbulent water surfaces over a water‐worked gravel bed in a laboratory flume, a dependency is demonstrated between the temporal variation in the reflected acoustic pressure and measured hydraulic resistance. It is shown that the standard deviation in acoustic pressure decreases with increasing hydraulic resistance. This is shown to apply for a range of relative submergences and bed slopes that are typical of gravel‐bed rivers. This remote sensing technique is both rapid and inexpensive, and has the potential to be applied to natural river channels and to other environmental turbulent flows, such as overland flows. A whole new class of low‐cost, remote and non‐intrusive instruments could be developed as a result and used in a wide range of hydraulic and hydrological applications. Copyright © 2006 John Wiley & Sons, Ltd.