基于压缩感知的高分辨率平面波分解方法研究

平面波分解方法是地震资料处理中的一项关键技术,它广泛地运用在平面波偏移、各类Beam偏移等成像方法中.平面波分解不仅可以提高偏移成像的效率,而且可以压制地震资料中的随机噪音,提高地震数据的信噪比.线性Radon变换（LRT）是一种常见的实现平面波分解的方法,但常规的LRT存在以下两个缺点：（1）分辨率受测不准原理限制;（2）变换结果存在很多噪音和空间假频.为了克服LRT的上述缺点,本文提出一种基于压缩感知的高分辨率平面波分解方法,并利用加权匹配追踪（WMP）技术实现了该方法.该方法将LRT视为一个参数估计问题,并将LRT结果的稀疏性作为约束条件,在压缩感知理论的指导下利用WMP方法得到高分辨率、高信噪比的平面波分解结果.另外,该方法还可以用于提取地震数据的线性信号、压制随机噪音、实现高维地震数据规则化等地震资料处理技术.数值实验结果证明：WMP方法可以有效地提取地震数据中的线性信号,提高LRT的分辨率和信噪比,从而改善平面波分解的质量.     

亚洲地区ECMWF/NCEP资料计算 ZTD的精度分析

对流层延迟是卫星导航定位的主要误差源,气象观测的数值预报资料可用来计算对流层延迟改正量.本文通过分布于亚洲地区的49个GPS台站一年的实测ZTD资料,对利用欧洲中尺度天气预报中心(ECMWF)分析资料、美国国家环境预报中心(NCEP)再分析资料和NCEP预报资料,计算对流层天顶延迟(ZTD)改正的有效性和可能达到的精度进行了评估,分析了ECMWF和NCEP在亚洲地区的适用程度和其分辨率对计算ZTD精度的影响.研究结果表明:(1)相对于 GPS实测ZTD,用ECMWF资料计算ZTD的bias和rms分别为-1.0 cm 和2.7 cm,优于NCEP再分析资料,可用于高精度ZTD研究和应用;NCEP预报数据计算ZTD的bias和rms分别为2.4 cm 和 6.8 cm,足以满足广大GNSS实时导航定位用户对流层延迟改正的需要.(2)bias和rms呈现明显的季节性变化,总体上夏季大,冬季小;在空间分布上随着纬度的变化不明显,但随高度的增加rms总体上有递减趋势.另外还发现,亚洲东部地区夏季日平均bias和rms和南部热带地区冬季的日平均bias和rms变化相对较大.(3)ECMWF2.5°和0.5°的资料进行了对比分析,发现0.5°分辨率资料的rms比2.5°减小1~5 mm.这些结果,为在亚洲地区的空间大地测量、导航定位和INSAR等工作中,应用ECMWF/NCEP的资料进行对流层大气延迟改正的有效性和可能达到的精度提供了重要参考.     

基于Maxwell型阻尼器的多重调谐质量阻尼器性能评价

研究了基于Maxwell型阻尼器的多重调谐质量阻尼器(MTD—MTMD)在控制结构地震反应方面的最优动力特性。利用建立的设置MTD-MTMD时结构的传递函数,定义了设置MTD—MTMD时结构的动力放大系数(DMF)。将MTD-MTMD的优化准则定义为结构最大动力放大系数的最小值的最小化(Min．Min．Max．DMF)。利用定义的优化准则,评价了Maxwell型阻尼器的松弛时间系数(RTC)对MTD—MTMD最优参数和有效性的影响。利用最大的MTD—MTMD动力放大系数(DMF),评价了RTC对MTD-MTMD冲程的影响。     

基于CMPCC二维面波的吸收散射综合分析法

多道瞬态面波方法(MASW)在铁路路基质量及岩溶注浆检测中有许多应用,但分辨率不高限制了MASW更广泛的应用.针对MASW存在的问题,共中心点互相关(CMPCC)二维面波方法显著地提高了横向分辨率.但在复杂地下结构的情况下,CMPCC通过对频散曲线一维反演重构获得的二维横波断面,有可能畸变或损失大量地下信息,造成分析解释错误.通过大量数据对比分析研究,首次定义了表征面波传播的动力学参数——吸收散射系数,提出了基于CMPCC二维面波的吸收散射综合分析法(ASCA):(1)将表征面波传播动力学性质的吸收散射系数(AS)和运动学性质的相速度c一并综合分析,提高了空间分辨率及对异常性质解释的准确性;(2)继承了CMPCC方法优点,放弃了对频散曲线进行一维反演,数据处理没有人为影响,图件简单直接,易于分析;(3)在地基沉降与塌陷、地基处理与注浆检测、黄土湿陷、活动断层位置确定等浅层勘察中有着广泛的应用前景.     

