Bayesian inference of plastosphere viscosities near the Kunlun Fault, northern Tibet

IEECAS SKLLQG

Analytical solution of a multidimensional Langevin equation at high friction limits and probability passing over a two-dimensional saddle

The analytical solution of a multidimensional Langevin equation at the overdamping limit is obtained and the probability of particles passing over a two-dimensional saddle point is discussed. These results may break a path for studying further the fusion in superheavy elements synthesis.

Probability distribution of the four-phase structure invariants of isomorphously related structure factors and its applications

The probability distribution of the four-phase structure invariants (4PSIs) involving four pairs of structure factors is derived by integrating the direct methods with isomorphous replacement (IR). A simple expression of the reliability parameter for 16 types of invariant is given in the case of a native protein and a heavy-atom derivative. Test calculations on a protein and its heavy-atom derivative using experimental diffraction data show that the reliability for 4PSI estimates is comparable with that for the three-phase structure invariants (3PSIs), and that a large-modulus invariants method can be used to improve the accuracy.

Probability distribution of random wave-current forces

Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave-wave and wave-current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan-Pierson-Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. (c) 2006 Elsevier Ltd. All rights reserved.

油气勘探风险评价与决策技术研究

PetroChina and other national petroleum incorporations need rigorous procedures and practical methods in risk evaluation and exploration decision at home and abroad to safeguard their international exploration practice in exploration licence bidding, finding appropriate ratio of risk sharing with partners, as well as avoiding high risk projects and other key exploration activities. However, due to historical reasons, we are only at the beginning of a full study and methodology development in exploration risk evaluation and decision. No rigorous procedure and practical methods are available in our exercises of international exploration. Completely adopting foreign procedure, methods and tools by our national incorporations are not practical because of the differences of the current economic and management systems in China. The objective of this study is to establish a risk evaluation and decision system with independent intellectual property right in oil and gas exploration so that a smooth transition from our current practice into international norm can take place. The system developed in this dissertation includes the following four components: 1. A set of quantitative criteria for risk evaluation is derived on the basis of an anatomy of the parameters from thirty calibration regions national wide as well as the characteristics and the geological factors controlling oil and gas occurrence in the major petroleum-bearing basins in China, which provides the technical support for the risk quantification in oil and gas exploration. 2. Through analysis of existing methodology, procedure and methods of exploration risk evaluation considering spatial information are proposed. The method, utilizing Mahalanobis Distance (MD) and fuzzy logic for data and information integration, provides probabilistic models on the basis of MD and fuzzy logic classification criteria, thus quantifying the exploration risk using Bayesian theory. A projection of the geological risk into spatial domain provides a probability map of oil and gas occurrence in the area under study. The application of this method to the Nanpu Sag shows that this method not only correctly predicted the oil and gas occurrence in the areas where Beibu and Laoyemiao oil fields are found in the northwest of the onshore area, but also predicted Laopu south, Nanpu south and Hatuo potential areas in the offshore part where exploration maturity was very low. The prediction of the potential areas are subsequently confirmed by 17 exploration wells in the offshore area with 81% success, indicating this method is very effective for exploration risk visualization and reduction. 3. On the basis of “Methods and parameters of economic evaluation for petroleum exploration and development projects in China”, a ”pyramid” method for sensitivity analysis was developed, which meets not only the need for exploration target evaluation and exploration decision at home, but also allows a transition from our current practice to international norm in exploration decision. This provides the foundation for the development of a software product “Exploration economic evaluation and decision system of PetroChina” (EDSys). 4. To solve problem in methodology of exploration decision, effort was made on the method of project portfolio management. A drilling decision method was developed employing the concept of geologically risked net present value. This method overcame the dilemma of handling simultaneously both geological risk and portfolio uncertainty, thus casting light into the application of modern portfolio theory to the evaluation of high risk petroleum exploration projects.

On the Dirichlet Prior and Bayesian Regularization

A common objective in learning a model from data is to recover its network structure, while the model parameters are of minor interest. For example, we may wish to recover regulatory networks from high-throughput data sources. In this paper we examine how Bayesian regularization using a Dirichlet prior over the model parameters affects the learned model structure in a domain with discrete variables. Surprisingly, a weak prior in the sense of smaller equivalent sample size leads to a strong regularization of the model structure (sparse graph) given a sufficiently large data set. In particular, the empty graph is obtained in the limit of a vanishing strength of prior belief. This is diametrically opposite to what one may expect in this limit, namely the complete graph from an (unregularized) maximum likelihood estimate. Since the prior affects the parameters as expected, the prior strength balances a "trade-off" between regularizing the parameters or the structure of the model. We demonstrate the benefits of optimizing this trade-off in the sense of predictive accuracy.

