复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

东营凹陷下第三系岩性油气藏形成条件与预测

The Dongying depression, located in the northern part of the jiyang Sag in the Buohaiwan Basin, comprises one of the major oil-producing bases of the Shengli oil-field. The prediction and exploration of subtle or litho1ogical oil traps in the oil-field has become the major confronted target. This is also one of the frontier study areas in the highly-explored oil-bearing basins in East China and abroad. Based on the integrated analysis of the geological, seismic and logging data and the theories of sequence stratigraphy, tectono-stratigraphy and petroleum system, the paper has attempted to document the characteristics of the sequence stratigraphic and structural frameworks of the low Tertiary, the syndepositional faults and their control on deposition, and then to investigate the forming conditions and distribution of the tithological oil traps in the depression. The study has set up a set of analysis methods, which can be used to effectively analysis the sequence stratigraphy of inland basins and predict the distribution of sandstone reservoirs in the basins. The major achievements of the study are as follows: 1. The low Tertiary can be divided into 4 second-order sequences and 13 third-order sequences, and the systems tracts in the third-order sequences have been also identified based on the examination and correction of well logging data and seismic profiles. At the same time, the parasequences and their stacking pattern in the deltaic systems of the third member of the Shahejie Formation have been recognized in the key study area. It has been documented that the genetic relation of different order sequences to tectonic, climatic and sediment supply changes. The study suggested that the formation of the second-order sequences was related to multiple rifting, while the activity of the syndepositional faults controlled the stacking pattern of parasequences of the axial deltaic system in the depression. 2. A number of depositional facies have been recognized in the low Tertiary on the basis of seismic facies and well logging analysis. They include alluvial fan, fan delta or braided delta, axial delta, lowstand fan, lacustrine and gravity flow deposits. The lacustrine lowstand fan deposits are firstly recognized in the depression, and their facies architecture and distribution have been investigated. The study has shown that the lowstand fan deposits are the important sandstone reservoirs as lithological oil traps in the depression. 3. The mapping of depositional systems within sequences has revealed the time and special distrbution of depositional systems developed in the basin. It is pointed out that major elastic systems comprise the northern marginal depositional systems consisting of alluvial fan, fan delta and offshore lowstand fan deposits, the southern gentle slope elastic deposits composed of shallow lacustrine, braided delta and lowstand fan deposits and the axial deltaic systems including those from eastern and western ends of the depression. 4. The genetic relationship between the syndepositional faults and the distribution of sandstones has been studied in the paper, upper on the analysis of structural framework and syndepositional fault systems in the depression. The concept of structural slope-break has been firstly introduced into the study and the role of syndepositional faults controlling the development of sequence architecture and distribution of sandstones along the hinged and faulted margins have been widely investigated. It is suggested that structural styles of the structural slope-break controlled the distribution of lowstand fan deposits and formed a favorable zone for the formation of lithological or structure-lithological oil traps in the basin. 5. The paper has made a deep investigation into the forming condition and processes of the lithological traps in the depression, based the analysis of composition of reservoir, seal and resource rocks. It is pointed out that there were two major oil pool-forming periods, namely the end of the Dongying and Guangtao periods, and the later one is the most important. 6. The study has finally predicted a number of favorable targets for exploration of lithologieal traps in the depression. Most of them have been drilled and made great succeed with new discovered thousands tons of raw oil reserves.

济阳地区深部结构与含油气性讨论

