三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

黄东海及朝鲜海峡冰后期陆架沉积物的环境磁学研究

There are four chapters in this dissertation. The first chapter briefly synthesizes the basic theories, methods and present-day applying situation of environmental magnetism. The second chapter probes into the magnetic mineral diagenesis in the post-glacial muddy sediments from the southeastern South Yellow Sea and its response to marine environmental changes, using the muddy sediment of Core YSDP103 formed in the shelf since about 13 ka BP. The third chapter illustrates the high-resolution early diagenetic processes by investigating the rapidly deposited muddy sediments during the last 6 ka in Cores SSDP-102 and SSDP-103 from the near-shore shelf of Korea Strait. The fourth chapter presents the results of detailed rock magnetic investigation of the surface sediments from the fine-grained depositional area on the outer shelf of the East China Sea in an attempt to provide environmental magnetic evidence for the provenance of the fine-grained deposit. Core YSDP103 was retrieved in the muddy deposit under the cold eddy of the southeastern South Yellow Sea, and the uppermost 29.79 m core represents the muddy sediments formed in the shelf since about 13 ka BP. The lower part from 29.79 to 13.35 m, called Unit A2, was deposited during the period from the post-glacial transgression to the middle Holocene (at about 6 ~(14)C ka BP) when the rising sea level reached its maximum, while the upper part above 13.35 m (called Unit Al) was deposited in a cold eddy associated with the formation of the Yellow Sea Warm Current just after the peak of post-glacial sea level rise. For the the uppermost 29.79 m core, detailed investigation of rock-magnetic properties and analyses of grain sizes and geochemistry were made. The experimental results indicate that the magnetic mineralogy of the core is dominated by magnetite, maghemite and hematite and that, except for the uppermost 2.35 m, the magnetic minerals were subject to reductive diagenesis leading to significant decline of magnetic mineral content and the proportion of low-coercivity component. More importantly, ferrimagnetic iron sulphide (greigite) is found in Unit A2 but absent in Unit Al, suggesting the control of marine environmental conditions on the magnetic mineral diagenesis. Magnetic parameters show abrupt changes across the boundary between the Unit Al and A2, which reflects a co-effect of environmental conditions and primary magnetic components of the sediments on the diagenesis. Alternating zones of high and low magnetic parameters are observed in Unit A2 of Core YSDP103, which is presumably due to periodic changes of the concentration and/or grain size of magnetic minerals carried into the study area. Cores SSDP-102 and SSDP-103, two studied sediment cores from the Korea Strait contain mud sequences (14 m and 32.62 m in thickness) that were deposited during the last 6,000 years. Analyses of grain sizes and geochemistry of the cores have demonstrated that the sediments have uniform lithology and geochemical properties, however, marked down-core changes in magnetic properties suggest that diagenesis has significantly impacted the magnetic properties. An expanded view of early diagenetic reactions that affect magnetic mineral assemblages is evident in these rapidly deposited continental shelf sediments compared to deep-sea sediments. The studied sediments can be divided into four descending intervals, based on magnetic property variations. Interval 1 is least affected by diagenesis and has the highest concentrations of detrital magnetite and hematite, and the lowest solid-phase sulfur contents. Interval 2 is characterized by the presence of paramagnetic pyrite and sharply decreasing magnetite and hematite concentrations, which suggest active reductive dissolution of detrital magnetic minerals, the absolute minimum abundance of magnetite is reached at the end of this interval. Interval 3 is marked by a progressive loss of hematite with depth, and at the base of this interval, 82% to 88% of the hematite component was depleted and the bulk magnetic mineral concentration was reduced to the lowest value in the entire studied mud section. Interval 4 has an increasing down-core enhancement of authigenic greigite, which is interpreted to have formed due to arrested pyritization resulting from consumption of pore water sulfate with depth. This is the first clear demonstration from an active depositional environment for a delay of thousands of years for acquisition of a magnetization carried by greigite. This detailed view of diagenetic processes in continental shelf sediments suggests that studies of geomagnetic field behavior from such sediments should be conducted with care. Paleoceanographic and paleoclimatic studies based on the magnetic properties of shelf sediments with high sedimentation rates like those in the Korea Strait are also unlikely to provide a meaningful signature associated with syn-depositional environmental processes. The rock magnetic properties of the surface sediments from the fine-grained depositional area on the outer shelf of the East China Sea, an area surrounded by sands, are investigated with a view to providing information on the sediment provenance. Multiple magnetic parameters such as magnetic susceptibility (%), anhysteretic remanent magnetization (ARM), saturation rernanent magnetization (SIRM), coercivities of SIRM (Her), and S ratios (relative abundance of low-coercivity magnetic minerals) are measured for all 179 surface samples, and partial representative samples are examined for their magnetic hysteresis parameters, temperature-dependence of magnetic susceptibility and x-ray diffraction spectra. Our research indicates that the magnetic mineralogy is dominated by magnetite with a small amount of hematite and is primarily of pseudo-single domain (PSD) to multidomain (MD) nature with a detrital origin. In the surface sediments, the granulometry of magnetic fractions is basically independent of grain sizes of the sediment containing the magnetic grains, and the composition of magnetic minerals remains almost homogeneous, that is, with a relatively constant ratio of low to high coercivity fraction throughout the area. The magnetic concentration in the study area generally decreases to the east or southeast accompanied by magnetic-particle fining to the east or to the northeast. The geographic pattern of magnetic properties is most reasonably explained by a major source of sediment jointly from the erosion of the old Huanghe River deposit and the discharge of the Changjiang River. The rock magnetic data facilitate understanding of the transport mechanism of fine-grained sediments in the outer shelf of the East China Sea.

