Estimation of Geoacoustic Properties of Marine Sediment Using a Hybrid Differential Evolution Inversion Method

This paper reports the inversion of midfrequency (1500–4500 Hz) chirps from a short-range transmission experiment conducted on the New Jersey Continental Shelf during the 2006 Shallow Water Experiment (SW06). The source was held at different depths and the sound signals were recorded at a vertical line array to investigate the interactions with the sea bottom at various grazing angles. Strong reflections from the sediment layer were seen in the data for all of the sources. Due to the presence of complex microstructures in the thermocline of the oceanic sound-speed profile, fluctuations both in amplitude and arrival time of the direct path arrivals were observed. Time variation of the water-column environment was also evident during the source transmissions. To mitigate the effects of the ocean environment on the seabed property estimation, a multistage optimization inversion was employed. The sound speed and the experimental geometry were inverted first using only the travel times of the water-column arrivals. The bottom sound speed and the sediment layer thickness were then inverted by matching the travel times of the bottom and sub-bottom reflections. The average of the estimated values for the sediment sound speed is 1598 m/s, consistent with in situ measurements from other experiments in the vicinity.      

Geoacoustic Inversion for the New Jersey Shelf: 3-D Sediment Model

A 3-D model of sediment sound speed for a 90-km $^{2}$ area on the New Jersey shelf was constructed by application of a geoacoustic inversion technique. This approach is based on a combination of seismic reflection measurements and a perturbative inversion scheme using horizontal wave number estimates. In a two-step process, seismic reflection measurements were used to identify depths at which discontinuities in the sound-speed profile (SSP) likely occur. Then, the perturbative inversion algorithm made use of this a priori information by employing qualitative regularization, an optimal method for addressing stability and uniqueness issues associated with solving the ill-posed inversion problem that provides for resolution of the layered seabed structure.      

Estimation of Transmission Loss in the Presence of Geoacoustic Inversion Uncertainty

A common problem in sonar system prediction is that the ocean environment is not well known. Utilizing probabilistic based results from geoacoustic inversions we characterize parameters relevant to sonar performance. This paper describes the estimation of transmission loss and its statistical properties based on posterior parameter probabilities obtained from inversion of ocean acoustic array data. This problem is solved by first finding an ensemble of relevant environmental model parameters and the associated posterior probability using a likelihood based inversion of the acoustic array data. In a second step, each realization of these model parameters is weighted with their posterior probability to map into the transmission loss domain. This approach is illustrated using vertical-array data from a recent benchmark data set and from data acquired during the Asian Seas International Acoustics Experiment (ASIAEX) 2001 in the East China Sea. The environmental parameters are first estimated using a probabilistic-based geoacoustic inversion technique. Based on the posterior probability that each of these environmental models fits the ocean acoustic array data, each model is mapped into transmission loss. This enables us to compute a full probability distribution for the transmission loss at selected frequencies, ranges, and depths, which potentially could be used for sonar performance prediction.     

Results from the Geoacoustic Inversion Techniques Workshop using modal inverse methods

This paper presents the results of a perturbative inverse approach applied to the range-dependent acoustic data provided as part of the Geoacoustic Inversion Techniques Workshop. The method is based on the Hankel transform relationship between complex-pressure field data measured on a horizontal array and the depth-dependent Green's function for a horizontally stratified medium. The input data to the inversion algorithm are discrete values of horizontal wavenumbers estimated for locally range-independent subapertures of the acoustic data. Inversion results are presented with emphasis placed on inverting for compressional wave speed as a function of both range and depth in the bottom.     

Modal Inversion Analysis for Geoacoustic Properties of the New Jersey Continental Shelf in the SWAT Experiments

