GRACE和SLR观测的地球动力学扁率最大熵谱及小波相关分析

卫星重力测量技术的实现为测定地球动力学扁率提供了新的方式和途径,GRACE卫星是目前最新的重力测量卫星,据其恢复的低阶重力场较以往精度得到大大提高,然而其观测地球动力学扁率（二阶项）却与卫星激光测距（SLR）结果相差较大.本文采用最大熵谱和小波分析方法对GRACE和SLR观测的地球动力学扁率时间序列信号进行定量比较分析,结果表明：GRACE观测的地球动力学扁率年际周期变化振幅仅为SLR观测结果的25%,并且目前GRACE观测的地球动力学扁率数据中含有系统输入信息和相位差,但前者较后者包含有较强的短周期（2～6月）信息.造成这种差异的主要原因可能来自于GRACE与SLR全球观测数据时空分布不同.     

叠前地震反演技术的进展及其在岩性油气藏勘探中的应用

与叠后地震反演相比,叠前地震反演可得到更丰富的储层信息,可提高储层的描述精度.目前,叠前地震反演主要包括弹性阻抗反演、叠前P波阻抗和S波阻抗联合反演、叠前地震波形反演.文中概述了叠前地震反演各项技术的主要进展及其典型应用实例,提出了叠前地震反演技术在岩性油气藏勘探开发中的应用策略.叠前地震反演技术在未来的油气勘探开发中将有广阔的应用前景.     

黄海典型漂浮大型藻类光谱特征分析

利用实测高光谱和多光谱相机数据对在黄海海域采集的浒苔和马尾藻样本进行光谱特征分析。结果显示:漂浮大型藻类的高光谱反射率均值(Rmean_H)与多光谱反射率均值(Rmean_M)、高光谱藻类指数(VB-FAHH)与多光谱藻类指数(VB-FAHM)之间的相对误差小于15.6%,说明这两种数据具有较好的一致性;基于光谱特征波长的反射率建立比值光谱指数,其中浒苔和马尾藻的差异最为明显的是R640/R554指数;基于多光谱相机计算反射谷深度指数(T-depth),采用阈值法来区分浒苔和马尾藻,T-depth阈值的区分精度分别为69.6%±30.0%和69.6%±24.3%,稍优于虚拟基线漂浮藻类指数(VB-FAH)阈值的区分精度。     

长江口南汇湿地植被的光谱吸收特征研究

植被的光谱吸收特征与植被生长状况密切相关,同时也受土壤湿度等外界因素影响。本文以长江口南汇湿地典型植被芦苇、互花米草和海三棱藨草为研究对象,在计算其光谱特征参数的基础上,分析了这3种植被的光谱吸收特征形变PSDI值。最后结合实测的土壤湿度,分析了3种植被的PSDI值和土壤湿度的相关性。研究结果发现:互花米草的PSDI值最大,海三棱藨草的PSDI值最小,芦苇介于两者之间,表明互花米草受外界的环境影响较大;互花米草的PSDI值与土壤湿度的相关性最大,芦苇的PSDI值与土壤湿度的相关性最小,海三棱藨草介于两者之间,表明互花米草适应湿地环境能力较强。     

基于并行计算的海洋平台谱分析疲劳筛选技术研究

为了得到海洋平台疲劳失效风险最大的焊接区域,提出了以有限元网格单元应力为基础的谱分析筛选方法,并开发了程序。程序采用并行计算架构,单次完成数万个单元的累积损伤度/疲劳寿命计算,并通过有限元软件以云图方式显示疲劳寿命。采用该程序计算双浮筒半潜式平台的疲劳寿命。结果表明,该型平台疲劳失效风险较高的区域为立柱与上壳体下浮体连接对角线转角区域、浮筒中纵舱壁与立柱连接区域。并行计算方法的引入,大大缩减疲劳损伤度计算的分析时间,采用22核心的CPU,时间缩短为原来的1/18.5,17小时完成全平台外壳单元谱分析计算。     

基于遥感光谱反射率反演黄河口海域Ⅱ类水体悬浮泥沙浓度模型

依据黄河口海域含沙水体野外实测光谱反射率为参考实验数据,结合水体表面光谱观测数据建立了经验回归模型,从而定量估测黄河口海域水体表层悬浮泥沙质量浓度,最终得出黄河口海域不同年份水体的悬浮泥沙质量浓度图,这对分析河口海岸带冲淤变化过程具有重要的意义。     

基于南海北部海区浮游植物吸收光谱斜率变化的粒级结构反演