元谋地区晚新生代河湖相沉积物磁性地层学及古环境研究

Yuanmou area lies on the southeastern edge of the Tibetan Plateau, the middlesegment of Yunnan-Sichuan North-South Extending Tectonic Belt and the upper reaches of the Yangztze River, which is renowned for its thick late Cenozoic fluvial-lacustrine sequences that yield rich mammalian fossils including hominoid and early human. The sediments provides great potentials for understanding the relationships between uplift of the Tibetan Plateau, evolution of hominoids and other mammalian and evolution and formation of basins in Southwest China since late Miocene. However, lithostratigrphic and chronologic views on them are controversial and hinder further discussion of the relationships of them. To this end, we selected the Baozidongqing section and the Dapoqing section to carry out systematic lithostratigraphic, magnetostratlgraphic and environmental magnetism researches in this area.The Baozidongqing section was dated to about 10.95-7.17 Ma. The age estimation of the topmost hominoid-bearing layer was about 7.43-7.17 Ma. Rock magnetic results show that the dominant magnetic carrier is hematite, with minor amount of magnetite. Both the composition and concentration of magnetic minerals strongly correlate with the lithostratigraphy, indicating that Yuanmou basin is characterized by alternating of long-term torrid-humid climate and short-term dry-hot climate. But the pattern of these short-term hot-dry events, including both the lasting time and the frequency of their occurrence dramatically changed since -8.1 Ma. Our results infer that the drying process of the Asian west interior and a significant uplift of the Tibetan Plateau would have probably caused jointly the extinction of hominoids, or the emigration of hominoids from Yuanmou to adjacent relatively torrid-humid areas.The strata between the upper of the Dapoqing section, the Niujianbao Hill and Shangnabang area can be linked by three mark layers of conglomerate, which is rather continuous and coherent than physical disturbance by new tectonic activities. Rock magnetic studies indicate that hematite is the main magnetic carriers. The section is dated back to about 2.8-1.37 Ma. Its paleocurrent flowed northeastward, which was a close and stagnant river and swamp environment about 2.2 Ma ago. Then it ran northwestward and turned into an open overflown and braid river sedimentary face during 2.2 to 1.57 Ma. Since 1.57 Ma, the paleocurrent flowed intensely northwestern and about 1.37 Ma ago, it ended the basically continuous fluvial-lacustrine deposition.

汉字形音义加工时间进程的事件相关电位（ERP）研究

