补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.     

位场向下延拓的波数域广义逆算法

位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果.     

位场向下延拓的最小曲率方法

针对位场向下延拓的不适定性,我们将位场向下延拓视为向上延拓的反问题,提出以位场最小曲率作为约束条件来求解稳定的下延位场.我们将剖面位场向上延拓表达式用傅里叶矩阵的形式表示,以矩阵乘法形式给出延拓的表达式,同时向待反演的下延位场引入最小曲率约束,得到向下延拓的最小曲率解,并利用正交变换给出了更为简洁的频率域解.随后,利用Kronecker积将上述全部结果拓展至三维位场,给出了三维位场向下延拓的最小曲率解.此外,我们将位场数据的填充、扩充问题与向下延拓问题统筹考虑,提出一种新的向下延拓迭代格式,该算法面向实际资料处理需求、无须预扩充或填补数据.下延迭代时,对原始数据直接向下延拓,而空白部分利用上一次下延位场估计的上延值替代其空白值并对其向下延拓,直至获得最小曲率约束下稳定的向下延拓结果.同时,我们也讨论了利用改进L曲线和广义交叉验证(GCV)计算正则参数最优估计的问题.对理论模型和实际航空重力资料进行了向下延拓检验,处理结果表明位场向下延拓的最小曲率方法解能满足实际位场资料对向下延拓处理的需求,具有较高的下延精度.     

位场向下延拓三种迭代方法之比较

位场向下延拓在重磁资料解释和用于位场导航的基准数据库构建中发挥着重要作用.本文针对第一类Fredholm积分方程的三种空间域迭代解法:迭代Tikhonov正则化法、Landweber正则化迭代法和积分选代法,基于算子理论和不适定问题的正则化处理方法,首先利用傅里叶变换将空间域迭代法变换到波数域,然后由数学归纳法推导得到这三种迭代法对应的波数域位场向下延拓算子；由Landweber迭代法和积分迭代法在迭代形式上的相似性,探讨了它们在位场向下延拓中的异同及各自优势.模型对比分析表明:(1)两种迭代正则化方法在正则化参数选择合适的条件下,其向下延拓的效果要明显优于积分迭代法,且当收敛到相同误差水平时,迭代Tikhonov正则化法在迭代次数上要远远小于Landweber选代法,但迭代Tikhonov正则化方法存在对正则化参数敏感的问题；(2)从实际应用上讲,由于积分迭代法不存在正则化参数的选择问题,所以该迭代法具有较强的实用性,但需考虑其波数域向下延拓算子时噪声的放大效应.     

位场向下延拓组合滤波器的设计和应用

一 重、磁异常向下延拓能将水平方向迭加的异常分开,并能突出异常细节,从而丰富异常的解释内容,提高异常解释的可靠性.本文从对引起向下延拓不稳定的各种因素的讨论入手,结合对一些向下延拓方法在频率域的分析,提出了一种与已知资料频谱特征相“匹配”的组合滤波器。经理论模型和实例考核,说明该滤波器既能保证向下延拓有较高的精度,又有较强的抗干扰能力。     

位场向下延拓的波数域迭代法及其收敛性

提出了位场向下延拓的波数域迭代法. 对水平面上的位场观测值进行Fourier变换,得到其波谱. 根据第一类Fredholm积分方程的空间域迭代解法,推导出计算向下延拓水平面上位场波谱的波数域迭代公式. 在波数域中进行迭代,一直进行到相继两次迭代近似解的差值最大绝对值小于给定的精度,或迭代达到给定的最大迭代次数. 对这种迭代近似解进行Fourier逆变换,得到向下延拓的位场. 数值计算结果表明：与空间域迭代法比较,这种波数域迭代法简单、快速,并有同样好的向下延拓效果. 本文还证明了这种迭代法是收敛的,并给出了它的收敛特性和滤波特性.     

基于位场分离与延拓的视密度反演