Oil and gas exploration of marine strata in China's Pre-Cenozoic residual basins is regarded as a worldwide puzzle because of existent problems and cruxes. Objectively speaking, the subsurface geologic structure is complicated, and the surface conditions of some areas are tough. On the other hand, there are still many problems to be solved in oil and gas exploration technologies of Pre-Cenozoic marine fades, and theoretic cognition about petroleum geology is not profound yet. Therefore, it is principal to explore integrated geophysical research ways of Pre-Cenozoic residual basins. Seismic prospecting and geophysical integrated interpretation technologies aimed at middle Paleozoic marine facies with deeper burial and complicated geologic conditions have not formed due to bad quality of deep strata data. Pre-Cenozoic strata, and especially extension, thickness and internal structure of Paleozoic strata can not be recognized from seismic profiles, so it is hard to systematically cognize structural features and oil-gas resources prospect of Pre-Cenozoic basins. To further investigation of fabric and structural features, basin prototype, formation and evolution pattern of Pre-Cenozoic basins, and also their control over formation, migration and aggregation of oil and gas, will play a guiding and promotive role in developing new surveying areas, selecting advantageous zones and predicting oil-gas resources.This paper follows the modem macrocontrol theory of "Region constrains local, deep strata controls shallow ones", and uses the integrated geophysical method of "One guide, two hinges, three combinations and multi feedbacks'*. Based on several years of geological and geophysical results of the Shengli Oilfield, and 14 newly-joint regional seismic profiles, deep structure and oil-gas bearing capacity of the Jiyang area are discussed and new cognitions are drawn as below.Seismic identification marks Tr, Tg, Tgl and Tg2 are established for importantPre-Cenozoic geological interfaces, and promoted to the whole Jiyang area.Through area-wide tracking and clogging of important seismic reflection marker,the isochronic framework of pre-Tertiary basin is set up in the Jiyang area for the firsttime, which is vital for basin research.Integrated with geological and geophysical research results, the Jiyang area isdivided into four first-order tectonic sequences- basement, lower tectonic layer,upper tectonic layer, and top tectonic layer. The basement and lower tectonic sequence which are related to Pre-Cenozoic are studied with emphasis.Through the research of regional seismic profiles, the point of view is given thatthe Kongdian Formation of Jiyang is structural transition period. The top-bottomunconformable interface of the Kongdian Formation is found out for the first time,and the basin model is determined primarily, which lay a basis for prototype basinresearch of the Jiyang Kongdian Formation.The distribution status of Middle-Paleozoic is delineated in the Jiyang area.The maximum thickness of Paleozoic lies in the top of the south declivity of half-graben. The thickness gets thinner towards the center of Mesozoic and Cenozoic half-graben basin, and even disappears. Structural action in the west-north affects the distribution of Paleozoic residual strata.6. The features of second-order tectonic sequence of the Jiyang depression isstudied and its evolution history of is rebuilt.Combined with the 5-stage evolution history of the China continent and structure evolution features of the Jiyang area, the structure sedimentary process since Paleozoic is divided into 5 periods - basement forming , Indosinian orogenic, Yanshan negative reversal, Himalayan extension and Neogene subsidence period.Combined with the research results of gravity, magnetic surveying and regionalprofiles, this paper brings forward the idea for the first time that the western boundaryof the Jiyang depression is the Ningjin-Yangpan fracture zone, and forms aside-column assemblage with the Wudi fracture zone.The opinion that under Middle-Cenozoic basins in the middle Jiyang area theremight superimpose an old residual basin is given for the first time. And if it is provedto be true, a new exploration space will be pioneered for Jiyang and even north China.There exists many types of tectonic-stratigraphic traps formed under piezotropy,extension and compound action in Pre-Cenozoic Jiyang. Therein all kinds of burialhills are the most important oil-gas trap type of Pre-Cenozoic, which should besurveyed layeredly according to the layout of oil sources.As such a new challenging project and field, the paper systematically analyses different geophysical responses of the Jiyang area, frames the deep structure of the area, and preliminarily recognizes the Pre-Cenozoic residual basins. It breaks through to a certain extent in both theory and practice, and is expected to provide new geophysical and geotectonic clues for deep exploration in Shengli.

Pump-dump control and the related transient absorption spectroscopies

We combine theories of optimal pump-dump control and the related transient probe absorption spectroscopy in order to elucidate the relation between these two optical processes and the possibility of experimental realization. In the weak response regime, we identify the globally optimal pair of pump-dump control fields, and further propose a second-order difference detection scheme to monitor the wave packets dynamics that is jointly controlled by both the pump and dump fields. The globally optimal solution serves also as the initial input for the iterative search for the optimal control fields in the strong response regime. We use a model I-2 molecule to demonstrate numerically the pump-dump control and the detection of a highly vibrationally excited wave packet focusing dynamics on the ground X surface in both the weak and strong response regimes. The I2B surface serves as the intermediate to assist the pump-dump control and the optical detection processes. Demonstrated in the strong response regime are the optimal pair of pump-dump molecular-pi pulses that invert nearly total population onto the predefined target region within a half period of vibration motion. (C) 1999 American Institute of Physics. [S0021-9606(99)00115-4].

A framework for the transformation of the creative industries in a digital age.