The time-courses of orthographic, phonological and semantic processing of Chinese characters were investigated systematically with multi-channel event-related potentials (ERPs). New evidences concerning whether phonology or semantics is processed first and whether phonology mediates semantic access were obtained, supporting and developing the new concept of repetition, overlapping, and alternating processing in Chinese character recognition. Statistic parameter mapping based on physiological double dissociation has been developed. Seven experiments were conducted: I) deciding which type of structure, left-right or non-left-right, the character displayed on the screen was; 2) deciding whether or not there was a vowel/a/in the pronunciation of the character; 3) deciding which classification, natural object or non-natural object, the character was; 4) deciding which color, red or green, the character was; 5) deciding which color, red or green, the non-character was; 6) fixing on the non-character; 7) fixing on the crosslet. The main results are: 1. N240 and P240:N240 and P240 localized at occipital and prefrontal respectively were found in experiments 1, 2, 3, and 4, but not in experiments 5, 6, or 7. The difference between the former 4 and the latter 3 experiments was only their stimuli: the former's were true Chinese characters while the latter's were non-characters or crosslet. Thus Chinese characters were related to these two components, which reflected unique processing of Chinese characters peaking at about 240 msec. 2. Basic visual feature analysis: In comparison with experiment 7 there was a common cognitive process in experiments 1, 2, 4, and 6 - basic visual feature analysis. The corresponding ERP amplitude increase in most sites started from about 60 msec. 3. Orthography: The ERP differences located at the main processing area of orthography (occipital) between experiments 1, 2, 3, 4 and experiment 5 started from about 130 msec. This was the category difference between Chinese characters and non-characters, which revealed that orthographic processing started from about 130 msec. The ERP differences between the experiments 1, 2, 3 and the experiment 4 occurred in 210-250, 230-240, and 190-250 msec respectively, suggesting orthography was processed again. These were the differences between language and non-language tasks, which revealed a higher level processing than that in the above mentioned 130 msec. All the phenomena imply that the orthographic processing does not finished in one time of processing; the second time of processing is not a simple repetition, but a higher level one. 4. Phonology: The ERPs of experiment 2 (phonological task) were significantly stronger than those of experiment 3 (semantic task) at the main processing areas of phonology (temporal and left prefrontal) starting from about 270 msec, which revealed phonologic processing. The ERP differences at left frontal between experiment 2 and experiment 1 (orthographic task) started from about 250 msec. When comparing phonological task with experiment 4 (character color decision), the ERP differences at left temporal and prefrontal started from about 220 msec. Thus phonological processing may start before 220 msec. 5. Semantic: The ERPs of experiment 3 (semantic task) were significantly stronger than those of experiment 2 (phonological task) at the main processing areas of semantics (parietal and occipital) starting from about 290 msec, which revealed semantic processing. The ERP differences at these areas between experiment 3 and experiment 4 (character color decision) started from about 270 msec. The ERP differences between experiment 3 and experiment 1 (orthographic task) started from about 260 msec. Thus semantic processing may start before 260 msec. 6. Overlapping of phonological and semantic processing: From about 270 to 350 msec, the ERPs of experiment 2 (phonological task) were significantly larger than those of experiment 3 (semantic task) at the main processing areas of phonology (temporal and left prefrontal); while from about 290-360 msec, the ERPs of experiment 3 were significantly larger than those of experiment 2 at the main processing areas of semantics (frontal, parietal, and occipital). Thus phonological processing may start earlier than semantic and their time-courses may alternate, which reveals parallel processing. 7. Semantic processing needs part phonology: When experiment 1 (orthographic task) served as baseline, the ERPs of experiment 2 and 3 (phonological and semantic tasks) significantly increased at the main processing areas of phonology (left temporal and frontal) starting from about 250 msec. The ERPs of experiment 3, besides, increased significantly at the main processing areas of semantics (parietal and frontal) starting from about 260 msec. When experiment 4 (character color decision) served as baseline, the ERPs of experiment 2 and 3 significantly increased at phonological areas (left temporal and frontal) starting from about 220 msec. The ERPs of experiment 3, similarly, increased significantly at semantic areas (parietal and frontal) starting from about270 msec. Hence, before semantic processing, a part of phonological information may be required. The conclusion could be got from above results in the present experimental conditions: 1. The basic visual feature processing starts from about 60 msec; 2. Orthographic processing starts from about 130 msec, and repeats at about 240 msec. The second processing is not simple repetition of the first one, but a higher level processing; 3. Phonological processing begins earlier than semantic, and their time-courses overlap; 4. Before semantic processing, a part of phonological information may be required; 5. The repetition, overlapping, and alternating of the orthographic, phonological and semantic processing of Chinese characters could exist in cognition. Thus the problem of whether phonology mediates semantics access is not a simple, but a complicated issue.