由于目前计算机内存和速度的限制,在对大面积重力资料进行三维密度反演时,已有的反演方法很难奏效.文中提出了一种基于位场分离与延拓的视密度快速反演方法：首先应用场分离的切割法对平面上的重力场进行不同深度层源的切割分离;然后运用大深度向下延拓方法将各层的场延拓至相应的深度;最后反演出各深度层的密度.反演得到的密度是各深度层密度的近似分布,称为视密度反演.本反演方法克服了传统已有方法计算时间特别长、解稳定性差的缺点,在主频1.99 GHz的微机上,反演128×128×10个密度值的计算时间小于20 s,理论模型和实际资料的应用表明反演具有较好的效果.     

海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.     

位场向下延拓的改进迭代维纳滤波法

根据维纳滤波理论导出的位场向下延拓滤波器为最佳下延滤波器,但因其实现需要已知待求位场和噪声的功率谱而在实际应用中受到限制.针对该问题,本文首先提出一种基于位场径向平均功率谱的位场噪声水平估计方法,进而利用偏差准则求取正则化参数,实现位场正则化向下延拓;然后将位场正则化下延结果的功率谱作为待求位场功率谱的估计初值,采用带修正项的迭代维纳滤波方法来更新对待求位场功率谱的估计,最后提出本文的位场向下延拓改进迭代维纳滤波方法.基于理论重力模型数据及航磁实测数据进行了向下延拓对比试验,结果表明,改进迭代法具有较好的收敛性,且下延精度优于Tikhonov正则化法和递增型维纳滤波法.     

Changes in sediment transport in the Kuye River in the Loess Plateau in China

In this paper, the changes in sediment transport over 51 years from 1955 to 2006 in the Kuye River in the Loess Plateau in China are assessed. Key factors affecting sediment yield and sediment transport, such as precipitation depth, discharge, and human activities are studied. To investigate the changes in sediment yield in this watershed, a trend analysis on sediment concentration, precipitation depth, and discharge is conducted. Precipitation depths at 2 Climate Stations （CSs）, as well as discharge and sediment transport at 3 Gauging Stations （GSs） are used to assess the features of sediment transport in the Kuye River. The rtmoff modulus （defined as the annual average discharge per unit area, L/（s·km^2）） and the sediment transport modulus （defined as the annual suspended sediment transport per unit area, t/（yr km^2）） are introduced in this study to assess the changes in runoff and sediment yield for this watershed. The results show that the highest average monthly discharge during the study period in the Kuye River is 66.23 m^3/s in August with an average monthly sediment concentration of 88.9 kg/m^3. However, the highest average monthly sediment concentration during the study period in the Kuye River is 125.34 kg/m^3 and occurs in July, which has an average discharge of 42.6 m^3/s that is much less than the average monthly discharge in August. It is found that both the runoff modulus and sediment transport modulus at Wenjiachuan GS on the Kuye River has a clear downward trend. During the summer season from July to August, the sediment transport modulus at Wenjiachuan GS is much higher than those at Toudaoguai and Longmen GSs on the Yellow River. The easily erodible loess in the Kuye River watershed and the sparse vegetation are responsible for the extremely high sediment yield from the Kuye River watershed. The analyses of the grain size distribution of suspended load in the Kuye River are presented. The average monthly median grain size of suspended load in the Kuye River is largest in February and then decreases until June. In July, the average monthly median grain size of suspended load approaches another peak and decreases until September. Then, the median grain size of suspended load starts to increase until February of the following year. However, the average monthly median grain size of suspended load in the Yellow River at Toudaoguai and Longmen GSs is the smallest between early summer and late fall The median grain size in the Yellow River starts to increase in November and approaches the largest size in January.     

Variations in trace gas concentrations in different environments in India

Observations of trace gases (SO₂, NH₃, NO₂ and O₃) were made during the period 1981 to 1984 at 6 different locations representative of urban industrial, urban, nonurban, thermal power plant and marine environment. Diurnal variations of the trace gases were studied in an urban environment. Except in the urban industrial environment, the concentration of NH₃ was found in the range of background values. Also, the average concentrations of NO₂ and O₃ at the different environments were in the order of background values. However, the concentrations of SO₂ were substantially higher by about 7 times, in urban industrial and thermal power plant environments. The diurnal variations of SO₂, NH₃ and NO₂ showed anitphase relationship with surface temperature at the urban environment station which is relatively free of industrial pollution. Discussion is centred on trace gas variations in different environments in India together with the values reported for various countries in the world.     