考虑卷积完全匹配层数值色散的井间电磁三维正演

本文将以卷积完全匹配层为吸收边界条件的时域有限差分法应用到井间电磁的三维正演模拟中.证明了卷积完全匹配层中的数值色散会因为有效延伸因子而产生,列举常规有效延伸因子和网格间距对电磁波相速度各向异性的影响,并通过波场快照验证卷积完全匹配层中数值色散的存在;进而推广三维卷积完全匹配层中最大有效延伸因子、最大网格间距与激励源主频之间的约束,以此完善卷积完全匹配层的最优参数设置方案.在此基础上,通过正演结果中二次场的垂直磁场分量和走时来展示灵敏度分布,以此划定井间电磁勘探的优势区域,提出井间电磁正反演所需要的约束测井系列和最佳取井方案.     

三种河流相对高程模型生成方法对比

相对高程模型（REM,Relative Elevation Model的缩写）表示相对于河流水面或活性河道的海拔高度,衍生于数字高程模型（DEM,Digital Elevation Model的缩写）,能够消除因河流纵向的地势造成的DEM表征不清楚,高差过大造成细节模糊地势的影响,从而在平面图上更为清晰地表征河流地形的细微变化,因而研究REM的生成方法及其适用性具有十分重要的意义。通过文献综述,选取美国加利福尼亚州内华达山脉的卡森河某河段为实验区,利用其高精度（米级分辨率）DEM数据,介绍核密度法、反距离加权法、横截面插值法三种REM生成方法。通过REM和DEM的对比,以及三种REM的效果对比,得出：当高程差过大时,DEM可视化方案往往不能很好地在河流地貌刻画中发挥作用,尤其是在解译时,利用REM可解决高差过大造成的细节模糊问题;横截面插值法主要的优势表现在可根据用户需要做适用性调整,相对而言鲁棒性强,较为灵活,但其自动化程度较低,且花费的时间较长,因此需较多的人工干预;核密度法耗时相对较少,REM结果很少有伪影或错误的相对高程值,但核密度法导致一些潜在的伪影或者错误的相对高度,体现在...     

页岩气储层地球物理测井研究进展

页岩气是储存在页岩中以游离态和吸附态存在的天然气,是一种非常规油气资源.地球物理测井技术具有方法多、纵向分辨率高、信息量大、连续、方便、灵敏等优势,能够获得沿井眼剖面的岩石物理参数,是常规油气资源和非常规油气储层评价的重要手段.文章综述了页岩气储层的地质特点,总结了国内外页岩气测井系列和页岩气的测井响应特征.针对测井信息的处理,给出了基于测井有机碳含量(TOC)计算和成熟度指数(MI)的研究进展,讨论了页岩岩性识别、矿物计算和裂缝的评价方法,并对美国页岩气测井评价实例做了介绍.最后,指出目前页岩气测井技术存在的问题和未来的突破方向和思路,对指导页岩气勘探和经济评价具有重要指导意义.     

土木工程结构鲁棒控制的发展

评述了结构控制的发展,指出发展结构鲁棒控制策略的重要性。重点评述了结构双重调谐质量阻尼器(DTMD)和多重双重调谐质量阻尼器(MDTMD)的控制策略,提出了需进一步发展主动双重调谐质量阻尼器(ADTMD)和主动多重双重调谐质量阻尼器(AMDTMD)控制策略、此外,评述了结构鲁棒控制的设计准则与高层建筑和大跨桥梁在风与地震作用下的统一自适应主动鲁棒控制策略。     

半周期半对称延拓的EMD边界处理方法的实现

环境脉动下结构非平稳输出信号处理和复杂结构模态参数的识别方法中,希尔伯特-黄变换是种有效的方法.该方法的核心是经验模态分解和希尔伯特-黄变换谱分析.针对何立志(2006)提出的半周期半对称延拓经验模态分解边界的处理方法,本文以EL Centro地震波的筛分为例,对该方法的实现进行了详细说明和具体运用,证明该方法具有一定的实用价值.     

基于SPT-APD-DDA的砂土液化评价方法研究

传统的规范法在建立经验判别准则过程中存在大量不确定性因素,工程中亟需针对复杂场地液化判别方法的综合性研究。通过对上海市包头路段沉管抢修项目49组砂土土样的研究,在标准贯入试验法(SPT)定性评价的基础上采用提出的绝对差值百分比法(APD),划分出35组判别结果精确及14组存在判别误差的土样。借鉴判别分析理论思想,将平均粒径D₅₀、不均匀系数C_u、比贯入阻力P_s、标准贯入点深度d_s、地下水埋深d_w和标准贯入锤击数N_63.5等物理力学参数加入液化评价,利用精确土样建立适宜研究区工程地质特性的距离判别分析模型(DDA),对14组存在判别误差或失误的土层进行定量评价。研究结果表明：建立的砂土液化综合评价体系对研究区具有较强的针对性与适应性,从多角度全面、客观地评价了研究区砂土液化情况。最终通过合理的抗液化沉陷措施,保障了后续工程的稳步进行,为类似地质条件场区的砂土液化稳定性评价提供参考和借鉴。     