On the Physiology of Bistable Percepts

Binocular rivalry refers to the alternating perceptions experienced when two dissimilar patterns are stereoscopically viewed. To study the neural mechanism that underlies such competitive interactions, single cells were recorded in the visual areas V1, V2, and V4, while monkeys reported the perceived orientation of rivaling sinusoidal grating patterns. A number of neurons in all areas showed alternating periods of excitation and inhibition that correlated with the perceptual dominance and suppression of the cell"s preferred orientation. The remaining population of cells were not influenced by whether or not the optimal stimulus orientation was perceptually suppressed. Response modulation during rivalry was not correlated with cell attributes such as monocularity, binocularity, or disparity tuning. These results suggest that the awareness of a visual pattern during binocular rivalry arises through interactions between neurons at different levels of visual pathways, and that the site of suppression is unlikely to correspond to a particular visual area, as often hypothesized on the basis of psychophysical observations. The cell-types of modulating neurons and their overwhelming preponderance in higher rather than in early visual areas also suggests -- together with earlier psychophysical evidence -- the possibility of a common mechanism underlying rivalry as well as other bistable percepts, such as those experienced with ambiguous figures.

Preparation and a time-resolved fluoroimmunoassay application of new europium fluorescent nanoparticles

New silica-based europium fluorescent nanoparticles having surface amino groups were prepared by a covalent binding-copolymerization technique. In the nanoparticles, the fluorescent Eu3+ chelate molecules were covalently bound to silicon atoms to protect the nanoparticles from dye leaking in bio-applications. The amino groups on the surface of nanoparticles made the surface modification and bioconjugation of nanoparticles easier. The nanoparticles were characterized and developed as a new type of fluorescence probe for a highly sensitive time-resolved fluoroimmunoassay (TR-FIA) of human hepatitis B surface antigen (HBsAg).

Theory of Handwriting

Handwriting production is viewed as a constrained modulation of an underlying oscillatory process. Coupled oscillations in horizontal and vertical directions produce letter forms, and when superimposed on a rightward constant velocity horizontal sweep result in spatially separated letters. Modulation of the vertical oscillation is responsible for control of letter height, either through altering the frequency or altering the acceleration amplitude. Modulation of the horizontal oscillation is responsible for control of corner shape through altering phase or amplitude. The vertical velocity zero crossing in the velocity space diagram is important from the standpoint of control. Changing the horizontal velocity value at this zero crossing controls corner shape, and such changes can be effected through modifying the horizontal oscillation amplitude and phase. Changing the slope at this zero crossing controls writing slant; this slope depends on the horizontal and vertical velocity zero amplitudes and on the relative phase difference. Letter height modulation is also best applied at the vertical velocity zero crossing to preserve an even baseline. The corner shape and slant constraints completely determine the amplitude and phase relations between the two oscillations. Under these constraints interletter separation is not an independent parameter. This theory applies generally to a number of acceleration oscillation patterns such as sinusoidal, rectangular and trapezoidal oscillations. The oscillation theory also provides an explanation for how handwriting might degenerate with speed. An implementation of the theory in the context of the spring muscle model is developed. Here sinusoidal oscillations arise from a purely mechanical sources; orthogonal antagonistic spring pairs generate particular cycloids depending on the initial conditions. Modulating between cycloids can be achieved by changing the spring zero settings at the appropriate times. Frequency can be modulated either by shifting between coactivation and alternating activation of the antagonistic springs or by presuming variable spring constant springs. An acceleration and position measuring apparatus was developed for measurements of human handwriting. Measurements of human writing are consistent with the oscillation theory. It is shown that the minimum energy movement for the spring muscle is bang-coast-bang. For certain parameter values a singular arc solution can be shown to be minimizing. Experimental measurements however indicate that handwriting is not a minimum energy movement.

Boundary layer growth and advection of heat over snow and soil patches: Modelling and parametrization

Essery, RLH, RJ Granger and JW Pomeroy, 2006. Boundary layer growth and advection of heat over snow and soil patches: Modelling and parametrization. Hydrological Processes, 20, 953 - 967.

