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

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

南海北部海底沉积物在温度变化下的三种声速类型

研究南海北部海底沉积物温度变化状态下声速性质,得出以下结论:(1) 南海北部海底沉积物具有声速温度正增长(STPIK)、声速温度负增长(STNIK)和声速温度波动(STWK)三种类型,后两种类型在世界其他范围内海域未见报道.(2) 声速温度正增长类型和声速温度负增长类型沉积物的温度变化对声速影响都具有十分显著的线性关系,但是原状样品由于组成不均匀性,增大了声速的非线性变化.(3)南海北部海底沉积物的孔隙度、密度等主要物理参数差异不明显,难以直接解释三类样品的温度-声速性质的不同.(4)对于STPIK类型沉积物,理论分析与实验结果的统一性,可以运用海底沉积物与海水声速比进行校正不同温度状态下的海底沉积物的声速.(5)对于STNIK和STWK类型沉积物,需要深入研究,从理论和实验角度揭开其机理和成因.海底沉积物声速-温度特性研究将为提高南海北部海域天然气水合物声学探测精度和准确度提供声速性质依据.     

海底大地基准网建设及其关键技术

海底大地基准网建设涉及基准信标方舱研制、网型设计、勘选布放、观测策略、观测模型建立与优化、数据处理策略等.文章试图较全面地描述海底大地基准网建设的主要技术问题,并探索其技术途径;分别描述了浅海海底方舱和深海海底方舱建设的基本构想,提出了适用于水下定位导航的海底大地基准网整体建设准则和用于地壳运动监测的整体网型设计准则.为了提高海底大地基准网的应用效能,海底网型以对称网型为主;海底基准点定位的海面观测也以对称观测构型为主,建议采用圆形观测和交叉十字(或井字型)观测相结合的海底定位模式; GNSS相位中心与声学换能器的偏移改正,建议采用外部测定与模型参数估计相结合的方法进行;声速误差影响改正建议采用观测值改正与参数化模型改正相结合的方式,并以模型改正为主削弱声线误差对海底定位的影响.利用文章的基本设计,分别进行了浅海和深海定位试验.在浅海充分验证海底方舱性能的基础上,建立了3000m水深的长期海底基准点,该基准点初步定位结果内符精度优于5cm.     

基于GNSS-A的海洋声速变化估计及其对定位的影响

声速时空变化是影响水下定位精确度的重要因素.首先，根据声线跟踪原理推导出了包含声速扰动的观测方程，并采用三次B样条曲线对声速扰动时空模型进行建模；之后，提出利用最小范数条件选择最优超参数，实现海底应答器坐标与声速时空模型的联合估计.采用2019年7月3000 m水深实测声学数据计算了海底站点的坐标、声速的连续时间变化与二维空间梯度.结果表明，3组数据集之间E/N/U分量的最大差异分别为0.391,0.035和0.120 m,海底坐标解算内符精度达到cm级；声速时间变化小于1 m·s^-1,呈现明显的日变化和内波引起的短周期波动特性；声速空间变化小于0.1 m·s^-1,空间梯度小于0.04 m·s^-1·km^-1;不同数据集定位结果的差异是由于声速的周期性时间变化引起的.     

南海海底钙质土声学物理性质及其工程地质意义

分析和讨论了南海热带海域海底钙质土基本特征和声学物理性质 ,它们与以陆源物质为主的沉积物在声学物理性质方面截然不同 .一般而言 ,钙质土具有高含水量、高孔隙度、低湿密度、较高声速以及抗压强度变化较大的特征 ,主要是由碳酸盐颗粒的海洋生物残骸等所组成     

利用地震波速研究青藏高原东北缘地壳组成及其动力学