TritiogenicHe in groundwater in Takaoka

³He/⁴He ratios up to 3.5 times the ratio of atmospheric He were found in groundwater samples. The³He enrichment can be attributed to radiogenic³He produced by in-situ beta-decay of³H. This shows that tritiogenic³He is accumulating in confined waters. From tritiogenic³He and³H concentrations, ages of groundwaters can be calculated. Detection of tritiogenic³He gives a tool to trace a tritium contamination which occurred in the past and cannot be assessed only by the³H counting method.     

2008年云南地震灾害概况

列出了2008年云南省破坏性地震目录.在此基础上,结合云南省地震局的地震灾害评估资料,总结了2008年云南地震灾害的主要数据和特性.最后列出1990～2008年云南地震灾害的相应数据,并进行了简要对比.     

Applied volcanology in geothermal exploration in Iceland

In an active spreading area like Iceland, where the regional geothermal gradient is in the range 50–150°C/km, it is normally not a problem to find high enough temperature with deep drilling, but the difficulties arise with finding permeable layers at depth within the strata. Various volcanological methods can be applied in the search for aquifers and geothermal reservoir rocks. The flow pattern (as deduced from deuterium studies) indicates that the thermal water flows preferentially along high porosity stratiform horizons and dyke swarms from the recharge areas in the highlands to the hot spring areas in the lowlands. The primary porosity of the volcanic strata is dependent on the chemical composition and the mode of eruption of the volcanic units. Both the reservoir rocks and the flow channels forming the geothermal plumbing system are thought to vary from the Tertiary to the Quaternary provinces due to environmental conditions at the eruptive sites.     

2007年中国大陆地震灾害损失述评

在列出2007年中国5级以上地震目录的基础上,结合有关省(自治区、直辖市)地震局的地震灾害评估资料,总结出2007年中国大陆地震灾害的主要数据和特性。最后列出1990—2007年中国大陆地震灾害的相应数据,并进行了简要的比较。     

Projected changes in flood indices in selected catchments in Poland in the 21st century

The aim of this study is to estimate likely changes in flood indices under a future climate and to assess the uncertainty in these estimates for selected catchments in Poland. Precipitation and temperature time series from climate simulations from the EURO-CORDEX initiative for the periods 1971–2000, 2021–2050 and 2071–2100 following the RCP4.5 and RCP8.5 emission scenarios have been used to produce hydrological simulations based on the HBV hydrological model. As the climate model outputs for Poland are highly biased, post processing in the form of bias correction was first performed so that the climate time series could be applied in hydrological simulations at a catchment-scale. The results indicate that bias correction significantly improves flow simulations and estimated flood indices based on comparisons with simulations from observed climate data for the control period. The estimated changes in the mean annual flood and in flood quantiles under a future climate indicate a large spread in the estimates both within and between the catchments. An ANOVA analysis was used to assess the relative contributions of the 2 emission scenarios, the 7 climate models and the 4 bias correction methods to the total spread in the projected changes in extreme river flow indices for each catchment. The analysis indicates that the differences between climate models generally make the largest contribution to the spread in the ensemble of the three factors considered. The results for bias corrected data show small differences between the four bias correction methods considered, and, in contrast with the results for uncorrected simulations, project increases in flood indices for most catchments under a future climate.     

冬季亚洲中强震与夏季我国东部主雨带的准同纬性

根据1954～2003年50年亚洲中强震与我国地温资料,分析了冬季震中纬度—震涡中心纬度—冬季我国东部强地热涡中心纬度—夏季我国东部主雨带纬度4者之间的相互关系,发现冬季中强震与夏季主雨带之间具有较好的准同纬性. 这是每年汛期降水预报的一个指标. 文中也讨论了应用此预报指标时的一些特例.     

Advances in ground motion studies in China

This paper briefly summarizes the works in the processing of strong ground motion data, the factors affecting strong ground motion, the modeling of strong ground motion and the calculating of broad-band response spectrum which have been done recent years by engineering seismologists and seismologists of China. In addition, we think back to the international cooperation in strong ground motion of the recent years and make some expectations for the future.