Dynamical Instabilities and Deterministic Chaos in Ballistic Electron Motion in Semiconductor Superlattices

We consider the motion of ballistic electrons within a superlattice miniband under the influence of an alternating electric field. We show that the interaction of electrons with the self-consistent electromagnetic field generated by the electron current may lead to the transition from regular to chaotic dynamics. We estimate the conditions for the experimental observation of this deterministic chaos and discuss the similarities of the superlattice system with the other condensed matter and quantum optical systems.

How Does Binocular Rivalry Emerge from Cortical Mechanisms of 3D Vision?

Under natural viewing conditions, a single depthful percept of the world is consciously seen. When dissimilar images are presented to corresponding regions of the two eyes, binocular rivalyr may occur, during which the brain consciously perceives alternating percepts through time. How do the same brain mechanisms that generate a single depthful percept of the world also cause perceptual bistability, notably binocular rivalry? What properties of brain representations correspond to consciously seen percepts? A laminar cortical model of how cortical areas V1, V2, and V4 generate depthful percepts is developed to explain and quantitatively simulate binocualr rivalry data. The model proposes how mechanisms of cortical developement, perceptual grouping, and figure-ground perception lead to signle and rivalrous percepts. Quantitative model simulations include influences of contrast changes that are synchronized with switches in the dominant eye percept, gamma distribution of dominant phase durations, piecemeal percepts, and coexistence of eye-based and stimulus-based rivalry. The model also quantitatively explains data about multiple brain regions involved in rivalry, effects of object attention on switching between superimposed transparent surfaces, and monocular rivalry. These data explanations are linked to brain mechanisms that assure non-rivalrous conscious percepts. To our knowledge, no existing model can explain all of these phenomena.

Does Binocular Rivalry Emerge from Cortical Mechanisms of 3D Vision

Under natural viewing conditions, a single depthful percept of the world is consciously seen. When dissimilar images are presented to corresponding regions of the two eyes, binocular rivalry may occur, during which the brain consciously perceives alternating percepts through time. Perceptual bistability can also occur in response to a single ambiguous figure. These percepts raise basic questions: What brain mechanisms generate a single depthful percept of the world? How do the same mechanisms cause perceptual bistability, notably binocular rivalry? What properties of brain representations correspond to consciously seen percepts? How do the dynamics of the layered circuits of visual cortex generate single and bistable percepts? A laminar cortical model of how cortical areas V1, V2, and V4 generate depthful percepts is developed to explain and quantitatively simulate binocular rivalry data. The model proposes how mechanisms of cortical development, perceptual grouping, and figure-ground perception lead to single and rivalrous percepts.

Optimisation of Smart Grid performance using centralised and distributed control techniques

A massive change is currently taking place in the manner in which power networks are operated. Traditionally, power networks consisted of large power stations which were controlled from centralised locations. The trend in modern power networks is for generated power to be produced by a diverse array of energy sources which are spread over a large geographical area. As a result, controlling these systems from a centralised controller is impractical. Thus, future power networks will be controlled by a large number of intelligent distributed controllers which must work together to coordinate their actions. The term Smart Grid is the umbrella term used to denote this combination of power systems, artificial intelligence, and communications engineering. This thesis focuses on the application of optimal control techniques to Smart Grids with a focus in particular on iterative distributed MPC. A novel convergence and stability proof for iterative distributed MPC based on the Alternating Direction Method of Multipliers is derived. Distributed and centralised MPC, and an optimised PID controllers' performance are then compared when applied to a highly interconnected, nonlinear, MIMO testbed based on a part of the Nordic power grid. Finally, a novel tuning algorithm is proposed for iterative distributed MPC which simultaneously optimises both the closed loop performance and the communication overhead associated with the desired control.

Multi-scale modelling of atomic layer deposition

