C-RESPONSE OF GEOMAGNETIC DEPTH SOUNDING ON A 1D THIN SHELL MODEL

This paper tries to formulate the C-response of geomagnetic depth sounding(GDS)on an Earth model with finite electrical conductivity. The computation is performed in a spherical coordinate system. The Earth is divided into a series of thin spherical shells. The source is approximated by a single spherical harmonic P₁⁰ due to the spatial structure of electrical currents in the magnetosphere. The whole solution space is separated into inner and external parts by the Earth surface. Omitting displacement current, the magnetic field in the external space obeys Laplacian equation, while in the inner part, due to the finite conductivity, the electromagnetic fields obey Helmholtz equation. To connect the magnetic fields in the inner and external space, the continuity condition of magnetic fields is used on the Earth surface. The external magnetic fields are expressed by the inner and external source coefficients, from which a new parameter called C-response is computed from the inner coefficient divided by the external coefficient, thus normalizing the actual source strength. The inner magnetic fields in each layer can be recursively derived by the continuity boundary condition of both normal and tangential components of the magnetic field from the initial boundary condition at core-mantle-boundary. The consistency of our C-responses with that from a typical 1-D global model validates the accuracy of the proposed algorithm. Numerical results also show that the C-response estimated from the geomagnetic transfer function method will deviate exceeding 5%from the actual response at longer periods than about 10⁶s, which means that ignoring the curvature of the Earth at extreme long periods will make inversion result unreliable. Therefore, an accurate C-response should be computed in order to lay a solid foundation for reliable inversion.     

基于多项式拟合均匀荷载面曲率曲线的结构损伤识别方法

为仅利用结构损伤状态的柔度矩阵对结构进行损伤程度识别,先对损伤状态的均匀荷载面曲率曲线进行最小二乘法拟合。根据曲率曲线差判断结构的损伤位置,对损伤位置的点进行剔除后,再利用未损伤位置上的点进行局部最小二乘法拟合,代替损伤前的均匀荷载面曲率曲线,用于结构的损伤定位与定量。通过一简支梁数值算例,先以理论的二次多项式进行拟合,考虑单损伤和多损伤的情况,进行损伤识别分析,再分析多项式次数、测点数目以及不同噪声水平对损伤定量精度的影响。结果表明:在一定范围内,次数越高拟合误差越小,但差别不明显,采用理论的二次多项式拟合即可满足结构损伤识别要求,无噪声的情况下,测点数目减半不影响损伤识别的精度,该方法具有一定的抗噪声能力。     

全球地磁感应测深数据三维反演

全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.     

二类水体表观光学特性的测量与分析--水面之上法方法研究

水体光谱特性的测量与分析是海洋水色遥感的基础工作之一。水体光谱特性包括两个方面：表观光学特性和固有光学特性。现场直接获取表观光学特性有剖面法和水面之上法。在多次现场试验的基础上,主要针对二类水体从误差来源的角度来对水面之上法进行方法研究。同时,借鉴一类水体水面之上法的多种数据处理方法,对二类水体水面之上法的数据处理方法进行探讨。     

大气散射对采用归一化植被指数进行赤潮遥感监测的影响研究

1　引言近年来赤潮灾害在我国海域频繁发生,对沿海环境、经济和人民生产生活都造成一定的影响,因此对赤潮进行实时监测并业务化运行就显得尤其重要.赤潮遥感监测是目前常用手段之一,在实践中取得了一定的效果,但到目前为止还没有一种成熟的算法能用于业务化监测赤潮的发生.很多研究人员借鉴归一化植被指数(NDVI)阈值法、海表温度(SST)异常等来判断赤潮的发生,但归一化植被指数阈值法经常会出现误判或无法区分正常与异常水体的情况.该阈值范围由先验知识得到,没有定量的理论依据.归一化植被指数是陆地遥感中用以表征地表植物覆盖程度和健…     

Seafloor topography and landform characteristics of Zengmu Shoal waters and their evolvement

The topography data of the Zengmu Shoal and its adjacent area are acquired by the SeaBat-8111 multibeam system.The full survey covering over the Zengmu Shoal area gives us new knowledge about the area. It is a flat underwater platform consisting ofcoral reefs, its fathom line is like a trap in near SN direction and its shape is like an isosceles triangle with the acute angle pointing to the north. The east and west sides ofthe shoal are steep while its south side is more gentle. As viewed from geological evolution, the Zengmu Shoal is a typical biologic coral landform caused by tectonic uplift.     