The creative industries sector faces a constantly changing context characterised by the speed of the development and deployment of digital information systems and Information Communications Technologies (ICT) on a global scale. This continuous digital disruption has had significant impact on the whole value chain of the sector: creation and production; discovery and distribution; and consumption of cultural goods and services. As a result, creative enterprises must evolve business and operational models and practices to be sustainable. Enterprises of all scales, type, and operational model are affected, and all sectors face ongoing digital disruption. Management consultancy practitioners and business strategy academics have called for new strategy development frameworks and toolkits, fit for a continuously changing world. This thesis investigates a novel approach to organisational change appropriate to the digital age, in the context of the creative sector in Scotland. A set of concepts, methods, tools, and processes to generate theoretical learning and practical knowing was created to support enterprises to digitally adapt through undertaking journeys of change and organisational development. The framework is called The AmbITion Approach. It was developed by blending participatory action research (PAR) methods and modern management consultancy, design, and creative practices. Empirical work also introduced to the framework Coghlan and Rashford’s change categories. These enabled the definition and description of the extent to which organisations developed: whether they experienced first order (change), second order (adaptation) or third order (transformation) change. Digital research tools for inquiry were tested by a pilot study, and then embedded in a longitudinal study over two years of twentyone participant organisations from Scotland’s creative sector. The author applied and investigated the novel approach in a national digital development programme for Scotland’s creative industries. The programme was designed and delivered by the author and ran nationally between 2012-14. Detailed grounded thematic analysis of the data corpus was undertaken, along with analysis of rich media case studies produced by the organisations about their change journeys. The results of studies on participants, and validation criteria applied to the results, demonstrated that the framework triggers second (adaptation) and third order change (transformation) in creative industry enterprises. The AmbITion Approach framework is suitable for the continuing landscape of digital disruption within the creative sector. The thesis contributes to practice the concepts, methods, tools, and processes of The AmbITion Approach, which have been empirically tested in the field, and validated as a new framework for business transformation in a digital age. The thesis contributes to knowledge a theoretical and conceptual framework with a specific set of constructs and criteria that define first, second, and third order change in creative enterprises, and a robust research and action framework for the analysis of the quality, validity and change achieved by action research based development programmes. The thesis additionally contributes to the practice of research, adding to our understanding of the value of PAR and design thinking approaches and creative practices as methods for change.

Optimization of physical and optical properties of biodegradable edible films based on pea starch and guar gum

The influence of process variables (pea starch, guar gum and glycerol) on the viscosity (V), solubility (SOL), moisture content (MC), transparency (TR), Hunter parameters (L, a, and b), total color difference (ΔE), yellowness index (YI), and whiteness index (WI) of the pea starch based edible films was studied using three factors with three level Box–Behnken response surface design. The individual linear effect of pea starch, guar and glycerol was significant (p < 0.05) on all the responses. However, a value was only significantly (p < 0.05) affected by pea starch and guar gum in a positive and negative linear term, respectively. The effect of interaction of starch × glycerol was also significant (p < 0.05) on TR of edible films. Interaction between independent variables starch × guar gum had a significant impact on the b and YI values. The quadratic regression coefficient of pea starch showed a significant effect (p < 0.05) on V, MC, L, b, ΔE, YI, and WI; glycerol level on ΔE and WI; and guar gum on ΔE and SOL value. The results were analyzed by Pareto analysis of variance (ANOVA) and the second order polynomial models were developed from the experimental design with reliable and satisfactory fit with the corresponding experimental data and high coefficient of determination (R2) values (>0.93). Three-dimensional response surface plots were established to investigate the relationship between process variables and the responses. The optimized conditions with the goal of maximizing TR and minimizing SOL, YI and MC were 2.5 g pea starch, 25% glycerol and 0.3 g guar gum. Results revealed that pea starch/guar gum edible films with appropriate physical and optical characteristics can be effectively produced and successfully applied in the food packaging industry.

Skin Color-Based Video Segmentation under Time-Varying Illumination

A novel approach for real-time skin segmentation in video sequences is described. The approach enables reliable skin segmentation despite wide variation in illumination during tracking. An explicit second order Markov model is used to predict evolution of the skin-color (HSV) histogram over time. Histograms are dynamically updated based on feedback from the current segmentation and predictions of the Markov model. The evolution of the skin-color distribution at each frame is parameterized by translation, scaling and rotation in color space. Consequent changes in geometric parameterization of the distribution are propagated by warping and resampling the histogram. The parameters of the discrete-time dynamic Markov model are estimated using Maximum Likelihood Estimation, and also evolve over time. The accuracy of the new dynamic skin color segmentation algorithm is compared to that obtained via a static color model. Segmentation accuracy is evaluated using labeled ground-truth video sequences taken from staged experiments and popular movies. An overall increase in segmentation accuracy of up to 24% is observed in 17 out of 21 test sequences. In all but one case the skin-color classification rates for our system were higher, with background classification rates comparable to those of the static segmentation.

Beta-Reduction As Unification