High-permittivity ("high-k") dielectric materials are used in the transistor gate stack in integrated circuits. As the thickness of silicon oxide dielectric reduces below 2 nm with continued downscaling, the leakage current because of tunnelling increases, leading to high power consumption and reduced device reliability. Hence, research concentrates on finding materials with high dielectric constant that can be easily integrated into a manufacturing process and show the desired properties as a thin film. Atomic layer deposition (ALD) is used practically to deposit high-k materials like HfO2, ZrO2, and Al2O3 as gate oxides. ALD is a technique for producing conformal layers of material with nanometer-scale thickness, used commercially in non-planar electronics and increasingly in other areas of science and technology. ALD is a type of chemical vapor deposition that depends on self-limiting surface chemistry. In ALD, gaseous precursors are allowed individually into the reactor chamber in alternating pulses. Between each pulse, inert gas is admitted to prevent gas phase reactions. This thesis provides a profound understanding of the ALD of oxides such as HfO2, showing how the chemistry affects the properties of the deposited film. Using multi-scale modelling of ALD, the kinetics of reactions at the growing surface is connected to experimental data. In this thesis, we use density functional theory (DFT) method to simulate more realistic models for the growth of HfO2 from Hf(N(CH3)2)4/H2O and HfCl4/H2O and for Al2O3 from Al(CH3)3/H2O.Three major breakthroughs are discovered. First, a new reaction pathway, ’multiple proton diffusion’, is proposed for the growth of HfO2 from Hf(N(CH3)2)4/H2O.1 As a second major breakthrough, a ’cooperative’ action between adsorbed precursors is shown to play an important role in ALD. By this we mean that previously-inert fragments can become reactive once sufficient molecules adsorb in their neighbourhood during either precursor pulse. As a third breakthrough, the ALD of HfO2 from Hf(N(CH3)2)4 and H2O is implemented for the first time into 3D on-lattice kinetic Monte-Carlo (KMC).2 In this integrated approach (DFT+KMC), retaining the accuracy of the atomistic model in the higher-scale model leads to remarkable breakthroughs in our understanding. The resulting atomistic model allows direct comparison with experimental techniques such as X-ray photoelectron spectroscopy and quartz crystal microbalance.

Quaternary geology of the German North Sea and Western Irish Sea Mud Belt: revised stratigraphies, geotechnical properties, sedimentology and anthropogenic impacts

The climatic development of the Mid to Late Quaternary (last 400,000 years) is characterised by fluctuation between glacial and interglacial periods leading to the present interglacial, the Holocene. In comparison to preceding periods it was believed the Holocene represented a time of relative climatic stability. However, recent work has shown that the Holocene can be divided into cooler periods such as the Little Ice Age alternating with time intervals where climatic conditions ameliorated i.e. Medieval Warm Period, Holocene Thermal Optimum and the present Modern Optimum. In addition, the Holocene is recognised as a period with increasing anthropogenic influence on the environment. Onshore records recording glacial/interglacial cycles as well as anthropogenic effects are limited. However, sites of sediment accumulation on the shallow continental shelf offer the potential to reconstruct these events. Such sites include tunnel valleys and low energy, depositional settings. In this study we interrogated the sediment stratigraphy at such sites in the North Sea and Irish Sea using traditional techniques, as well as novel applications of geotechnical data, to reconstruct the palaeoenvironmental record. Within the German North Sea sector a combination of core, seismic and in-situ Cone Penetration Testing (CPT) data was used to identify sedimentary units, place them within a morphological context, relate them to glacial or interglacial periods stratigraphically, and correlate them across the German North Sea. Subsequently, we were able to revise the Mid to Late Quaternary stratigraphy for the North Sea using this new and novel data. Similarly, Holocene environmental changes were investigated within the Irish Sea at a depositional site with active anthropogenic influence. The methods used included analyses on grain-size distribution, foraminifera, gamma spectrometry, AMS 14C and physical core logging. The investigation revealed a strong fluctuating climatic signal early in the areas history before anthropogenic influence affects the record through trawling.

Pulsating tandem microbubble for localized and directional single-cell membrane poration.

The interaction of laser-generated tandem microbubble (maximum diameter of about 50  μm) with single (rat mammary carcinoma) cells is investigated in a 25-μm liquid layer. Antiphase and coupled oscillation of the tandem microbubble leads to the formation of alternating, directional microjets (with max microstreaming velocity of 10  m/s) and vortices (max vorticity of 350 000  s{-1}) in opposite directions. Localized and directional membrane poration (200 nm to 2  μm in pore size) can be produced by the tandem microbubble in an orientation and proximity-dependent manner, which is absent from a single oscillating microbubble of comparable size and at the same stand-off distance.