In this paper, inversion for bottom sediment properties at a site on the New Jersey continental shelf is studied as part of the Shallow Water Acoustic Technology (SWAT) project. A source towed at a constant water depth over a range of some tens of kilometers transmitted low-frequency continuous wave (cw) signals, which were measured on a bottom-moored vertical line array of receivers. For the along-shelf geometry, the zeroth-order asymptotic Hankel transform is then applied to the acoustic field at 50 Hz measured on the resulting synthetic aperture horizontal array created at each receiver depth. The resulting horizontal wave number spectra, which have peaks corresponding to the mode eigenvalues, are observed to have slightly different values at different receiver depths, and therefore, stochastic mode inversion is exploited to utilize all of the observed peak position information. The estimated sound-speed profile (SSP) for the upper 10 m of sediment is then compared with an inversion result obtained using midfrequency (2–16 kHz) chirp sonar pulses reflected at normal incidence from the sediment. Although obtained using totally different inversion techniques, both estimated profiles are shown to be in good agreement in the top 10 m of sediment. The acoustic field simulated using the inverted SSP also agrees well with the measured acoustic field at each receiver depth. Furthermore, simulated sound fields which use this profile as input data are shown to be effective in predicting the measurements obtained at a different frequency (125 Hz) and for a different (cross-shelf) geometry.      

聊城市人工增雨效果统计方法及计算结果分析

利用聊城市和菏泽市1963～2007年的降水观测资料和天气形势图,采用历史区域回归试验和天气系统分类回归方案效果检验法,对2004～2007年4年间聊城市人工增雨作业的效果进行了计算分析,为聊城市人工增雨作业选择有利的作业时机和天气条件提供一定的依据.结果表明:1)采用天气系统分类统计的方法可以有效提高区域雨量的相关性,从而提高了统计分析的效率;2)聊城市实施人工增雨作业比较有利的时节是6月,催化效果最理想的天气系统是南方气旋;3)人工增雨效果的最大值在不同的天气系统条件下对应着不同的自然雨强区间.     

Inversion Method for Sound Velocity Profile of Eddy in Deep Ocean

The modal wave number tomography approach is used to obtain sound speed profile of water column in deep ocean. The approach consists of estimation of the local modal eigenvalues from complex pressure field and use of these data as input to modal perturbative inversion method for obtaining the local sound speed profile. The empirical orthonormal function (EOF) is applied to reduce the parameter search space. The ocean environment used for numerical simulations includes the Munk profile as the unperturbed background speed profile and a weak Gaussian eddy as the sound speed profile perturbation. The results of numerical simulations show the method is capable of monitoring the oceanic interior structure.     

一种改进的Curvelet变换遥感影像融合方法

提出了一种改进的Curvelet变换的融合算法,以IKONOS影像为样本数据,融合过程中将原全色波段与多光谱影像中未参与HIS分解的波段进行加权组合,形成新全色分量,再将此分量与明度分量分别进行第二代Curvelet变换,以形成新明度分量,HIS逆变换后形成融合影像。并与改进的HIS变换融合、传统的HIS变换与Curvelet变换结合的融合算法相比较,结果表明,改进的方法能增加融合影像的信息量,降低光谱扭曲度,提高融合影像与原多光谱影像的相似度。     

利用海岸带遥感图像提取岸线的小波变换方法

小波变换作为一门新兴的数学分支已被广泛应用于图像处理领域。本文将小波变换应用于黄河三角洲遥感图像的边缘提取 ,检测出图像的边缘 ,从而得到了三角洲岸线信息。实验结果表明基于小波变换的图像边缘提取要优于经典边缘算子的提取。此方法对于把握河口三角洲的冲淤演变规律和海岸带开发具有重要意义     

Argo浮标CTD温度传感器实验室中测试评估方法

文中主要研究Argo用CTD在实验室内进行海洋工作状态的模拟。由于Argo用CTD温度传感器布放在深海区域,在高压低温的状态下,在剖面测量中曾经出现低温温度跃迁的现象,因此在CTD实验室对其温度性能测试评估方法加以改进。具体方法如下:先以常规测试评估方法检验,将已经确定达到指标要求温度传感器在压力罐内反复加压、保压、泄压,对其进行长期温度性能测试。根据实验数据总结出Argo浮标用CTD温度传感器的海洋环境温度测量特性及具体测试评估方法,使实验室温度测试水平得到提升。     

