基于非结构化网格的任意复杂2D RMT有限元模拟

射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.     

日本东北M_W9.0地震海啸在港池及邻近区域诱发的涡流危险性计算与评估分析

海啸造成的灾害与损失并非都与淹没有关,特别是港口中海啸诱导的强流会对船只及海事设施产生重要的影响及损害.由于海啸流观测数据稀缺及海啸诱导涡流机制的不确定性,过去60年海啸科学主要集中于对海啸波特征及淹没过程的研究与分析,海啸流模拟及验证工作开展较少,导致对海啸流基本特征及其造成灾害现象的曲解.开展海啸诱导的涡流研究及预警服务显得尤为重要及紧迫.考虑快速海啸预警需要,综合对比海啸诱导涡流的物理框架及模型方法,探索兼顾效率与计算精度的海啸流模拟方法是本文的核心工作及出发点.通过分析浅层湍流相干结构(TCS)产生的主要物理耗散机制,确定了考虑2D水平耗散机制的非线性浅水方程可用于海啸涡流的模拟分析.基于高精、高分辨率有限体积模型Geoclaw建立了三个精细化的港口海啸流模型,模型分辨率为5m.利用基于海啸浮标反演的海啸源模型作为初始条件,模拟分析了日本东北地震海啸在远场的海啸波流特征.海啸波流特征模拟结果与观测吻合较好,结果可信.对比发现:波驱动的自由表面流,小的位相或波幅误差就会导致大的流速误差,流的模拟和预报相对波幅来说更具挑战性.研究了海啸波流能量在港池中的分布特征,得到:港池入口及防波堤两端常被强流控制,具有极高的危险性;相对于波幅的空间变化,海啸流具有更强的空间敏感性;所建立的高分辨率海啸模型模拟再现了日本海啸在近场的涡旋结构,给出了与观测基本一致的涡流特征.最后,引入海啸流危险等级标准,分析了港口海啸流危险性等级分布、船只疏散的安全深度及回港的时间周期.针对港口、海湾同时考虑海啸波流特征的海啸预警与评估对于港口应急管理者科学决策具有重要意义.     

西沙地块地壳结构及其构造属性

西沙地块作为在南海形成演化过程中形成的微陆块,记录了南海演化历史的重要信息,其地壳结构、物质组成及构造属性是探讨南海形成演化的关键.基于采集到的OBS2013-3测线海底地震仪数据,用射线追踪和正演走时拟合方法,获得了西沙地块的二维纵波速度模型.模型显示沉积层速度为2.2~3.2km·s-1,厚度为0.8~3.0km,局部基底面起伏较大,上地壳顶部速度为5.0~5.5km·s-1,下地壳底部速度为6.9km·s-1,上地幔顶部速度为8.0km·s-1.西沙地块的地壳厚度平均为23km,上地壳厚度约为9km,下地壳厚度约为14km,莫霍面埋深为23~27km.从穿过西沙地块的纵、横两条大剖面推算,块体大小约为9.2×105 km3,与华南陆缘相比,表现为整体减薄的陆壳特征.西沙地块与南沙地块垂直于西南次海盆扩张脊分布,根据二者地壳结构的特征对比,二者互为共轭关系.     

OBS广角地震探测中二次反射Pg震相特征及在地壳结构成像中的作用

在海底地震仪(OBS,Ocean Bottom Seismometer)广角地震探测中,经常可以见到能量很强、连续性很好的多次波震相,但对于如何确定这些多次波的属性以及怎样充分利用它们来约束地下的结构特征,至今还缺乏详细的研究.本文利用南海北部已获取的OBS探测数据,对与地壳折射Pg震相近乎平行并紧随其后的二次反射Pg震相进行了系统分析研究,发现在实测OBS地震记录剖面上,二次反射Pg震相具有连续、清晰、可追踪偏移距较远(约60km以上)等特点,其地震波形和粒子运动轨迹与初至Pg震相相似,但二次反射Pg震相波形最大振幅值和粒子震动能量明显比初至Pg震相大.通过理论模型模拟以及对实测地震剖面三种不同反射层路径的走时计算,确定了二次反射Pg震相主要来自沉积层的反射.在此基础上,对二次反射Pg震相在地壳结构成像中的作用进行分析,发现加入二次反射Pg震相前后,由于反射震相的增加,沉积层界面的约束程度得到极大的提高;另外,通过对理论模型和实测剖面OBS2010地壳结构进行加入二次反射Pg震相前后的分辨率测试,结果发现加入二次反射Pg震相数据后,由于射线密度的增加,沉积层和上地壳结构的成像分辨率有显著改善.     

Features and formation processes of multiple deposition layers from the 2004 Indian Ocean Tsunami at Ban Nam Kem,southern Thailand

