Matched-field inversion of sound speed profile in shallow water using a parallel genetic algorithm

A sound speed profile plays an important role in shallow water sound propagation. Concurrent with in-situ measurements, many inversion methods, such as matched-field inversion, have been put forward to invert the sound speed profile from acoustic signals. However, the time cost of matched-field inversion may be very high in replica field calculations. We studied the feasibility and robustness of an acoustic tomography scheme with matched-field processing in shallow water, and described the sound speed profile by empirical orthogonal functions. We analyzed the acoustic signals from a vertical line array in ASIAEX2001 in the East China Sea to invert sound speed profiles with estimated empirical orthogonal functions and a parallel genetic algorithm to speed up the inversion. The results show that the inverted sound speed profiles are in good agreement with conductivity-temperature-depth measurements. Moreover, a posteriori probability analysis is carried out to verify the inversion results.     

Geoacoustic inversion for sediments in the South China Sea based on a hybrid inversion scheme

We performed a long range acoustic propagation experiment in the South China Sea (SCS) in November 2004. The environment of the experiment was with an isothermal sound speed profile, where influence of water volume fluctuation was small, meaning that bottom parameters can be well estimated from acoustic signals. We inverted the acoustic parameters of sediment by using a hybrid inversion scheme that combines the matched field processing inversion with Hamilton sediment empirical relationship and transmission loss data. The numerical results show excellent agreement with the experiment data, indicating validity of the inverted parameters.     

Model of wind-generated ambient noise in stratified shallow water

In order to build a rapid ocean ambient noise model adapted for a stratified shallow water, a hybrid model of normal mode method (for far field) and ray method (for near field) is suggested which combines the advantages of both methods. Since the near field of wind-generated noise is not sensitive to the sound speed profile, the sound speed profile is regarded as a constant; which makes the model rapid and accurate. The simulation results are in agreement with those of the wave model.     

近地表弹性介质瑞雷波勘探研究进展与展望

瑞雷波法已在研究地球内部结构、近地表地球物理工程和超声无损检测等领域中获得了广泛应用,尤其是近年来瑞雷波法作为近地表场地表征新兴领域的前沿技术已成为国际学术研究与应用的热点。对国内外近地表弹性介质瑞雷波勘探的主要研究成果与进展进行了综述,通过对瑞雷波勘探现有的研究成果和进展密切追踪发现当前近地表瑞雷波勘探主要基于水平地表弹性水平层状介质模型,利用单分量瑞雷波相速度频散曲线单目标反演获得一维横波速度剖面和相关岩土力学参数。但是,该方法也存在着现有瑞雷波频散曲线反演极易出现模式误识别、现有单分量单目标瑞雷波反演未充分利用多分量信息、现有瑞雷波相速度反演未充分利用群速度传播特性、现有瑞雷波反演未充分利用质点椭圆极化振动特性等挑战性学术难题和不足。基于上述问题,建议未来近地表弹性介质瑞雷波勘探重点研究方向应集中在进行多模式表面波全速度谱反演研究、多站多分量表面波相速度多目标全速度谱反演研究、单站多分量表面波群速度多目标全速度谱反演研究和单站多分量表面波椭圆极化振动特性多目标反演研究。由此构建新的近地表多分量瑞雷波多目标全速度谱反演理论,引领多分量瑞雷波多目标反演学科前沿,拓展现有单分量瑞雷波单目标反演理论范畴,推动多分量瑞雷波高精度实用勘探技术的发展。      

THE EFFECTS OF A SLOPING BOTTOM ON THE EQUATORIAL WAVE-INDUCED DEEP RESIDUAL CURRENTS

ImODUcnONThedeepequatorialoceanhasobvioussignilicantflowapparentlycarryinghacelsa1ongandacrosstheequator(WissCtal.,l985).RmtfloatmsurementSshoWedthattheflowishigh1yvariable(Richardsonetal.,l993).ThomPsonandKawase(l993)pro-posedthatthelargeinstantaneousandfloatvelocitiesasWellasthevariabilitysuggestthatthetracersignaIsreflCCtreCtificationoftimedependentmotionsandniinginsteadofrneanEulerianflow.TheresultSoftheirstudyonthegenerationofmeancurrentSbyperi-odicfordngintheequatorialoceaninasir…     