深基坑内支撑等效刚度数值计算影响因素分析

深基坑平面内支撑结构等效刚度的计算常采用有限元等数值计算方法,通过分析深基坑支护结构的4种典型平面内支撑结构型式,探讨了数值计算模型中内支撑结构的单元划分长度、约束弹簧刚度、约束弹簧间距、约束弹簧布置型式等因素对内支撑等效刚度计算结果的影响.经计算分析可知,在其他影响因素不变而只改变某个影响因素时,平面内支撑结构的等效刚度随腰梁单元长度的增加、约束弹簧刚度的增加呈增加的趋势,随约束弹簧间距的增加呈降低的趋势,但基本不受约束弹簧布置型式、除腰梁之外的其他内支撑构件单元长度的影响.在此基础上,得出了数值计算中上述各影响因素的建议取值范围和约束弹簧的布置型式.     

基于等效声速剖面法的多波束测深系统声线折射改正技术

多波束测深技术是近几年比较流行的海洋勘测手段之一,其具有全覆盖、无遗漏、高精度和高效率等特点。多波束测深系统的声学原理和海水的不均匀性,使得声波在传播过程中发生声线的折射现象,对波束测点的最终位置的归算带来较大的误差。针对多波束数据后处理方法,介绍了利用等效声速剖面法对多波束测深系统的声线折射进行改正的一种新方法,并通过试验对比,证明了该方法可以提高多波束测深系统的整体测量精度。     

基于条件价值评估法的围海养殖外部成本测算 ——以浙江舟山为例

我国海水养殖污染的一项重要原因在于对海水养殖生态环境影响的认识不清,作为我国仅次于开放式养殖的第二大海水养殖方式,围海养殖对海洋生态环境的破坏更大。为提高围海养殖效率,降低围海养殖的外部成本,提升公民海洋生态环境保护意识,文章以浙江省舟山市为例,基于391份随机调查样本,采用条件价值评估法对其围海养殖外部成本进行了评估。研究结果表明:舟山市围海养殖的外部成本约为42.97万元/月;相较于一般居民,围海养殖户对围海养殖外部成本的认知更为清晰和深刻;强化海洋生态环境保护宣传,改进围海养殖技术,出台更严格的围海养殖废水排放标准,将海域使用金征收标准与围海养殖废水处理水平挂钩等措施将有助于降低围海养殖的外部成本。     

海面水包油乳化液溢油量等效估算方法研究

激光诱导荧光(LIF)是一种主动光学探测技术,该技术已在海面未乳化溢油油种鉴别及油膜厚度评估方面取得了一定的研究成果,但乳化溢油的监测理论基础和探测方法尚未成熟。溢油乳化物的量化会直接影响溢油污染的应急处理和灾害评估,但目前对溢油乳化液溢油量的计算尚未有研究报道。本文以水包油乳化液为研究对象,基于等效思想,根据油量相等建立了水包油乳化液与油膜的荧光检测等效模型,并依据光的辐射传输机理推导出等效油膜厚度的估算公式,最后采用仿真实例对估算方法的适用性和有效性进行了验证。实验结果表明,当水包油乳化液含油率小于12×10^-6,厚度小于5.98 cm时,利用所提出的等效方法能有效计算出水包油实际溢油厚度。该方法可为海面乳化溢油量的估算提供一种创新性的方法,具有良好的借鉴意义和实际应用价值。     

2.5维海洋非线性共轭梯度反演应用效果研究

为了研究频率域2.5维非线性共轭梯度反演在海洋天然气水合物探测中的实际应用效果,利用美国Scripps研究所在South Hydrate海域采集的可控源电磁探测数据进行2.5维非线性共轭梯度反演计算,结合其他地球物理方法在同一海区的研究结果,综合分析2.5维反演计算的准确性和有效性。结果表明,非线性共轭梯度反演能够清晰地反映出该区域海底面以下1~2 km地层的电导率结构分布,其局部反演结果与地震反射、地震层序、测井分析结果基本一致,且对小尺寸异常、浅层高阻体分辨力较强,尤其在海洋天然气水合物探测方面具有较大潜力。     