浮游植物的粒级结构是一个重要的生物参数。基于南海北部海区不同水体环境下测量的生物光学数据, 作者深入研究了粒级结构对浮游植物吸收光谱的影响。结果表明, 选择443和510nm波段计算得到的浮游植物光谱斜率S对粒级结构的变化具有较高的敏感性, 其随着小型浮游植物比例的增大呈不断增加的趋势。S与水体叶绿素a浓度、浮游植物吸收系数(aph(443))之间表现出明显的正相关特征。以40%为界对不同粒级浮游植物的优势进行定义, 发现在S与叶绿素a浓度、aph(443)的关系分布中小型(Micro)和微微型(Pico)浮游植物占据优势的水体表现出较为明显的分界, 叶绿素a浓度和aph(443)分别在0.70mg•m-3和0.05m-1附近, 相应的S在0.0004(m•nm)-1左右。基于实测数据建立的遥感反射率蓝绿波段比值与S之间的统计关系, 决定系数高达0.91, 为从水色遥感数据反演浮游植物粒级结构提供了重要手段。     

用H/V谱比法计算云南区域数字地震台站的场地响应

用H/V谱比法,计算给出了云南区域数字地震台网23个子台S波随频率变化的台站场地响应特征。结果表明:在1～10Hz频段内,23个子台S波的场地响应相对较平坦,在1.41～2.91之间变化,平均2.08;在大于10Hz的高频段,部分台站的场地响应有较明显的放大。     

Attenuation of compressional waves in peridotite measured as a function of temperature at 200 MPa

A technique has been developed to determine attenuation in rocks at high temperature using a gas-media, high-pressure apparatus. A pulse transmission technique and a spectral ratio method are used to study compressional seismic properties of rocks. Seismic waves are transmitted to and from the sample through buffer rods of mullite. The effect of seismic wave reflections within the sample assembly are cancelled out by taking ratios of the spectra measured at different temperatures. In order to obtain good signal-to-noise ratio for resolving the attenuation at high pressure and temperature, special care is taken in the sample assembly and the ultrasonic coupling between the sample, buffer rods and transducers. A very tight connection of the sample-buffer rod-transducer is essential for obtaining high frequency signals (>300 kHz) at high temperature. A small mass is attached to each outside end of the transducer to drive low frequency signals (<250 kHz) into the sample. Before attenuation measurements, the sample and the buffer rods are tightly compacted in a platinum tube at high pressure and room temperature to ensure pressure seal of the sample assembly. The frequency range of measurement covers 50 to 450 kHz for the sample. Attenuation is very small in the buffer rod compared to the sample for the entire temperature range of the study. Because of the small attenuation, a wide frequency band of 50 kHz to 3.2 MHz can be covered for investigating the attenuation in the buffer rod. The technique has been used to measure attenuation at high confining pressure, and temperatures including sub- and hyper-solidus of upper mantle rocks. Therefore, effects of partial melting on attenuation can be studied.The method is applied to the attenuation measurement in a peridotite as a function of temperature to 1225°C at 200 MPa confining pressure. At high temperature, signal amplitude decays more rapidly at high frequency than at low frequency, from which attenuation (andQ) can be determined using a spectral ratio method. No frequency dependence ofQ is resolved for both the sample and the buffer rod over the entire temperature and frequency ranges of the measurement. The results show thatQ decreases rapidly with increasing temperature even in the temperature range below the solidus of peridotites. Such temperature sensitivity ofQ is probably more useful to probe thermal structure in the upper mantle than that of conductivity at temperatures below the solidus. The results in this study are compared with available seismic velocity, electrical conductivity and solidus data for peridotites, suggesting that there is no discontinuous change in both mechanical and electrical properties of peridotites at the solidus temperature. Even at hypersolidus temperatures, it appears that velocity drops and conductivity increases continuously (not abruptly) with increasing melt fraction. This implies that mechanical and electrical properties of the upper mantle will gradually change at the boundary where the geotherm crosses the solidus.     

Global spectra of energy and enstrophy and their fluxes during July 1979

Transient and stationary spectra of kinetic energy (KE), available potential energy (APE) and enstrophy (EN), and their spectral fluxes as a function of the two-dimensional wavenumbern were computed for July 1979. Triangular truncation at zonal wavenumber 42 was used for computation. The slopes of various spectra in the wavenumber range 14≤n≤25 were obtained by fitting a straight line in log-log scale by the least square method. The transientKE, APE andEN spectra in the lower (upper) troposphere had slopes −2·21 (−2·30), −2·65 (−2·64) and −0·36 (−0·46), respectively. The effect of stationary and divergent motion on the slope values was investigated. The possible correlation between the slope and percentage of transient component in the combined energy and enstrophy was examined to identify the transient motion of the atmosphere with the two-dimensional homogeneous isotropic turbulence. The vertically averaged slope of kinetic energy and enstrophy in the lower (upper) troposphere was close to the value at 700 (200) hPa level. The spectral fluxes of kinetic energy and enstrophy in the wavenumber range 14≤n≤25 satisfied, to a very rough approximation, the criteria of inertial subrange. The stationary fluxes were small. The estimated stationary-transient component of flux was larger, comparable and less than the corresponding transient flux of APE, KE and EN. Representative levels for computation of energy and enstrophy spectra and their fluxes in the lower and upper troposphere were identified.