Modal wave number spectrum for mesoscale eddies

The variations of ocean environmental parameters invariably result in variations of local modal wave numbers of a sound pressure field. The asymptotic Hankel transform with a short sliding window is applied to the complex sound pressure field in the water containing a mesoscale eddy to examine the variation of local modal wave numbers in such a range-dependent environment. The numerical simulation results show that modal wave number spectra obtained by this method can reflect the location and strength of a mesoscale eddy, therefore it can be used to monitor the strength and spatial scale of ocean mesoscale eddies.     

利用SH横波反射资料中的勒夫面波反演浅部横波速度结构

由于勒夫波速度与纵波速度无关,勒夫波多道分析方法具有频散曲线简单、能量清晰、反演过程中未知参数少、对初始模型依赖小、反演结果稳定、分辨率高等特点。实验结果表明,利用浅层SH横波反射地震数据中的勒夫面波信息,反演得到近地表横波速度结构,在浅层地质构造、断层位置、上断点埋深等方面与浅层SH横波反射地震剖面以及地质钻孔数据吻合较好。该方法可为地下构造解释提供更多依据,提高探测成果的可靠性。     

The relation of the average sound intensity to range in homogeneous shallow-water

Using a “three parameters” bottom model, the average sound intensity in homogeneous shallow water with constant depth is found by a better expression to be finite for normal mode number of water layer. The space structure of the sound field and some factors controlling the near, middle and far fields are discussed with this expression. This paper was published in Chinese inOceanologia et Limnologia Sinica 14 (2): 198–202, 1983.     

大兴安岭隆起北段最小一维速度模型反演和地震重新定位

使用2009-01-01～2017-03-30内蒙古地震台网和邻省台网记录的大兴安岭北段地区593个天然地震的3 848条P波到时数据,应用VELEST方法反演大兴安岭北段地壳一维P波速度模型、台站走时校正值和重新定位结果。联合反演获得的台站走时校正值反映了地壳速度的横向不均匀性及台站下方波速异常,大兴安岭北段东南侧速度高于西北侧;重新定位结果中,震中位置在空间上更加收敛;震源深度剖面图显示,大兴安岭北段西北侧震源深度较浅,东南侧较深,大兴安岭山脊区域震源深度分布下限较深。     

A New Multigrid 3D-VAR Optimization Method for Bottom Friction Using HF Radar Current Observation

This paper proposes a new method for data assimilation of the surface radial current observed by High Frequency ground wave radar and optimization of the bottom friction coefficient. In this method, the shallow water wave equation is introduced into the cost function of the multigrid three-dimensional variation data assimilation method as the weak constraint term, the surface current and the bottom friction coefficient are defined as the analytical variables, and the high spatiotemporal resolution surface radial flow observed by the high-frequency ground wave radar is used to optimize the surface current and bottom friction coefficient. This method can effectively consider the spatiotemporal correlation of radar data and extract multiscale information from surface radial flow data from long waves to short waves. Introducing the shallow water wave equation into the cost function as a weak constraint condition can adjust both the momentum and mass fields simultaneously to obtain more reasonable analysis information. The optimized bottom friction coefficient is introduced into the regional ocean numerical model to carry out numerical experiments. The test results show that the bottom friction coefficient obtained by this method can effectively improve the accuracy of the numerical simulation of sea surface height in the offshore area and reduce the simulation error.     

Numerical simulation and preliminary analysis of typhoon waves during three typhoons in the Yellow Sea and East China Sea

In this study, typhoon waves generated during three typhoons(Damrey(1210), Fung-wong(1416), and Chan-hom(1509)) in the Yellow Sea and East China Sea were simulated in a simulating waves nearshore(SWAN) model, and the wind forcing was constructed by combining reanalyzed wind data with a Holland typhoon wind model. Various parameters, such as the Holland fitting parameter(B) and the maximum wind radius(R), were investigated in sensitivity experiments in the Holland model that affect the wind field construction. Six different formulations were considered and the parameters determined by comparing the simulated wind results with in-situ wind measurements. The key factors affecting wave growth and dissipation processes from deep to shallow waters were studied, including wind input, whitecapping, and bottom friction. Comparison with in-situ wave measurements suggested that the KOMEN scheme(wind input exponential growth and whitecapping energy dissipation) and the JONSWAP scheme(dissipation of bottom friction) resulted in good reproduction of the significant wave height of typhoon waves. A preliminary analysis of the wave characteristics in terms of wind-sea and swell wave revealed that swell waves dominated with the distance of R to the eye of the typhoon, while wind-sea prevailed in the outer region up to six to eight times the R values despite a clear misalignment between wind and waves. The results support the hypothesis that nonlinear wave-wave interactions may play a key role in the formation of wave characteristics.     