Probabilistic Solution of Inverse Problems

In this thesis we study the general problem of reconstructing a function, defined on a finite lattice from a set of incomplete, noisy and/or ambiguous observations. The goal of this work is to demonstrate the generality and practical value of a probabilistic (in particular, Bayesian) approach to this problem, particularly in the context of Computer Vision. In this approach, the prior knowledge about the solution is expressed in the form of a Gibbsian probability distribution on the space of all possible functions, so that the reconstruction task is formulated as an estimation problem. Our main contributions are the following: (1) We introduce the use of specific error criteria for the design of the optimal Bayesian estimators for several classes of problems, and propose a general (Monte Carlo) procedure for approximating them. This new approach leads to a substantial improvement over the existing schemes, both regarding the quality of the results (particularly for low signal to noise ratios) and the computational efficiency. (2) We apply the Bayesian appraoch to the solution of several problems, some of which are formulated and solved in these terms for the first time. Specifically, these applications are: teh reconstruction of piecewise constant surfaces from sparse and noisy observationsl; the reconstruction of depth from stereoscopic pairs of images and the formation of perceptual clusters. (3) For each one of these applications, we develop fast, deterministic algorithms that approximate the optimal estimators, and illustrate their performance on both synthetic and real data. (4) We propose a new method, based on the analysis of the residual process, for estimating the parameters of the probabilistic models directly from the noisy observations. This scheme leads to an algorithm, which has no free parameters, for the restoration of piecewise uniform images. (5) We analyze the implementation of the algorithms that we develop in non-conventional hardware, such as massively parallel digital machines, and analog and hybrid networks.

An in silico approach for screening flavonoids as P-glycoprotein inhibitors based on a Bayesian-regularized neural network

P-glycoprotein (P-gp), an ATP-binding cassette (ABC) transporter, functions as a biological barrier by extruding cytotoxic agents out of cells, resulting in an obstacle in chemotherapeutic treatment of cancer. In order to aid in the development of potential P-gp inhibitors, we constructed a quantitative structure-activity relationship (QSAR) model of flavonoids as P-gp inhibitors based on Bayesian-regularized neural network (BRNN). A dataset of 57 flavonoids collected from a literature binding to the C-terminal nucleotide-binding domain of mouse P-gp was compiled. The predictive ability of the model was assessed using a test set that was independent of the training set, which showed a standard error of prediction of 0.146 +/- 0.006 (data scaled from 0 to 1). Meanwhile, two other mathematical tools, back-propagation neural network (BPNN) and partial least squares (PLS) were also attempted to build QSAR models. The BRNN provided slightly better results for the test set compared to BPNN, but the difference was not significant according to F-statistic at p = 0.05. The PLS failed to build a reliable model in the present study. Our study indicates that the BRNN-based in silico model has good potential in facilitating the prediction of P-gp flavonoid inhibitors and might be applied in further drug design.

Methods to accelerate the learning of bayesian network structures

R. Daly and Q. Shen. Methods to accelerate the learning of bayesian network structures. Proceedings of the Proceedings of the 2007 UK Workshop on Computational Intelligence.

A scenario-driven decision support system for serious crime investigation

Q. Shen, J. Keppens, C. Aitken, B. Schafer, and M. Lee. A scenario driven decision support system for serious crime investigation. Law, Probability and Risk, 5(2):87-117, 2006. Sponsorship: UK Engineering and Physical Sciences Research Council grant GR/S63267; partially supported by grant GR/S98603

Speeding up the learning of equivalence classes of Bayesian network structures

R. Daly, Q. Shen and S. Aitken. Speeding up the learning of equivalence classes of Bayesian network structures. Proceedings of the 10th International Conference on Artificial Intelligence and Soft Computing, pages 34-39.

Using ant colony optimisation in learning Bayesian network equivalence classes

R. Daly, Q. Shen and S. Aitken. Using ant colony optimisation in learning Bayesian network equivalence classes. Proceedings of the 2006 UK Workshop on Computational Intelligence, pages 111-118.