依据穿过巴颜喀拉地块的北部、秦岭地块、祁连地块、海原弧形构造区和鄂尔多斯地块的玛沁-兰州-靖边人工地震剖面的P波、S波的速度结构和泊松比结构,对青藏高原东北缘的地壳组成进行研究,并探讨其动力学过程. 首先,系统地归纳总结出一套将地震测深得到的原位P波速度校正到实验室温压条件下波速的具体可行的方法,利用大地热流值求取地壳不同深度的温度是该方法的关键. 然后,将上述剖面的原位P波速度校正到600 MPa和室温条件下,结合泊松比与相同温压条件下的实验室岩石波速测量结果进行对比,确定研究区的岩性组成. 结果表明,青藏高原东北缘地壳平均P波校正波速为6.43 km/s,地壳整体像上地壳一样呈酸性. 巴颜喀拉地块和秦岭地块南部的下地壳底部缺失校正速度V_p>6.9 km/s的基性岩,下地壳中酸性互层,下地壳整体呈酸性. 其他地块下地壳底部有2～10 km厚的校正速度V_p>6.9 km/s的基性岩,下地壳整体呈中性. 最后,根据青藏高原东北缘地壳结构和组成的研究成果,支持地壳增厚主要发生在下地壳的观点;提出巴颜喀拉地块和秦岭地块南部曾发生过下地壳拆沉作用,并导致高原的加速隆升.     

基于改进遗传算法的BISQ模型多参数反演方法研究

BISQ模型反演在油气资源与海洋地震勘探中具有广阔的应用前景与实用价值.为解决具高维参数与多峰值的BISQ模型参数反演问题,本文在传统的遗传算法基础上添加插入了全局保优机制、Boltzmann生存机制、小生境之DC漂移技术和群体灾变等方法与技术,使改进组合的遗传算法能在高维参数与多峰值的反问题中搜索得到全局最优解.结合海底多频原位沉积物声学测量数据,应用本文方法,得到了与声波传播速度相似性最优模型的各参数值,效果明显.     

基于弯曲元技术的含水合物松散沉积物声学特性研究

含水合物松散沉积物的声学特性对海上天然气水合物地球物理勘探和资源评价具有重要意义.研制了适用于高压条件下含水合物沉积物声学特性探测的纵横波一体化新型弯曲元换能器,提出利用频谱分析(FFT)和小波分析(WT)相结合的方法获取纵横波速度,并进行了多个轮次的水合物声学特性模拟实验研究.结果表明,新型弯曲元技术可以灵敏探测松散沉积物中水合物的生成和分解,随着水合物饱和度(Sh)的增大,纵横波速度呈规律性增长:当Sh<25%时,纵横波速度增长较快,水合物可能胶结沉积物颗粒生成;25%~60%之间,声速增长较为缓慢,水合物可能与沉积物颗粒呈接触关系;在Sh>60%时声速随着水合物饱和度增加又快速增长,表明水合物可能重新胶结沉积物颗粒生成.     

被动源海底地震仪数据预处理技术构建与应用

被动源海底地震仪(OBS)探测是海洋深部结构研究的重要方法之一.受投放区洋流、海底地形崎岖,以及海底温压环境等多种因素影响,OBS数据预处理面临时钟漂移,姿态及水平方位角偏差等问题.本文基于国产OBS硬件结构特点与其观测数据特性,针对仪器姿态校正、时间校正和水平方位角校正三个关键环节构建技术解决方案.通过硬件架构特征构建仪器姿态校正方案;结合仪器记录的石英晶振特点和背景噪声互相关方法提出时间校正技术;基于最小化P波切向分量能量和P波主成分分析两个技术方法进行水平方位角校正.上述预处理技术体系在苏拉威西海域被动源OBS数据预处理中取得了良好的应用效果,综合背景噪声水平与地震信号的时频分析结果,表明该预处理技术能有效提高被动源OBS观测数据质量,为进一步的科学研究奠定坚实基础.     

ROV专用的海底热流探测技术与应用研究

地热流资料是用于研究油气、天然气水合物矿产资源的控制与分布,评价油气矿产资源的主要依据,为了了解水合物勘探区地温场的结构、状态,需要测量水合物勘探区更精细的热流.本文通过介绍ROV专用的海底热流探测技术的特点、海底地温梯度测量作业的技术与方法,并结合我国深海探查的工作特点,讨论该项技术的适用范围、应用领域和前景.     

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.     

海底天然气渗漏的地球物理特征及识别方法

海底天然气渗漏是在海洋环境中广泛分布的自然现象.海底天然气渗漏可以指示沉积层中的烃类聚集带,渗漏出的大量气体（主要是甲烷）可能影响全球的气候变化,此外,与海底渗漏相关的浅层气改变了海底沉积物的土工性质,可能对海底工程构成威胁.因此海底渗漏研究意义重大.海底天然气渗漏不仅影响海底沉积物的物理性质,而且还极大地改变海底地形地貌,它能在海底形成麻坑、泥火山、冷泉碳酸盐岩以及化学自养生物群落等现象.在海底渗漏发生的地方,地形地貌特征可以在海洋测深和逆向散射数据上得到反映;沉积物的声学特征可以在地震剖面上得到反映,如形成声混浊、声空白、亮点、多次波、速度下拉等;有些渗漏在海面形成油渍膜,油渍膜可以在合成孔径雷达图像上得到反映.根据海底渗漏的上述地球物理特征,可以识别出可能渗漏区域,海底渗漏的证实需要用到海底观测和取样分析资料.     