Application of artificial bee colony algorithm in Rayleigh wave inversion

In order to solve the problems of multi-parameter,multi-extreme and multi-solution in the nonlinear iterative optimization process of Rayleigh wave inversion,the artificial bee colony(ABC)algorithm is selected for global nonlinear inversion.The global nonlinear inversion method does not rely on a strict initial model and does not need to calculate the derivative of the objective function.The ABC algorithm uses the local optimization behavior of each individual artificial bee to finally highlight the global optimal value in the colony,and the convergence speed is faster.While searching for the global optimal solution,an effective local search can also be performed to ensure the reliability of the inversion results.This paper uses the ABC algorithm to perform Rayleigh wave dispersion inversion on the actual seismic data to obtain a clear undergrounding of shear wave velocity profile and accurately identify the location of the high-velocity interlayer.It is verified that the ABC algorithm used in the inversion of the Rayleigh wave dispersion curve is stable and converges quickly.     

On Long-waves breaking II. in a half infinite channel with bottom slope

The propagation of long-waves, such as tidal waves from the coastal ocean into shallow estuarine waters, often produces asymmetries of velocity and water level in time series resulting in long-wave breaking. Tian (1994) studied the mechanism of long-wave breaking in an infinite channel with constant depth, consideringnth power bottom friction. This study is for the case of a half infinite channel with bottom slope, taking linear bottom friction into account. The wave breaking time and wave breaking location are estimated and the criteria for long-wave breaking in this particular case are obtained. The results obtained can also be easily applied to the case considering wind stress.     

基于卫星高光谱遥感影像的浅海水深反演方法

遥感水深反演具有非接触测量和省时省力等优点,能够为航海、岛礁工程与珊瑚礁生态调查等活动提供重要参考。随着高光谱遥感卫星数量的增长,基于高光谱遥感影像的水深反演具有良好的发展与应用潜力。HOPE（Hyperspectral Optimization Process Exemplar）算法是比较常用的高光谱水深反演算法。鉴于HOPE算法在低遥感反射率海域会出现水深被高估的问题,本文基于Hyperion高光谱遥感影像提出一种改进的水深反演算法。该算法针对危险或难以到达海域往往具有水体光学性质较为均一的特点,利用深水区遥感反射率的观测值来估计整个研究区域内的水体光学性质参数并将其固定,以便减少未知参数数量,解决水深被高估的问题,最终达到提高水深反演整体精度的目的。塞班岛和中业岛的实验结果表明,改进算法能够有效克服常规HOPE算法在低遥感反射率水域高估水深的问题。改进算法能够将平均遥感反射率小于0.0075sr^-1（塞班岛）和0.001 sr^-1（中业岛）范围内的水域的水深反演平均绝对误差从常规HOPE算法的2.94 m和6.44 m分别降低至2.56 m和4.99 m,从而能够相应地将整体的均方根误差从3.18 m和5.39 m分别降低至2.30 m和3.32 m,而将整体的平均相对误差从32.4%和27.1%分别降低至30.6%和23.9%。因此,改进算法在提高卫星高光谱遥感影像水深反演效果方面具有可行性和有效性。     

Study of the Sensitive Properties of Marine Gas Hydrate Based on the Prestack Elastic Inversion