We define a unification problem ^UP with the property that, given a pure lambda-term M, we can derive an instance Gamma(M) of ^UP from M such that Gamma(M) has a solution if and only if M is beta-strongly normalizable. There is a type discipline for pure lambda-terms that characterizes beta-strong normalization; this is the system of intersection types (without a "top" type that can be assigned to every lambda-term). In this report, we use a lean version LAMBDA of the usual system of intersection types. Hence, ^UP is also an appropriate unification problem to characterize typability of lambda-terms in LAMBDA. It also follows that ^UP is an undecidable problem, which can in turn be related to semi-unification and second-order unification (both known to be undecidable).

System F with Constraint Types

System F is a type system that can be seen as both a proof system for second-order propositional logic and as a polymorphic programming language. In this work we explore several extensions of System F by types which express subtyping constraints. These systems include terms which represent proofs of subtyping relationships between types. Given a proof that one type is a subtype of another, one may use a coercion term constructor to coerce terms from the first type to the second. The ability to manipulate type constraints as first-class entities gives these systems a lot of expressive power, including the ability to encode generalized algebraic data types and intensional type analysis. The main contributions of this work are in the formulation of constraint types and a proof of strong normalization for an extension of System F with constraint types.

An Information Theoretic Framework for Field Monitoring Using Autonomously Mobile Sensors

We consider a mobile sensor network monitoring a spatio-temporal field. Given limited cache sizes at the sensor nodes, the goal is to develop a distributed cache management algorithm to efficiently answer queries with a known probability distribution over the spatial dimension. First, we propose a novel distributed information theoretic approach in which the nodes locally update their caches based on full knowledge of the space-time distribution of the monitored phenomenon. At each time instant, local decisions are made at the mobile nodes concerning which samples to keep and whether or not a new sample should be acquired at the current location. These decisions account for minimizing an entropic utility function that captures the average amount of uncertainty in queries given the probability distribution of query locations. Second, we propose a different correlation-based technique, which only requires knowledge of the second-order statistics, thus relaxing the stringent constraint of having a priori knowledge of the query distribution, while significantly reducing the computational overhead. It is shown that the proposed approaches considerably improve the average field estimation error by maintaining efficient cache content. It is further shown that the correlation-based technique is robust to model mismatch in case of imperfect knowledge of the underlying generative correlation structure.

Estimation and Prediction of Evolving Color Distributions for Skin Segmentation Under Varying Illumination

A novel approach for real-time skin segmentation in video sequences is described. The approach enables reliable skin segmentation despite wide variation in illumination during tracking. An explicit second order Markov model is used to predict evolution of the skin color (HSV) histogram over time. Histograms are dynamically updated based on feedback from the current segmentation and based on predictions of the Markov model. The evolution of the skin color distribution at each frame is parameterized by translation, scaling and rotation in color space. Consequent changes in geometric parameterization of the distribution are propagated by warping and re-sampling the histogram. The parameters of the discrete-time dynamic Markov model are estimated using Maximum Likelihood Estimation, and also evolve over time. Quantitative evaluation of the method was conducted on labeled ground-truth video sequences taken from popular movies.

Type Inference for Variant Object Types

Existing type systems for object calculi are based on invariant subtyping. Subtyping invariance is required for soundness of static typing in the presence of method overrides, but it is often in the way of the expressive power of the type system. Flexibility of static typing can be recovered in different ways: in first-order systems, by the adoption of object types with variance annotations, in second-order systems by resorting to Self types. Type inference is known to be P-complete for first-order systems of finite and recursive object types, and NP-complete for a restricted version of Self types. The complexity of type inference for systems with variance annotations is yet unknown. This paper presents a new object type system based on the notion of Split types, a form of object types where every method is assigned two types, namely, an update type and a select type. The subtyping relation that arises for Split types is variant and, as a result, subtyping can be performed both in width and in depth. The new type system generalizes all the existing first-order type systems for objects, including systems based on variance annotations. Interestingly, the additional expressive power does not affect the complexity of the type inference problem, as we show by presenting an O(n^3) inference algorithm.

Modeling the Possible Influences of Eye Movements on the Refinement of Cortical Direction Selectivity

The second-order statistics of neural activity was examined in a model of the cat LGN and V1 during free-viewing of natural images. In the model, the specific patterns of thalamocortical activity required for a Bebbian maturation of direction-selective cells in VI were found during the periods of visual fixation, when small eye movements occurred, but not when natural images were examined in the absence of fixational eye movements. In addition, simulations of stroboscopic reming that replicated the abnormal pattern of eye movements observed in kittens chronically exposed to stroboscopic illumination produced results consistent with the reported loss of direction selectivity and preservation of orientation selectivity. These results suggest the involvement of the oculomotor activity of visual fixation in the maturation of cortical direction selectivity.