Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

角度格式转换过程中可靠算法研究

提出了一种可靠算法是将原始度.分秒加上一个与之同符号,且不影响计算精度的有效微小量ｅ,然后再进行格式转换,能有效减少工程返工和损失。     

C1相干算法及其在三维地震资料处理解释中的应用

C1相干算法是计算每道的横测线和相邻纵测线的互相关，依据相干值来判断地层和岩性的横向不均匀性。根据互相关原理，提出了相干体算法中二维和三维C1算法的内容及其具体实现。通过对实际三维地震资料的处理，结合已知地质资料，证明C1相干技术在三维地震资料的断层解释方面有着明显的效果。     

A recursive decomposition algorithm for network seismic reliability evaluation

A new probabilistic analytical approach to evaluate seismic system reliability of large lifeline systems is presented in this paper. The algorithm takes the shortest path from the source to the terminal of a node weight or edge weight network as decomposition policy, using the Boolean laws of set operation and probabilistic operation principal, a recursive decomposition process then could be constructed. For a general weight network, the modified Torrieri method (NTR/T method) is introduced to combine with the suggested algorithm. Therefore, the recursive decomposition algorithm may be applied to evaluate the seismic reliability of general lifeline systems. A series of case studies, including a practical district electric power network system and a large urban water supply system, show that the suggested algorithm supplies a useful probabilistic analysis means for the seismic reliability evaluation of large lifeline systems. Copyright © 2002 John Wiley & Sons, Ltd.     

应用CP网络进行岩性识别

为通过测井解决岩性识别问题，引入了具有分类准确、算法简练等优点的CP（Counter-Propagation)网络。在详细介绍CP网络的网络模型和算法的基础上，结合某油田的实际测井资料，进行了CP网络识别研究。应用结果表明：CP网络训练周期短、识别准确率高、不存在收敛问题。通过试验研究得出结论：CP网络完全可以用于解决岩性识别等问题，具有广阔的应用前景。     

遗传算法在边坡数值计算中的应用

改进了进化方向的遗传算法与有限元数值法，结合并研制了相应的软件。应用该软件对多类型岩土边坡进行弹性模量、内聚力、内摩擦角等参数反演分析，显示误差很小，收敛速度也很快，这说明改进进化方向遗传算法这种新型的优化算法在多类型岩土参数优化估计中具有独特的优势。     

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.     

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.