三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

Seismic structural control using semi-active tuned mass dampers

This paper focuses on how to determine the instantaneous damping of the semi-active tuned mass damper (SATMD) with continuously variable damping. An off-and-towards-equilibrium (OTE) algorithm is employed to examine the control performance of the structure/SATMD system by considering the damping as an assumptive control action. The damping modification of the SATMD is carried out according to the proposed OTE algorithm, which is formulated based on analysis of the structural movement under external excitations, and the measured responses of the structure at every time instant. As examples two numerical simulations of a five-storey and a ten-storey shear structures with a SATMD on the roof are conducted. The effectiveness on vibration reduction of MDOF systems subjected to seismic excitations is discussed. Analysis results show that the behavior of the structure with a SATMD is significantly improved and the feasibility of applying the OTE algorithm to the structural control design of SATMD is also verified.     

一种求气井二项式方程系数的新方法

针对开发早期气田缺乏产能试井资料的实际情况，由于无法用传统的压力降法计算产能方程，故而基于最优化算法，提出了一种利用稳定生产数据和井控储量求取气井产能方程系数的新方法。在只知道稳定生产数据情况下便可进行计算分析。该方法在气田进行了实际应用，并和实际资料以及传统试凑法计算的结果进行了比较，表明方法的计算结果较为可靠。     

基坑支护结构水平变形预测的遗传神经网络方法

采用遗传算法和误差反向传播算法相结合的混合算法来训练前馈人工神经网络,先用遗传学习算法进行全局训练,再用BP算法进行精确训练。就遗传算法过程中的选择、变异进行了探索,提出了用BP网络训练产生变异的遗传算法。作为实例,将该方法应用于预测基坑支护结构水平变形中。结果表明,该方法有收敛速度较快、预测精度高等优点。     

等离子体发射光谱法直接测定海洋沉积物中的微量稀土元素

应用全谱直读电感耦合等离子体发射光谱法对海洋沉积物中微量稀土元素的测定进行了研究，选择了最佳的仪器工作条件，采用多谱线拟合技术进行校正，样品的加标回收率为94．0％-110．0％。方法采用海底沉积物国家一级标准物质进行验证，测定结果与标准值较吻合，各元素10次测定的RSD≤11．3％。     

单道扫描电感耦合等离子体发射光谱法测定珊瑚礁中主量和微量元素

建立了用单道扫描电感耦合等离子体发射光谱法测定珊瑚礁中主量和微量元素的方法。珊瑚礁样品经酸溶解后,直接进行Ca、Mg、K、Na、V、Sr、Ba、Co、Ni、Pb、Li、Rb的测定,方法检出限为0.0004～0.1mg/L,精密度RSD<5%(n=6),样品的加标回收率为97%～105%。测定w(CaO)为50%的样品,结果与传统的EDTA容量法相符。方法经碳酸盐岩石国家一级标准物质验证可行,已用于大批量的珊瑚礁样品的测定。     

浅水方程数值计算方法的研究

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.     

加卸载响应比(LURR)理论的实验研究

利用大尺度岩石试件破坏声发射实验，对加卸载响应比(Load/Unload Response Ratio)理论进行了实验研究。实验采用的加载条件力求能够模拟地下岩石的复杂受力过程。在加载过程中，岩石试件内部出现的损伤用声发射技术进行了记录。实验数据经过分析之后所得的结论与加卸载响应比理论预测的结果相符合。把加卸载响应比值Y的变化情况与声发射事件定位的过程相对比，也可以进一步看出岩石材料内部的损伤发展过程与加卸载响应比Y值的变化具有一致性。而高Y值的异常情况出现，又是在材料发生最终破坏之前，所以，可以用加卸载响应比理论定量地分析岩石试件损伤演化情况，从而对岩石材料的破坏及地震做出预测。实验结果再一次验证了加卸载响应比理论。     

The inversion of anelastic coefficient,source parameters and site respond using genetic algorithm

It gradually becomes a common work using large seismic wave data to obtain source parameters, such as seismic moment, break radius, stress drop, with completingof digital seismic network in China (Hough, et al, 1999; Bindi, et al, 2001). These parameters are useful on earthquake prediction and seismic hazard analysis.Although the computation methods of source parameters are simple in principle and the many research works have been done, it is not easy to obtain the parameters accurately. There are two factors affecting the stability of computation results. The first one is the effect of spread path and site respond on signal. According to the research results, there are different geometrical spreading coefficients on different epicenter distance. The better method is to introduce trilinear geometrical spreading model (Atkinson, Mereu, 1992; Atkinson, Boore, 1995; WONG, et al, 2002). In addition, traditional site respond is estimated by comparing with rock station, such as linear inversion method (Andrews, 1982), but the comparative estimation will introduce some errors when selecting different stations. Some recent research results show that site respond is not flat for rock station (Moya, et al, 2000; ZHANG,. et al, 2001; JIN, et al, 2000; Dutta, et al, 2001). The second factor is to obtain low-frequency level and corner frequency fromdisplacement spectrum. Because the source spectrum model is nonlinear function,these values are obtained by eye. The subjectivity is strong. The small change of corner frequency will affect significantly the result of stress drop.     

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.