Modeling selecting models

Stage 5 (Building the Model). This involves the classical steps of defining the conceptual model, selecting model code, calibrating the model against field data and... 

The guidance material available via WISF-RTD covers more general topics such as guidance on the use of tools (e.g. model selection, model linking, model calibration and validation, model sensitivity and uncertainty analysis, etc.), the monitoring process, the stakeholder participation, etc. Such technical guidances are set up by various RTD... 

The results of this case study clearly showed that none of the packages investigated contained all the advanced statistical tools to perform the requested tasks concerning model selection, model discrimination and experimental design. 

The next step is to perform model selection. Models may be rejected for three different reasons (i) because the differences between the experimental data and the data calculated with the fitted model are much larger than the measurement error (the model is then qualified as inadequate ), (ii) because the fit of the model is significantly worse than an alternative model, and (iii) because one or more parameters in the kinetic model cannot be estimated accurately and independently, which usually indicates that the model contains too many parameters. Although there were large differences between the options available in the packages investigated, none of the packages were capable to perform all these checks. 

The most reliable estimates of the parameters are obtained from multiple measurements, usually a series of vapor-liquid equilibrium data (T, P, x and y). Because the number of data points exceeds the number of parameters to be estimated, the equilibrium equations are not exactly satisfied for all experimental measurements. Exact agreement between the model and experiment is not achieved due to random and systematic errors in the data and due to inadequacies of the model. The optimum parameters should, therefore, be found by satisfaction of some selected statistical criterion, as discussed in Chapter 6. However, regardless of statistical sophistication, there is no substitute for reliable experimental data. 

Two generally accepted models for the vapor phase were discussed in Chapter 3 and one particular model for the liquid phase (UNIQUAC) was discussed in Chapter 4. Unfortunately, these, and all other presently available models, are only approximate when used to calculate equilibrium properties of dense fluid mixtures. Therefore, any such model must contain a number of adjustable parameters, which can only be obtained from experimental measurements. The predictions of the model may be sensitive to the values selected for model parameters, and the data available may contain significant measurement errors. Thus, it is of major importance that serious consideration be given to the proper treatment of experimental measurements for mixtures to obtain the most appropriate values for parameters in models such as UNIQUAC. 

Table 5.4 gives the specific energies of selected organic liquid compounds. Compared with the isooctane chosen as the base reference, the variations from one compound to another are relatively small, on the order of 1 to 5%, with the exception of some particular chemical structures such as those of the short chain nitroparaffins (nitromethane, nitroethane, nitropropane) that are found to be energetic . That is why nitromethane, for example, is recommended for very small motors such as model airplanes it was also used in the past for competitive auto racing, for example in the Formula 1 at Le Mans before being forbidden for safety reasons. 

The reservoir model will usually be a computer based simulation model, such as the 3D model described in Section 8. As production continues, the monitoring programme generates a data base containing information on the performance of the field. The reservoir model is used to check whether the initial assumptions and description of the reservoir were correct. Where inconsistencies between the predicted and observed behaviour occur, the model is reviewed and adjusted until a new match (a so-called history match ) is achieved. The updated model is then used to predict future performance of the field, and as such is a very useful tool for generating production forecasts. In addition, the model is used to predict the outcome of alternative future development plans. The criterion used for selection is typically profitability (or any other stated objective of the operating company). 

Reservoir pressure is measured in selected wells using either permanent or nonpermanent bottom hole pressure gauges or wireline tools in new wells (RFT, MDT, see Section 5.3.5) to determine the profile of the pressure depletion in the reservoir. The pressures indicate the continuity of the reservoir, and the connectivity of sand layers and are used in material balance calculations and in the reservoir simulation model to confirm the volume of the fluids in the reservoir and the natural influx of water from the aquifer. The following example shows an RFT pressure plot from a development well in a field which has been producing for some time. 

To prepare simulated free-of-noise 3D images of a complex body For this target the image of an internal pore in the real welding joint with extracted noise was used. The ray tracer model was applied for the simulation of five projections of the selected image. 

Micellar structure has been a subject of much discussion [104]. Early proposals for spherical [159] and lamellar [160] micelles may both have merit. A schematic of a spherical micelle and a unilamellar vesicle is shown in Fig. Xni-11. In addition to the most common spherical micelles, scattering and microscopy experiments have shown the existence of rodlike [161, 162], disklike [163], threadlike [132] and even quadmple-helix [164] structures. Lattice models (see Fig. XIII-12) by Leermakers and Scheutjens have confirmed and characterized the properties of spherical and membrane like micelles [165]. Similar analyses exist for micelles formed by diblock copolymers in a selective solvent [166]. Other shapes proposed include ellipsoidal [167] and a sphere-to-cylinder transition [168]. Fluorescence depolarization and NMR studies both point to a rather fluid micellar core consistent with the disorder implied by Fig. Xm-12. 