Contribution to the knowledge of the magmatic evolution by the study of the variation of the coefficient of viscosity

The determination of the coefficient of viscosity of eruptive products gives useful elements to the knowledge of possible variations of composition and physical conditions of the magmas.     

Pictorial Series of the Manifestations of the Dynamics of the Earth

Studia Geophysica et Geodaetica -     

Evolution of the Source of the Active Region of the Sun

Geomagnetism and Aeronomy - The magnetic evolution of the active region of the Sun at the level of the photosphere has been studied. Magnetic synoptic maps of the 2007 Kitt Peak Observatory were...     

Magnetic stratigraphy of the Ordovician in the lower reach of the Kotuy River: the age of the Bysy-Yuryakh stratum and the rate of geomagnetic reversals on the eve of the superchron

Until recently, the existing data prevented the geophysicists from accurately dating the Bysy-Yuryakh stratum, which outcrops in the middle reach of the Kotuy River, constraining the time of its formation to a wide interval from the end of the Late Cambrian to the beginning of the Silurian. The obtained paleomagnetic data unambiguously correlate the Bysy-Yuryakh stratum to the Nyaian regional stage and constrain its formation, at least a considerable part of it, by the Tremadocian. This result perfectly agrees with the data on the Bysy-Yuryakh conodonts studied in this work and yields a spectacular example of the successful application of paleomagnetic studies in solving important tasks of stratigraphy and, correspondingly, petroleum geology. Within the Bysy-Yuryakh stratum, we revealed a large normal-polarity interval corresponding to the long (>1 Ma) period when the geomagnetic reversals were absent. This result, in combination with the data for the Tremadocian and Middle–Upper Cambrian sequences of the other regions, indicates that (1) the rate of occurrence of the geomagnetic reversals on the eve of the Ordovician Moyero superchron of reversed polarity was at most one reversal per Ma; (2) the superchron does not switch on instantaneously but is preceded by a certain gradual change in the operation conditions of the dynamo mechanism which, inter alia, manifests itself by the reduction of the frequency of geomagnetic reversals with the approach of the superchron. This finding supports the views according to which a process preparing the establishment of the superchrons takes place at the core–mantle boundary.     

The effect of the pollution of the ground layer of the atmosphere on the albedo

Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.     

On the solution of the reverse problem of geodesics with the help of the chord

Резюме Дается вывод формул для вычислений длины хорды, длины и азимутов нормаляных сечений и длины геодезической линии эллипсоида. Применяются тригонометрические функции данных величин и постоянные принятого ргфгргнц-эллипсоида, что является целесообразным при использовании вычислителQjных машин, особенно автоматов. Рабочие формулы рекомендуются в форме (8–11).

---

---

Dedicated to Professor František Fiala on His 85th Birthday

---

Address: Veveři 95, Brno.     

On the maintenance of the axisymmetric part of the flow in the atmosphere

Summary The maintenance of the axisymmetric component of the flow in the atmosphere is investigated by means of a steady-state, quasi-geostrophic formulation of the meteorological equations. It is shown that the meridional variations in the time-averaged axisymmetric variables can be expressed as the sum of three contributions, one being due to the eddy heat transport, another to the eddy momentum transport, and a third to the convective-radiative equilibrium temperature which enters the problem through the specification of a Newtonian form of diabatic heating. The contributions by the large scale eddies are evaluated through the use of observed values for the eddy heat and momentum transports.The contributions from each of the three forcing mechanisms to the temperature and zonal wind fields are invstigated individually and found to be of about equal importance. The sum of the three contributions are also presented for the temperature, the zonal wind, the stream function associated with the mean meridional circulation and the corresponding vertical motion. Although the results fail to reproduce the main observed features of the lower stratosphere, they are found to be in good agreement with observations in the middle latitude troposphere. At any pressure level, for example, the computed mean zonal wind has a jet-like profile and the axis of the jet is found to slope to the south with height, as observed in the atmosphere.Based in part on a thesis submitted by the first author as partial fulfillment of the requirements for the Ph.D. degree at the University of Michigan. — Publication No. 194 from the Department of Meteorology and Oceanography, The University of Michigan.     

On the determination of the non-linearity of the gravimeter scale

Summary The results of two different methods of the determination of the non-linearity of the gravimeter scale are compared; the economy is discussed as well.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.