Using a bottom simulating reflector(BSR) on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method. However, owing to the abundance differences between the gas hydrate and free gas in different regions, the BSR may be unremarkable on the seismic profile and invisible in certain cases. With the improvement of exploration precision, difficulty arises in meeting the requirements of distinguishing the abundance differences in the gas hydrate based on BSR. Hence, we studied other sensitive attributes to ascertain the existence of gas hydrate and its abundance variations, eventually improving the success rate of drilling and productivity. In this paper, we analyzed the contradiction between the seismic profile data and drilling sampling data from the Blake Ridge. We extracted different attributes and performed multi-parameter constraint analysis based on the prestack elastic wave impedance inversion. Then, we compared the analysis results with the drilling sampling data. Eventually, we determined five sensitive attributes that can better indicate the existence of gas hydrate and its abundance variations. This method overcomes the limitations of recognizing the gas hydrate methods based on BSR or single inversion attribute. Moreover, the conclusions can notably improve the identification accuracy of marine gas hydrate and provide excellent reference significance for the recognition of marine gas hydrate. Notably, the different geological features of reservoirs feature different sensitivities to the prestacking attributes when using the prestack elastic inversion in different areas.     

ON LONG-WAVES BREAKING Ⅱ.IN A HALF INFINITE CHANNEL WITH BOTTOM SLOPE

The propagation of long-waves, such as tidal waves from the coastal oceam into shallow estuarine waters,often produces asymmetries of veolcity and water level in time series resulting in long-wave breaking.Tian (1994) studied the mechanism of long-wave breaking in an infinite channel with constant depth,considering nth power bottom friction. This study is for the case of a half infinite channel with bottomslope, taking linear bottom friction into account. The wave breaking time and wave breaking location areestimated and the criteria for long-wave breaking in this particular case are obtained. The results obtainedcan also be e asily applied to the case considering wind stress.     

Dynamic process simulation with a Savage-Hutter type model for the intrusion of landslide into river

Natural damming of rivers by mass movements is a very common and potentially dangerous phenomena which has been documented all over the world. In this paper, a two-layer model of Savage-Hutter type is presented to simulate the dynamic procedure for the intrusion of landslide into rivers. The two-layer shallow water system is derived by depth averaging the incompressible Navier-Stokes equations with the hydrostatic assumption. A high order accuracy scheme based on the finite volume method is proposed to solve the presented model equations. Several numerical tests are performed to verify the realiability and feasibility of the proposed model. The numerical results indicate that the proposed method can be competent for simulating the dynamic process of landslide intrusion into the river. The interaction effect between both layers has a significant impact on the landslide movement, water fluctuation and wave propagation.     

瑞雷面波勘探的过去、现在和未来

回顾了瑞雷面波勘探发展的历史,剖析了面波频散曲线正反演的进展及存在的问题,重点指出在石油和煤炭反射地震中,应从剔除面波转变为利用面波,并借助现有石油“大炮”反射地震资料中“有害的”瑞雷面波提取浅层地质结构信息,为反射地震静校正、近地表工程地质评价等提供物性参数,把反射地震的面波干扰变为对勘探十分有利的参数提取源,实现以强干扰面波为主要数据源的浅层结构探测,以促进瑞雷面波勘探技术的进一步发展。     

A THIRD GENERATION SHALLOW WATER WAVE NUMERICAL MODEL-YE-WAM

This paper pnesents a third gneration shallow Whter disode spedtal wave nbotal medeIYE-WAM based on the spedtal action balance equation. The mode accounts for all edevan effectsof currents on waves, incuding tmpotally and spatialy varying depth and current inded refraction,sttalning and fequency shift and also explidtly takeS into aanunt all source terms, speclally adePth-limited breaking dheipation. In addition, an energy forcing scheme is propond and applied to themode's open boundaries to areUn for the propagution of sedIs into the study spstem The upwinddiffeIenng scheme and a standard hybrid diffdrencing scheme for the propagaion terrn and a simpleEuler method for the source teme are employed.     

基于归一化磁源强度的聚焦反演方法

针对剩磁条件下铁磁物质反演中存在的问题,提出基于归一化磁源强度的聚焦反演方法。首先,利用归一化磁源强度作为实测数据对磁性目标进行反演,减弱剩余磁化对反演结果的影响;然后,利用深度加权矩阵和最小支撑矩阵对经典Tikhonov正则化理论框架下的反演模型进行约束得到目标函数,并有效解决了核函数随深度增大而快速衰减的问题;最后,通过对目标函数进行迭代奇异值分解获得最佳物性参数,并根据Morozov偏差原则自适应地确定目标函数在迭代过程中的正则化参数,提高了迭代速度和求解精度。