Figure Al.6.26. Stereoscopic view of ground- and excited-state potential energy surfaces for a model collinear ABC system with the masses of HHD. The ground-state surface has a minimum, corresponding to the stable ABC molecule. This minimum is separated by saddle points from two distmct exit chaimels, one leading to AB + C the other to A + BC. The object is to use optical excitation and stimulated emission between the two surfaces to steer the wavepacket selectively out of one of the exit chaimels (reprinted from [54]).

These equations lead to fomis for the thermal rate constants that are perfectly similar to transition state theory, although the computations of the partition functions are different in detail. As described in figrne A3.4.7 various levels of the theory can be derived by successive approximations in this general state-selected fomr of the transition state theory in the framework of the statistical adiabatic chaimel model. We refer to the literature cited in the diagram for details. 

Wyatt R E, Hose G and Taylor H S 1983 Mode-selective multiphoton excitation in a model system Phys. Rev. A 28 815-28... 

Table Bl.13.1 Selected dynamic models used to calculate spectral densities.

It is difficult to observe tliese surface processes directly in CVD and MOCVD apparatus because tliey operate at pressures incompatible witli most teclmiques for surface analysis. Consequently, most fundamental studies have selected one or more of tliese steps for examination by molecular beam scattering, or in simplified model reactors from which samples can be transferred into UHV surface spectrometers witliout air exposure. Reference [4] describes many such studies. Additional tliemes and examples, illustrating botli progress achieved and remaining questions, are presented in section C2.18.4. 

At this stage, we would like to mention that the model, without the vector potential, is constructed in such a way that it obeys certain selection rules, namely, only the even —> even and the odd —> odd transitions are allowed. Thus any deviation in the results from these selection rules will be interpreted as a symmetry change due to non-adiabatic effects from upper electronic states. 

By using this approach, it is possible to calculate vibrational state-selected cross-sections from minimal END trajectories obtained with a classical description of the nuclei. We have studied vibrationally excited H2(v) molecules produced in collisions with 30-eV protons [42,43]. The relevant experiments were performed by Toennies et al. [46] with comparisons to theoretical studies using the trajectory surface hopping model [11,47] fTSHM). This system has also stimulated a quantum mechanical study [48] using diatomics-in-molecule (DIM) surfaces [49] and invoicing the infinite-onler sudden approximation (lOSA). 

Fig. 7. Snapshots of rupture taken (A) at the start of the simulation (zcant = 0), (li) at ZcB.nl = 2.8 A, (C) at Zcnm = 4.1 A, (D) at Zcnm = 7.1 A, and (E) at Zcant = 10.5 A. The biotin molecule is drawn as a ball-and-stick model within the binding )ocket (lines). The bold dashed lines show hydrogen bonds, the dotted lines show selected water bridges.

The LIN method (described below) was constructed on the premise of filtering out the high-frequency motion by NM analysis and using a large-timestep implicit method to resolve the remaining motion components. This technique turned out to work when properly implemented for up to moderate timesteps (e.g., 15 Is) [73] (each timestep interval is associated with a new linearization model). However, the CPU gain for biomolecules is modest even when substantial work is expanded on sparse matrix techniques, adaptive timestep selection, and fast minimization [73]. Still, LIN can be considered a true long-timestep method. 

The two sources of stochasticity are conceptually and computationally quite distinct. In (A) we do not know the exact equations of motion and we solve instead phenomenological equations. There is no systematic way in which we can approach the exact equations of motion. For example, rarely in the Langevin approach the friction and the random force are extracted from a microscopic model. This makes it necessary to use a rather arbitrary selection of parameters, such as the amplitude of the random force or the friction coefficient. On the other hand, the equations in (B) are based on atomic information and it is the solution that is approximate. For ejcample, to compute a trajectory we make the ad-hoc assumption of a Gaussian distribution of numerical errors. In the present article we also argue that because of practical reasons it is not possible to ignore the numerical errors, even in approach (A). 

Nevertheless, chemists have been planning their reactions for more than a century now, and each day they run hundreds of thousands of reactions with high degrees of selectivity and yield. The secret to success lies in the fact that chemists can build on a vast body of experience accumulated over more than a hundred years of performing millions of chemical reactions under carefully controlled conditions. Series of experiments were analyzed for the essential features determining the course of a reaction, and models were built to order the observations into a conceptual framework that could be used to make predictions by analogy. Furthermore, careful experiments were planned to analyze the individual steps of a reaction so as to elucidate its mechanism. 

---