Seafloor geodetic network establishment and key technologies

Seafloor geodetic network construction involves the development of geodetic station shelter, network configuration design, location selection and layout, surveying strategy, observation model establishment and optimization, data processing strategy and so on. This paper tries to present main technological problems involved in the seafloor geodetic network construction, and seek the technically feasible solutions. Basic conceptions of developing seafloor geodetic station shelters for shallow sea and deep-sea are described respectively. The overall criteria of seafloor geodetic network construction for submarine navigation and those of network design for crustal motion monitoring are both proposed. In order to enhance application performances of the seafloor geodetic network, the seafloor network configuration should prefer a symmetrical network structure. The sea surface tracking line measurements for determining the seafloor geodetic station position should also adopt an approximately symmetrical configuration, and we recommend circle tracking line observations combined with cross-shaped line(or double cross-shape line) observations for the seafloor positioning mode. As to the offset correction between the Global Navigation Satellite System antenna phase center and the acoustic transducer, it is recommended to combine the calibration through external measurements and model parameter estimation. Besides, it is suggested to correct the sound speed error with a combination of observation value correction and parameterized model correction, and to mainly use the model correction to reduce the influence of acoustic ray error on the seafloor positioning. Following the proposed basic designs, experiments are performed in shallow sea area and deep-sea area respectively. Based on the developed seafloor geodetic shelter and sufficient verification in the shallow sea experiment, a long-term seafloor geodetic station in the deep-sea area of 3000 m depth was established for the first time, and the preliminary positioning result shows that the internal precision of this station is better than 5 cm.     

The influence of gas bubbles on sediment acoustic properties: in situ, laboratory, and theoretical results from Eckernförde Bay, Baltic sea

Acoustic turbidity caused by the presence of gas bubbles in seafloor sediments is a common occurrence worldwide,but is as yet poorly understood. The Coastal Benthic Boundary Layer experiment in the Baltic off northern Germany was planned to better characterize the acoustic response of a bubbly sediment horizon. In this context, in situ measurements of compressional wave speed and attenuation were made over the frequency range of 5–400 kHz in gassy sediments of Eckernförde Bay. Dispersion of compressional speed data was used to determine the upper limit of the frequency of methane bubble resonance at between 20 and 25 kHz. These data, combined with bubble size distributions determined from CT scans of sediments in cores retained at ambient pressure, yield estimates of effective bubble sizes of 0.3–5.0 mm equivalent radius. The highly variable spatial distribution of bubble volume and bubble size distribution is used to reconcile the otherwise contradictory frequency-dependent speed and attenuation data with theory. At acoustic frequencies above resonance (>25 kHz) compressional speed is unaffected by bubbles and scattering from bubbles dominates attenuation. At frequencies below resonance (<1 kHz) ‘compressibility effects’ dominate, speed is much lower (250 m s^-1) than bubble-free sediments, and attenuation is dominated by scattering from impedance contrasts. Between 1.5 and 25 kHz bubble resonance greatly affects speed and attenuation. Compressional speed in gassy sediments (1100–1200 m s^-1) determined at 5–15 kHz is variable and higher than predicted by theory (<250 m s^-1). These higher measured speeds result from two factors: speeds are an average of lower speeds in gassy sediments and higher speeds in bubble-free sediments; and the volume of smaller-sized bubbles which contribute to the lower observed speeds is much lower than total gas volume. The frequency-dependent acoustic propagation is further complicated as the mixture of bubble sizes selectively strips energy near bubble resonance frequencies (very high attenuation) allowing lower and higher frequency energy to propagate. It was also demonstrated that acoustic characterization of gassy sediments can be used to define bubble size distribution and fractional volume.     

Acoustic parameters inversion and sediment properties in the Yellow River reservoir