一种UUV集群抵近目标侦测决策效率评估方法研究

为实现对不同 UUV 集群自主抵近目标侦测任务中决策效率的评估,分析了被测目标采取不同策略时 UUV 集群决策灵活度,提出了评估 UUV 集群决策效率的 5 个指标：任务准备时间、目标停车集群响应时间、目标转向集群响应时间、目标旋回集群响应时间和目标加速集群响应时间,并建立了评估指标体系。 参考已有试验测试结果,基于字典序法对 5 种类型的 UUV 集群决策效率进行综合评估,并分析了该方法的优缺点,为不同类型 UUV 集群决策效率评估提供了研究思路。     

Bathymetry Inversion with the Gravity-Geologic Method: A Study of Long-Wavelength Gravity Modeling Based on Adaptive Mesh

Modeling of long-wavelength gravity anomaly is crucial for bathymetry inversion with a gravity-geologic method. We propose a new method, named as iGGM, to approximate the long-wavelength gravity anomalies by using a finite element method based on an adaptive triangular mesh which is constructed by uneven control points. The mesh size is suitably controlled to ensure that there are several control points in each grid. By using iGGM, the bathymetry in the South China Sea (Test Area #1: 112°E–119°E, 12°N–20°N) and East China Sea (Test Area #2: 125°E–130°E, 25°N–30°N) is estimated. The performance of the method was evaluated by comparing the predictions with Earth topographical database 1 (ETOPO1) model and shipborne depths in the test points. Results show that the depths derived by iGGM have a strong correlation with the shipborne depths. In the test points, the mean values of their differences are smaller than 10 m. The standard deviations of their differences are smaller than 120 m and their correlation is stronger than 0.98. Meanwhile, the results provided by the iGGM model are comparable with that obtained by the ETOPO1 model, e.g., the differences between iGGM and ETOPO1 models in test points for Test Areas 1 and 2 are 116 and 70 m in standard deviation, respectively.     

Geoacoustic provinces of the northern South China Sea based on sound speed as predicted from sediment grain sizes

Marine Geophysical Research - The determination of geoacoustic provinces has important applications in assessing the responses of sedimentary environment and processes. Suitable geoacoustic...     

北部湾南部重力柱状样的MSCL地声学性质测量及分析

利用钻孔多参数连续记录仪(multi-sensorcorelogger,MSCL),对在南海北部湾南部海区所取得的6个孔重力柱状样进行了测量,获得了连续P波波速、湿密度和孔隙度数据。在室内实验室对柱状样分样之后的沉积物样品进行了孔隙度与湿密度的测定。利用不同的统计回归方法对所获得的6个柱状样的孔隙度与P波波速进行相关分析,并对比室内测量的P波波速、湿密度、孔隙度,建立了基于孔隙度数据对沉积物P波波速进行预测的方法,对MSCL测量方法的优缺点进行了讨论。结果表明,MSCL测量结果与实验室柱状样测定结果较为吻合。样品所在深度对声速的变化影响不大。孔隙度与声速的多项式回归,样条插值回归和GAM模型回归都获得了较高的相关系数,GAM模型可以提供一个较为接近测量值的声速预测方法。MSCL用来测量海底沉积物,可以获得大量的高密度、高精度的沉积物地声学及其他参数数据,但是,如果空气混入沉积物样品中则会导致MSCL测量结果失真。该研究为使用MSCL在区域海底地声学性质、恢复区域海洋沉积历史、海底地声学模型建立等研究提供了参考。     

Delaunay三角剖分的快速实现

主要探讨了以Delaunay三角剖分的逐点插入法为基础构建不规则三角网的方法,并在程序设计中对该算法进行了改进,极大地提高了Delaunay三角网的构建效率。