时变重力场球面模型反演算法和模拟实验     

张贝  陈石  李红蕾  杨锦玲  韩建成  卢红艳 《地震》2021,41(1):13-24

时变重力场是研究地球内部介质物性变化的重要手段.本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法,该方法采用球坐标系下的六面体单元来模拟场源介质,适合大尺度地震流动重力测量数据的等效源模型构建.通过引入重力时变信号的一阶光滑先验条件,压制了时变重力信号中的短周期高频分量,可用于提取与地震孕育相关的长...  相似文献   

胶州湾增养殖海域营养状况与赤潮形成的初步研究   总被引：11，自引：1，他引：11  

郝建华  霍文毅  俞志明 《海洋科学》2000,24(4):37-41

根据对胶州湾女姑山增养殖海域1998年5月～9月的连续监测资料,参照潜在性富营养化的概念,应用NQI指数对该海域的营养状况进行分析。认为该海域水质富营养化是7月3～8日Skeletonemacostatum和Biddulphiaaurita混合型赤潮形成的基础,磷、硅营养盐的消耗是赤潮消亡的主要原因;赤潮消亡之后浮游植物群落发生演替,水体叶绿素a仍保持较高含量,NQI指数也相应较高,水质表现为磷限制潜在性富营养化,由于磷酸盐的限制没有发展为赤潮。  相似文献   

Numerical assessment of tsunami attack on a rubble mound breakwater using OpenFOAM®     

《Applied Ocean Research》2018

A numerical assessment study of tsunami attack on the rubble mound breakwater of Haydarpasa Port, located at the southern entrance of the Istanbul Bosphorus Strait in the Sea of Marmara, Turkey, is carried out in this study using a Volume-Averaged Reynolds-Averaged Navier-Stokes solver, IHFOAM, developed in OpenFOAM^® environment. The numerical model is calibrated with and validated against the data from solitary wave and tsunami overflow experiments representing tsunami attack. Furthermore, attack of a potential tsunami near Haydarpasa Port is simulated to investigate effects of a more realistic tsunami that cannot be generated in a wave flume with the present state of the art technology. Discussions on practical engineering applications of this type of numerical modeling studies are given focusing on pressure distributions around the crown-wall of the rubble mound breakwater, and the forces acting on the single stone located behind the crown-wall at the rear side of the breakwater. Numerical modeling of stability/failure mechanism of the overall cross-section is studied throughout the paper.The present study shows that hydrodynamics along the wave flume and over the breakwater can be simulated properly for both solitary wave and tsunami overflow experiments. Stability of the overall cross-section can only be simulated qualitatively for solitary wave cases; on the other hand, the effect of the time elapsed during tsunami overflow cannot be reflected in the simulations using the present numerical tool. However, the stability of the overall cross-section under tsunami overflow is assessed by evaluating forces acting on the rear side armor unit supporting the crown-wall of the rubble mound breakwater as a practical engineering application in the present paper. Furthermore, two non-dimensional parameters are derived to discuss the stability of this armor unit; and thus, the stability condition of the overall cross-section. Approximate threshold values for these non-dimensional parameters are presented comparing experimental and numerical results as a starting point for engineers in practice. Finally, investigations on the solitary wave and tsunami overflow experiments/simulations are extended to the potential tsunami simulation in the scope of both representation of a realistic tsunami in a wave flume and stability of the rubble mound breakwater.  相似文献   

小振幅海洋内波的演变、破碎和所致混合   总被引：1，自引：1，他引：0  

李丙瑞  范海梅  田纪伟  孙波  张占海 《海洋与湖沼》2010,41(6):807-815

利用基于谱方法和MPI并行运算的数值模式SpectralModel,直接数值模拟了三维小振幅海洋内波的演变、破碎和所致湍流混合,指出导致其不稳定而破碎的为PSI(parametric subharmonic instability)机制;对于内波破碎所致的湍流混合过程,分析了跨等密度面扩散系数kρ、混合效率γ、浮力通量谱、动能谱以及势能谱等统计性质:内波破碎前,kρ和γ保持低值水平,浮力通量谱值为负,且集中在低波数段;内波破碎后,kρ和γ迅速增大,最大值分别约为0.9×103m2/s和0.18,浮力通量谱值在低波数段为负值,在高波数段为正值,这是因为层化湍流中势能向小尺度运动传递和动能向小尺度运动传递相比更为有效。在内波破碎、强湍流混合阶段,势能谱存在一谱段满足kz3律,P(kz)=0.2N2kz3。此外,与二维模拟结果相比较,导致内波不稳定而破碎的均为PSI机制,kρ、浮力通量谱、势能谱变化趋势大体一致;但三维数值实验中,内波破碎时间提前,湍流衰减加快;KE谱在高波数部分下降速度相对减小,更接近于kz3。  相似文献   

伶仃洋海域潜在性富营养化评价与研究     

卢楚谦  蔡伟叙  余汉生 《海洋通报》2010,29(6)

根据2006年7月至2007年12月的调查资料,一年四季分析珠江口伶仃洋海域营养盐的空间分布特征及季节变化,该海域DIP表层年平均值为0.054mg/L,底层年平均值为0.53mg/L,DIN表层平均值为1.20mg/L,底层平均值为1.03mg/L;采用近海海域水体潜在性富营养化等级的划分标准和评价方法进行评价,评价结果表明,伶仃洋海域营养化水平基本上为磷限制潜在性富营养化至III级富营养水平。  相似文献   

Velocity potential formulations of highly accurate Boussinesq-type models     

Harry B. Bingham  Per A. Madsen  David R. Fuhrman 《Coastal Engineering》2009

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.  相似文献   

The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary   总被引：1，自引：1，他引：0  

Gábor Papp  Erika Szeghy  Judit Benedek 《Journal of Geodesy》2009,83(6):509-521

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.  相似文献   

位场曲化平的插值-迭代法   总被引：15，自引：5，他引：10       下载免费PDF全文

徐世浙  余海龙 《地球物理学报》2007,50(6):1811-1815

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.  相似文献   

逐步逼近曲化平方法研究   总被引：3，自引：1，他引：2       下载免费PDF全文

刘金兰  王万银  于长春 《地球物理学报》2007,50(5):1551-1557

位场资料的常规处理解释方法是建立在平面数据理论之上的,然而实际的位场数据大多为曲面数据,若把曲面数据当成平面数据进行处理,必然导致很大的误差,因此有必要进行“曲化平”处理. 曲化平目前存在的主要问题是计算精度低、计算量大,因此,研究快速、精度高且适合大数据量处理的曲化平方法具有重要的价值. 本文在已有的泰勒级数曲化平方法基础上提出了逐步逼近技术和平均平面技术,使得曲面位场资料处理的精度得到了很大提高. 最后通过理论模型和实际资料的处理验证了该方法的应用效果.  相似文献   

A ‘string function’ for describing the propagation of large-scale potential energy anomalies in a rotating fluid     

R. H. Tyler  R. Käse 《地球物理与天体物理流体动力学》2013,107(1-2):31-64

  相似文献