The physical properties of silt in river reservoirs are important to river dynamics. Unfortunately, traditional techniques yield insufficient data. Based on porous media acoustic theory, we invert the acoustic parameters for the top river-bottom sediments. An explicit form of the acoustic reflection coefficient at the water–sediment interface is derived based on Biot’s theory. The choice of parameters in the Biot model is discussed and the relation between acoustic and geological parameters is studied, including that between the reflection coefficient and porosity and the attenuation coefficient and permeability. The attenuation coefficient of the sound wave in the sediments is obtained by analyzing the shift of the signal frequency. The acoustic reflection coefficient at the water–sediment interface is extracted from the sonar signal. Thus, an inversion method of the physical parameters of the riverbottom surface sediments is proposed. The results of an experiment at the Sanmenxia reservoir suggest that the estimated grain size is close to the actual data. This demonstrates the ability of the proposed method to determine the physical parameters of sediments and estimate the grain size.     

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.     

Sediment shear properties from seafloor compliance measurements: Faroes-Shetland basin case study

Shear properties provide important information about the lithology, fluid content and stability of sediments but are difficult to measure using conventional seismics in the marine environment. Seafloor compliance measurements are sensitive to subsurface shear properties but have only been used in the Pacific Ocean and on shallow coastal shelves, where the source wave energy is known to be strong. We show here that seafloor compliance measurements can provide useful information about shear properties of marine sediments in less energetic settings and under high noise conditions caused by strong seafloor currents. We measured compliance at three sites in the Faroes‐Shetland sedimentary basin north of the Atlantic ocean. The sites have 1000 times higher noise levels than quiet seafloor sites and the source wave power is highly variable, but the data still reveal significant differences in sediment properties between two sites down to 2 kilometres beneath the seafloor. The first site, at the northern end of the basin, has an average shear velocity of 400 m/s in the upper 0.6 kilometres beneath the seafloor, increasing to approximately 2100 m/s at 2 kilometres beneath the seafloor. The second site, further south and to the west of the basin axis, has an average shear velocity of 150 m/s in the upper 0.6 kilometres beneath the seafloor, increasing to 1400 m/s at 2 kilometres beneath the seafloor. The sediments are probably unconsolidated in the upper 0.6 kilometres beneath the seafloor at both sites, with a mean grain size of 1 μm at the southern site and 20 μm at the northern site. The southern site has higher porosity at all depths and a higher risk of borehole collapse during drilling.     

基于合成孔径声学深拖调查的海底浅表层流体活动研究——以SAMS DT6000深拖在琼东南海域调查为例

声学深拖作为一个声学设备搭载平台,主要功能是获取高分辨率的声学数据,精细刻画海底地形地貌特征以及浅层剖面结构,对于研究海底浅表层流体活动系统的类型、形成机制和演化模式有着重要作用.本文介绍的合成孔径声学深拖（Synthetic Aperture Sonar Deep-tow）搭载了合成孔径声呐、浅地层剖面仪以及多波束系统等声学设备,相比于传统的侧扫声呐,合成孔径声呐采用小物理孔径基阵通过信号处理虚拟合成大孔径基阵来获得方位向高分辨率,大大提高了测绘速率,同时结合高分辨率的浅地层剖面和多波束背散射数据,可实现海底浅表层特征的三维立体显示.为查明调查区海底浅层流体活动的声学特征,分析天然气水合物相关的流体渗漏活动性与浅层构造之间的关系,我们利用声学深拖对研究区进行了全覆盖的扫测,获得了高分辨率的合成孔径声呐图像、浅地层剖面资料以及多波束背散射数据,平面上识别出多个呈条带状的海底丘状体,火焰状的流体渗漏,新月形的麻坑构造等流体活动地质构造;浅层剖面上可见气体聚集的声学空白段落,凸起的活跃喷口,以及反射杂乱的柱状浑浊带.通过识别流体活动的特征,我们总结了浅层流体活动演化模式具有周期性：游离气体通过高渗透运移通道上升至海底,首先扩散聚集造成局部沉积物体积膨胀形成丘状体;然后受其各种外界因素影响丘状体崩塌而引起气体渗漏;最后流体逸散剥蚀海底松散沉积物而形成麻坑构造;随着流体排出,喷口重新闭合,流体在地层中再次聚集,聚集的气体又将沉积地层上拱,在麻坑底部又可能生成含气丘状体.海底浅表层蕴藏着丰富的地质信息,这对于研究海底复杂的流体活动有着重要意义.