Fully implicit black-oil simulations result in huge, often very-illconditioned, linear systems of equations for different unknowns (e.g., pressure and saturations). It is well-known that the underlying Jacobian matrices contain both hyperbolic and nearly elliptic subsystems (corresponding to saturations and pressure, respectively). Because a reservoir simulation is typically driven by the behavior of the pressure, constrained-pressure-residual (CPR)type two-stage preconditioning methods to solve the coupled linear systems are a natural choice and still belong to the most popular approaches. After a suitable extraction and decoupling, the computationally most costly step in such two-stage methods consists in solving the elliptic subsystems accurately enough. Algebraic multigrid (AMG) provides a technique to solve elliptic linear equations very efficiently. Hence, in recent years, corresponding CPR-AMG approaches have been extensively used in practice.Unfortunately, if applied in a straightforward manner, CPR-AMG does not always work as expected. In this paper, we discuss the reasons for the lack of robustness observed in practice, and present remedies. More precisely, we will propose a preconditioning strategy (based on a suitable combination of left and right preconditioning of the Jacobian matrix) that aims at a compromise between the solvability of the pressure subproblem by AMG and the needs of the outer CPR process. The robustness of this new preconditioning strategy will be demonstrated for several industrial test cases, some of which are very ill-conditioned. Furthermore, we will demonstrate that CPR-AMG can be interpreted in a natural way as a special AMG process applied directly to the coupled Jacobian systems.

Short-term metro passenger flow forecasting is an essential component of intelligent transportation systems (ITS) and can be applied to optimize the passenger flow organization of a station and offer data support for metro passenger flow early warning and system management. LSTM neural networks have recently achieved remarkable recent in the field of natural language processing (NLP) because they are well suited for learning from experience to predict time series. For this purpose, we propose an empirical mode decomposition (EMD)-based long short-term memory (LSTM) neural network model for predicting short-term metro inbound passenger flow. The EMD algorithm decomposes the original sequential passenger flow into several intrinsic mode functions (IMFs) and a residual. Selected IMFs that are strongly correlated with the original data can be obtained via feature selection. The selected IMFs and the original data are integrated into inputs for LSTM neural networks, and a single LSTM prediction model and an EMD-LSTM hybrid forecasting model are developed. Finally, historical real automatic fare collection (AFC) data from metro passengers are collected from Chengdu Metro to verify the validity of the proposed EMD-LSTM prediction model. The results indicate that the proposed EMD-LSTM hybrid forecasting model outperforms the LSTM, ARIMA and BPN models.

This paper studies on the dual-channel (the traditional channel and the E-commerce channel) supply chain network design (SCND) for fresh agri-product (FAP) under information uncertainty. The model is to solve the integration network design issues about production, supply, and sales of FAP, as well as to minimize the supply chain operation cost and maximizing the satisfaction degree of the logistics demand between supply chain (SC) nodes simultaneously. First, the triangular-fuzzy-number is used to depict the information uncertainty in the SCND. Then the handling cost of wasted FAP and quantitative traceability cost of FAP are considered in addition to the regular transportation cost and fixed cost for facility location, which makes it more closely and suitable for the real case when compared with the models in existing literatures. Due to there is always a lower and upper bound for delivery time of each logistics demand node in SC in real life. This problem is then formulated as an attenuation function, for depicting the satisfaction degree in each demand node. Afterwards, the multi-objective planning method is utilized to represent the mutual trade-off among supply-chain participants. The proposed models are solved with two-phase method to guarantee the optimality of the solution. Finally, the effectiveness and applicability of the proposed models and algorithm are validated with a real case. The result shows

No abstract

Fully implicit petroleum reservoir simulations result in huge, often very ill-conditioned linear systems of equations to solve for different unknowns, for example, pressure and saturations. It is well known that the full system matrix contains both hyperbolic as well as nearly elliptic sub-systems. Since the solution of the coupled system is mainly determined by the solution of their elliptic (typically pressure) components, (CPR-type) two-stage preconditioning methods still belong to the most popular approaches to tackle such coupled systems. After a suitable extraction and decoupling, the numerically most costly step in such two-stage methods consists in solving these elliptic sub-systems. It is known that algebraic multigrid (AMG) provides a technique to solve elliptic linear equations very efficiently. The main advantage of AMG-based solvers -their numerical scalability -makes them particularly efficient for solving huge linear systems.Depending on the application, the system's properties range from simple to highly indefinite. Unfortunately decoupling pressure and saturation related parts may introduce further difficulties. Consequently, in complex industrial simulations, the application of AMG to elliptic sub-systems might not be straightforward. In fact, an important goal in defining an efficient two-stage preconditioning strategy consists in extracting elliptic sub-systems that are suitable for an efficient AMG solution and, at the same time, ensure a fast overall convergence of the two-stage approach.The importance of this will be demonstrated for several industrial cases. In particular, some of these cases are very hard to solve by AMG if applied in a standard way.Preliminary results for a CPR-type coupling of SAMG to CMG's PARASOL, a variable degree variable ordering ILU preconditioner using FGMRES, are compared to using PARASOL by itself. Alternative preconditioning operators will be presented giving elliptic sub-systems which are not only more suitable for applying AMG efficiently but also help accelerate the CPR-type process. Comparisons with one-level iterative methods will show the acceleration by AMG is highly superior. Finally, a strategy is presented that combines all linear solver parts in one single AMG-iteration. In this sense CPR can be seen as a special case of AMG for systems. This, in turn, yields a -formally -very simple but simultaneously very flexible solution approach.