M-估计法广义变分同化FY-3B/IRAS通道亮温

采用FY-3B/IRAS亮温资料进行广义变分同化研究。广义变分同化结合了经典变分同化和稳健M-估计两者的优点。区别于经典变分同化依赖于先前的质量控制并要求误差服从高斯分布,把M-估计法耦合到经典变分同化框架中,得到广义变分同化,其弱化了同化前的质量控制和误差服从高斯分布这两个条件。目标能量泛函包含M-估计以保证对离群值具有稳健性,从而能够得到较好的同化结果。对比变分同化前后的FNL资料湿度与GDAS湿度相关系数作为同化结果检验评价。具体操作过程在FNL作为背景场的基础上分别采用经典和M-估计不同的权重因子变分同化FY3B/IRAS资料,把得到的分析场与GDAS进行相关性比较,由于湿度具有较强的非高斯性,文中首先评估了安徽省13个站GPS/PWV和积分相关湿度廓线得到大气可降水量(即GDAS/PWV和FNL/PWV资料)的相关性,进一步基于信息熵自由度思想进行了近一个月IRAS 20个通道对分析场的影响贡献率诊断研究。     

云贵川地区ZTD垂直剖面格网模型

利用ERA5大气再分析资料研究ZTD高程尺度因子的精细时间变化特征,构建顾及高程尺度因子精细时间变化的云贵川地区ZTD垂直剖面格网模型(YZTD-H模型)。以云贵川地区探空站分层ZTD数据作为参考值,检验YZTD-H模型的精度,并将其与GPT2w模型和GPT3模型进行比较。结果表明,顾及精细时间变化和垂直剖面变化的YZTD-H模型在时间维度和垂直剖面维度上均表现出较好的稳定性。     

液化场地桥梁群桩-土耦合体系强震反应分析

针对振动台试验,采用u-p形式控制方程表述饱和砂土的动力属性,选用土的多屈服面塑性本构模型刻画饱和砂土和黏土的力学特性,引入非线性梁-柱单元模拟桩,建立试验受控条件下液化场地群桩-土强震相互作用分析的三维有限元模型,并通过试验结果验证数值建模途径与模拟方法的正确性。以实际工程中常用的2×2群桩为例,建立桩-土-桥梁结构强震反应分析三维有限元模型。基于此,针对不同群桩基础配置对液化场地群桩-土强震相互作用影响展开具体分析。对比发现,桩的数量相同时,桩排列方向与地震波输入方向平行时比垂直时桩基受力减小5%~10%,而对场地液化情况无明显影响;相同排列形式下,三桩模型中土体出现液化的时间约比双桩模型延缓5s,桩上弯矩和剪力减小33%~38%。由此可见,桩基数量增加,桩-土体系整体刚度更大,场地抗液化性能显著,桩基对上部桥梁结构的承载性能明显增强,其安全性与可靠性更高。这对实际桥梁工程抗震设计具有一定的借鉴意义。     

青藏高原东南三江流域黄土类堆积体斜坡地震作用下的动力响应试验研究——以四川巴塘下归洼滑坡为例

青藏高原东南三江流域广泛存在堆积体斜坡,该类斜坡极易演化为滑坡灾害。以四川巴塘下归洼斜坡为例,通过振动台模型试验研究地震诱发的堆积体滑坡动力响应特征与破坏模式。试验结果表明,在地震作用初期,由于孔隙被挤密,斜坡的自振频率增大;在茂县波激励下,斜坡肩部动力响应最强烈。高幅值茂县波激励下,因斜坡整体刚度降低,放大效应减弱。Hibert谱描述了地震波在斜坡中的传播特征,向斜坡顶部传播时,地震波高频能量显著增强,地震波经过堆积层后,放大效应减弱。在足以致使堆积层滑动的强震作用下,坡表一定深度下动力响应呈现一致剧烈现象,斜坡上部会形成“脱离体”,“脱离体”在地震作用下上下颠簸运动;在临滑时,坡表出现最强烈的动力响应。斜坡的破坏模式为：地震初期,在重力和地震耦合作用下,坡表土体剥落;随地震幅值增加,坡肩与堆积层发生明显相对位移,堆积层表面出现裂缝;最终,坡脚发生明显相对位移,随后堆积层偏离基岩加速滑塌。基于地震信号特征的斜坡堆积层滑落过程可分为三阶段：(1)稳定阶段;(2)临滑阶段;(3)滑移阶段。