Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Environmental impacts of bench blasting at Hisarcik Boron open pit mine in Turkey

This paper presents the results of ground vibration measurements carried out in Hisarcik Boron open pit mine located on the west side of central Anatolia near Kütahya province in Turkey. Within the scope of this study to predict peak particle velocity (PPV) level for this site, ground vibration components were measured for 304 shots during bench blasting. In blasting operations, ANFO (blasting agent), gelatin dynamite (priming), and delay electric detonators (firing) were used as explosives. Parameters of scaled distance (charge quantity per delay and the distance between the source and the station) were recorded carefully and the ground vibration components were measured for all blast events using two different types of vibration monitors (one White Mini-Seis and one Instantel Minimate Plus Model). The absolute distances between shot points and monitor stations were determined using GPS. The equation of square root scaled distance extensively used in the literature was taken into consideration for the prediction of PPV. Then, the data pairs of scaled distance and particle velocity obtained from the 565 event records were analyzed statistically. At the end of statistical evaluation of the data pairs, an empirical relation which gives 50% prediction line with a reasonable correlation coefficient was established between PPV and scaled distance.     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

Early Estimation of Seismic Hazard for Strong Earthquakes in Taiwan

A shakemap system providing rapid estimates of strong ground shaking could be useful for emergency response providers in a damaging earthquake. A hybrid procedure, which combines site-dependent ground motion prediction models and the limited observations of the Real-Time Digital stream output system (RTD system operated by Central Weather Bureau, CWB), was set up to provide a high-resolution shakemap in a near-real-time manner after damaging earthquakes in Taiwan. One of the main factors that affect the result of ground motion prediction analysis is the existence of site effects. The purpose of this paper is to investigate the local site effects and their influence in the ground shaking and then establish an early estimation procedure of potential hazard for damaging earthquakes. Based on the attenuation law, the site effects of each TSMIP station are discussed in terms of a bias function that is site and intensity-level dependent function. The standard deviation of the site-dependent ground motion prediction model can be significantly reduced. The nonlinear behavior of ground soil is automatically taken into account in the intensity-level dependent bias function. Both the PGA and the spectral acceleration are studied in this study. Based on the RTD data, event correctors are calculated and applied to precisely estimate the shakemap of damaging earthquakes for emergency response.     

二维电子地图和三维虚拟场景的互响应

本文提出了基于二维电子地图与三维虚拟场景之间相互响应的思想来克服它们各自的局限性.简要的描述了相互响应的原理、表现形式、实现流程及实现方法.最后提出了一些结论和展望.     

一种层状大地瞬变电磁响应正演计算的改进方法

频谱法是瞬变电磁法正演计算中常用的计算方法之一,其关键在于内层含贝塞尔函数的积分计算。对于重叠回线装置的瞬变电磁法正演的内层积分,以往采用在积分区间内寻找贝塞尔函数的零点分布,依次在相邻零点之间采用一般的数值积分法在求得各自的积分值后再叠加的方法进行计算。这种方法精度较高,但效率低。利用贝塞尔函数的大宗量渐近特性,提出一种新的计算方法。计算结果表明,该方法计算效率高、方法简单、精度可靠。     

超声与振荡提取沉积物中松散结合态磷的比较

用超声技术提取沉积物样品中松散结合态磷,并与振荡提取结果进行比较。结果表明,超声提取结果与振荡提取结果大致相当。超声提取耗时短,操作简单易行,测试结果精密度为3.8%(RSD,n=7),可以替代振荡方法提取沉积物中的松散结合态磷。     

WAFS产品中GRIB资料中国区产品评估

为了给使用WAFS产品的用户提供量化参考依据，选取WAFS产品中常用区域和层次的GRIB数据，利用由国家气象中心提供的风、温客观分析网格点资料，与WAFS中同时刻的预报场产品（风、温网格点资料），用均方根误差进行数字化形式分析比较。结果表明：WAFS提供的风、温预报，通常是短时效的风、温预报比长时效的风、温预报更接近客观分析场；低层的预报比高层的预报更接近客观分析场；风的预报以v矢量的预报优于u矢量的预报；风的误差主要来源于u矢量的误差。     

冬季浙闽沿岸水分布的短期变动与风的关系初探

利用1999年7月至2003年5月期间的遥感数据,包括AVHRR海表层温度、QuikSCAT风场和风应力数据,在分析4年内月平均遥感温度场和风场特征与历年现场观测所获得的认识一致的基础上,选取2002年002-008天这一连续晴空的时段,尝试建立简单的沿岸冷水影响面积表征方法,初步探讨了冬季台湾海峡浙闽沿岸水分布的短期变动与风应力之间的关系。结果表明,风是决定冬季台湾海峡海表层温度逐日变动的关键因素,日平均SST与风应力的相关系数R2达到0·90。采用温度法(SST≤17℃)和温度空间距平法(≤-1℃)表征的浙闽沿岸水影响面积的变化趋势基本一致,而且影响面积的逐日变动与风应力显著相关,二者的相关系数R2分别达到0·90和0·91。     

Wind and Wave Data Analysis for the Aegean Sea - Preliminary Results

In this paper, a statistical analysis on wind and wave buoy measurements and wind and wave model forecasts obtained during a two-year period (1999-2001) is presented with reference to four characteristic near-shore sites of the Aegean Sea. The measurements are a main product of the "POSEIDON" system aiming at the monitoring and forecasting of the state of the Greek seas, operated by the National Centre for Marine Research (NCMR). Although the two-year period is rather short for a thorough analysis of the local wind and wave climate, yet the obtained results, presented herein for the first time, reveal some interesting features of the corresponding wave and wind characteristics. Comparisons between the measurements and the forecast results are also performed at the locations under consideration. It is found that (i) wind speeds obtained from the POSEIDON weather forecasting system are, in general, in agreement with the measurements, except for high wind speeds which are systematically underestimated, (ii) the WAM model can successfully follow the monthly and over year trend of the evolution of wind and wave characteristics, but face significant problems for efficient sea-state forecasting. Finally, the overall pattern of the wind/wave climate for the entire Aegean Sea as obtained from the models is presented by means of the spatial distribution of the mean annual wind and sea-state intensity.