Multiple‐layered tsunami deposits have been frequently reported from coastal stratigraphic sequences, but the formation processes of these layers remain uncertain. A terrestrial sandy deposit formed by the 2004 Indian Ocean Tsunami was investigated at Ban Nam Kem, southern Thailand. Four internal layers induced by two tsunami waves were identified in the tsunami deposit. Sedimentary structures indicated that two units were formed by run‐up currents caused by the tsunami and the other two units were deposited by the backwash flows. Graded bedding was common in the layers, but inverse grading was observed at limited intervals on the surveyed transects. The characteristics of the multiple‐layered tsunami deposit vary remarkably over a very short distance (<1 m) in response to the local topography. Remarkable asymmetries in thickness and grain‐size distribution are recognized between the run‐up and backwash flow deposits. On the basis of the interpretation of sedimentary structures, the formation process of the multiple‐layered tsunami deposit observed in this study can be explained in a schematic model as the modification of the ideal tsunami sequence by local erosion and the asymmetric hydraulic properties of tsunami waves, such as the maximum shear velocity and the heterogeneity of the flow velocity field.     

Morphodynamique d'une dune sous-marine du détroit du pas de CalaisMorphodynamics of an undersea sandwave of the Dover Straits

This paper investigates the dynamics of sandwaves in relationship with tidal currents and weather conditions. The studied sandwaves are located in the Dover Straits and are covered by megaripples. These megaripples have migration speeds of $\text{1 }\text{m}\text{h}{}^{\mathrm{?1}}$ and slopes of 34°, suggesting the potential for avalanches to occur along the flanks of the sandwave. Tidal cycles without reversing currents were observed during stormy weather. Wind-induced currents lead to a unidirectional migration of megaripples and sandwaves. Well-defined areas without megaripples were observed and correlated with an increase in sandwave height. We propose interpretations for understanding sandwave saturation and migration. To cite this article: D. Idier et al., C. R. Geoscience 334 (2002) 1079–1085.     

底拖网配纲设计计算

根据疏目拖网逐目和逐个编结周期或剪裁循环地装配的要求,提出了先计算口门每目和网翼或网缘的每编结周期或剪裁循环的配纲长度,然后再计算出整段口门和网翼或网缘的配纲长度之设计计算方法。     

Event-stratigraphy of a caldera-forming ignimbrite eruption on Tenerife: the 273 ka Poris Formation

The 273 ka Poris Formation in the Bandas del Sur Group records a complex, compositionally zoned explosive eruption at Las Cañadas caldera on Tenerife, Canary Islands. The eruption produced widespread pyroclastic density currents that devastated much of the SE of Tenerife, and deposited one of the most extensive ignimbrite sheets on the island. The sheet reaches ~ 40-m thick, and includes Plinian pumice fall layers, massive and diffuse-stratified pumiceous ignimbrite, widespread lithic breccias, and co-ignimbrite ashfall deposits. Several facies are fines-rich, and contain ash pellets and accretionary lapilli. Eight brief eruptive phases are represented within its lithostratigraphy. Phase 1 comprised a fluctuating Plinian eruption, in which column height increased and then stabilized with time and dispersed tephra over much of the southeastern part of the island. Phase 2 emplaced three geographically restricted ignimbrite flow-units and associated extensive thin co-ignimbrite ashfall layers, which contain abundant accretionary lapilli from moist co-ignimbrite ash plumes. A brief Plinian phase (Phase 3), again dispersing pumice lapilli over southeastern Tenerife, marked the onset of a large sustained pyroclastic density current (Phase 4), which then waxed (Phase 5), covering increasingly larger areas of the island, as vents widened and/or migrated along opening caldera faults. The climax of the Poris eruption (Phase 6) was marked by widespread emplacement of coarse lithic breccias, thought to record caldera subsidence. This is inferred to have disturbed the magma chamber, causing mingling and eruption of tephriphonolite magma, and it changed the proximal topography diverting the pyroclastic density current(s) down the Güimar valley (Phase 7). Phase 8 involved post-eruption erosion and sedimentary reworking, accompanied by minor down-slope sliding of ignimbrite. This was followed by slope stabilization and pedogenesis. The fines-rich lithofacies with abundant ash pellets and accretionary lapilli record agglomeration of ash in moist ash plumes. They resemble phreatomagmatic deposits, but a phreatomagmatic origin is difficult to establish because shards are of bubble-wall type, and the moisture may have arisen by condensation within ascending thermal co-ignimbrite ash plumes that contained atmospheric moisture enhanced by that derived from the evaporation of seawater where the hot pyroclastic currents crossed the coast. Ash pellets formed in co-ignimbrite ash-clouds and then fell through turbulent pyroclastic density currents where they accreted rims and evolved into accretionary lapilli.Editorial Responsibility: J. Stix