Bimolecular quadruplexes and their transitions to higher-order molecular structures detected by ESI-FTICR-MS

Four individual quadruplexes, which are self-assembled in ammonium acetate solution from telomeric sequences of closely related DNA strands - d(G(4)T(4)G(4)), d(G(3)T(4)G(4)), d(G(3)T(4)G(3)), and d(G(4)T(4)G(3)) - have been detected in the gas phase using electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS). The bimolecular quadruplexes associate with the same number of NH4+ in the gas phase as NMR shows that they do in solution. The quadruplex structures formed in solution are maintained in the gas phase. Furthermore, the mass spectra show that the bimolecular quadruplexes generated by the strands d(G(3)T(4)G(3)) and d(G(4)T(4)G(3)) are unstable, being converted into trimolecular and tetramolecular structures with increasing concentrations of NH4+ in the solution. Circular dichroism (CD) spectra reveal structural changes during the process of strand stoichiometric transitions, in which the relative orientation of strands in the quadruplexes changes from an antiparallel to a parallel arrangement. Such changes were observed for the strand d(G(4)T(4)G(3)), but not for the strand d(G(3)T(4)G(3)). The present work provides a significant insight into the formation of various DNA quadruplexes, especially the higher-order species.

2-cyanoprop-2-yl dithiobenzoate mediated reversible addition-fragmentation chain transfer polymerization of acrylonitrile targeting a polymer with a higher molecular weight

The reversible addition-fragmentation chain transfer (RAFT) polymerization of acrylonitrile (AN) mediated by 2-cyanoprop-2-yl dithiobenzoate was first applied to synthesize polyacrylonitrile (PAN) with a high molecular weight up to 32,800 and a polydispersity index as low as 1.29. The key to success was ascribed to the optimization of the experimental conditions to increase the fragmentation reaction efficiency of the intermediate radical. In accordance with the atom transfer radical polymerization of AN, ethylene carbonate was also a better solvent candidate for providing higher controlled/living RAFT polymerization behaviors than dimethylformamide and dimethyl sulfoxide. The various experimental parameters, including the temperature, the molar ratio of dithiobenzoate to the initiator, the molar ratio of the monomer to dithiobenzoate, the monomer concentration, and the addition of the comonomer, were varied to improve the control of the molecular weight and polydispersity index. The molecular weights of PANS were validated by gel permeation chromatography along with a universal calibration procedure and intrinsic viscosity measurements. H-1 NMR analysis confirmed the high chain-end functionality of the resultant polymers.

Monte Carlo simulation of the grafting of maleic anhydride onto polypropylene at higher temperature

The graft of maleic anhydride (MAH) onto isotactic polypropylene (iPP) initiated by dicumyl peroxide (DCP) at 190 degreesC was studied by means of the Monte Carlo method. The ceiling temperature theory, i.e., no possibility for the homopolymerization of MA-H to occur at higher temperatures, was used in this study. The simulation results show that most MAH monomers were grafted onto the radical chain ends arising from beta scission at a lower MAH concentration, whereas the amount of MAH monomers attached to the tertiary carbons was much larger than that grafted onto the radical chain ends at a higher MAH concentration for various DCP concentrations. This conclusion gives a good interpretation for the disagreement on the grafting sites along a PP chain. Moreover, it was found that the grafting degree increased considerably up to a peak value; thereafter, it decreased continuously with increasing MA-H concentration. The peak shifted in the lower MAH concentration direction and became lower and lower with increasing DCP concentration. When the DCP concentration was below 0.1 wt %, the peak was hardly observed. Those results are in good agreement with the experiments.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa's results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa's (1999) results, and the applicable scope of water depth is deeper.

Higher order response of nonlinear composite in external AC electric field

The fifth-order effective nonlinear responses at fundament frequency and higher-order harmonics are given for nonlinear composites, which obey a current-field relation of the form J = sigmaE + x\E\(2) E, if a sinusoidal alternating current (AC) external field with finite frequency omega is applied. As two examples, we have investigated the cylinder and spherical inclusion embedded in a host and, for larger volume fraction, also derived the formulae of effective nonlinear responses at higher-order harmonics by the aid of the general effective response definition. Furthermore, the relationships between effective nonlinear responses at harmonics are given. (C) 2003 Elsevier Science B.V. All rights reserved.

NONLINEAR GRAVITY-WAVES IN DEEP-WATER

In consideration of the problem on the boundary condition of nonlinear free water wave, coordinate transform is used to handle the free boundary. Supposing the solution form be the traveling wave, the ordinary differential equations of the one-order autonomous system with two variables are caused, then expanding the nonlinear terms at the equilibrium point with the Taylor expansion, we obtained the solution to traveling wave. The linear approximate equation near the equilibrium point is the small amplitude wave. A new nonlinear periodic traveling wave and nonlinear dispersion relation are shown when expanding to the second-order terms. A conclusion that the expansion of dispersion relation does not contain any odd-power terms of wave steepness and because of the nonlinear effort an oscillate structure is produced in the vertical direction is drawn.

Effective AC and DC responses of nonlinear composites at higher concentration

For higher concentration of inclusions, an effective medium approximation (EMA) method is used to investigate the bulk effective response of weakly nonlinear composites, which are subject to the constitutive relation of electric displacement and electric field, D-alpha = epsilon E-alpha + chi(alpha)|E|(2) E. As an example of three dimensions, under the external AC and DC electric fields E-app = E-a (1 + sin omega t), we have derived the general effective nonlinear response of composites by the EMA method for higher concentration of spherical inclusions. Furthermore, the effective nonlinear responses at harmonics are predicted.

细粒沉积物多组分粒度特征及其成因机制研究

By using high-resolution laser grain size instrument Mastersizer 2000, the grain size distribution of windblown depositions (loess and sandy dunes), aqueous sediments (lake, river, riverside and foreshore sand), weathering crust, sloping materials and other fine-grain sediments are systemically measured. The multimodal characteristics of grain size distribution of these sediments are carefully studied. The standard patterns and their grain size characteristics of various sediments are systemically summarized. The discrepancies of multimodal distribution among windblown depositions, aqueous sediments and other sediments are concluded and the physical mechanisms of grain size multimodal distribution of various sediments are also discussed in this paper. The major conclusions are followed: 1. The multimodal characteristic of grain size distribution is a common feature in all sediments and results from properties of transportation medium, dynamic intensity, transportation manner and other factors. 2. The windblown depositions are controlled by aerodynamic forcing, resulting in that the median size of the predominant mode gradually decreases form sandy dunes to loess. Similarly, the aqueous sediments are impacted by dynamic forces of water currents and the median grain size of the predominant mode decreases gradually from river to lake sediments. Because the kinetic viscidity of air is lower than of water, the grain size of modes of windblown depositions is usually finer than that of corresponding modes of aqueous sediments. Typical characteristics of sediments grain size distribution of various sediments have been summarized in the paper: (1) Suspended particles which diameters are less than 75μm are dominant in loess and dust. There are three modes in loess’ grain size distribution: fine, median and coarse (the median size is <1μm、1-10μm、10-75μm, respectively). The coarse mode which percentage is larger than that of others is controlled by source distance and aerodynamic intensity of dust source areas. Some samples also have a saltation mode which median size is about 300-500μm. Our analysis demonstrates that the interaction of wind, atmospheric turbulence, and dust grain gravity along the dust transportation path results in a multimodal grain size distribution for suspended dust. Changes in the median sizes of the coarse and medium modes are related to variation in aerodynamic forcing (lift force related to vertical wind and turbulence) during dust entrainment in the source area and turbulence intensity in the depositional area. (2) There is a predominant coarse saltation mode in grain size distribution of sandy dunes, which median size is about 100-300μm and the content is larger than that of other modes. The grain size distribution curve is near axis symmetric as a standard logarithm normal function. There are some suspended particles in some samples of sandy dunes, which distribution of the fine part is similar to that of loess. Comparing with sandy samples of river sediments, the sorting property of sandy dunes is better than of river samples although both they are the saltation mode. Thus, the sorting property is a criterion to distinguish dune sands and river sands. (3) There are 5~6 modes (median size are <1μm, 1-10μm, 10-70μm, 70-150μm, 150-400μm, >400μm respectively) in grain size distribution of lacustrine sediments. The former 4 modes are suspensive and others are saltated. Lacustrine sediments can be divided into three types: lake shore facies, transitional facies and central lake facies. The grain size distributions of the three facies are distinctly different and, at the same time, the transition among three modes is also clear. In all these modes, the third mode is a criteria to identify the windblown deposition in the watershed. In lake shore sediments, suspended particles are dominant, a saltation mode sometimes occurs and the fourth mode is the most important mode. In the transitional facies, the percentage of the fourth mode decreases and that of the second mode increases from lake shore to central lake. In the central lake facies, the second mode is dominant. A higher content of the second mode indicates its position more close to the central lake. (4) The grain size distribution of river sediments is the most complex. It consist of suspension, saltation and rolling modes. In most situations, the percentage of the saltation mode is larger than that of other modes. The percentage of suspension modes of river sediments is more than of sandy dunes. The grain size distribution of river sediments indicates dynamic strength of river currents. If the fourth mode is dominant, the dynamic forcing of river is weaker, such as in river floodplain. If the five or sixth mode is dominant, the water dynamic forcing of rivers is strong. (5) Sediments can be changed by later forcing in different degree to form some complicated deposition types. In the paper, the grain size distribution of aqueous sediments of windblown deposition, windblown sediments of aqueous deposition, weathering crust and slope materials are discussed and analyzed. 3. The grain size distribution characteristics of different sediments are concluded: (1) Modal difference: Usually there are suspended and saltation modes in the windblown deposition. The third mode is dominant in loess dust and the fifth mode is predominant in sandy dunes. There are suspended, saltation and rolling particles in aqueous sediments. In lacustrine sediments, the second and fourth mode are predominant for central lake facies or lake shore facies, respectively. In river sediments, the fourth, or fifth, or sixth mode is predominant. Suspended modes: the grain size of suspended particles of windblown depositions usually is less than 75μm. The content of suspended particles is lower or none in sandy dunes. However, suspended particles of aqueous sediments may reach 150μm. Difference in grain size of suspended modes represents difference between transitional mediums and the strength of dynamic forcing. Saltation modes: the median size of saltation mode of sandy dunes fluctuates less than that of river sediments. (2) Loess dust and lacustrine sediment: Their suspended particles are clearly different. There is an obvious pit between the second and the third modes in grain size distribution of lacustrine sediments. The phenomenon doesn’t occur in loess dust. In lacustrine sediments, the second mode can be a dominant mode, such as central lake facies, and contents of the second and the third modes change reversely. However, the percentage of the third mode is always the highest in loess dust. (3) Dune Sand and fluvial sand: In these two depositions, the saltation particles are dominant and the median sizes of their saltation modes overlay in distribution range. The fifth mode of dune sand fluctuates is sorted better than that of fluvial sand. (4) Lacustrine and fluvial sediments: In lacustrine sediments, there are 5-6 modes and suspended particles can be predominant. The second mode is dominant in central lake facies and the third mode is dominant in lake shore facies. Saltation or roll modes occurred in central lake facies may indicate strong precipitation events. In fluvial sediments, saltation particles (or rolling particles) usually dominant. 4. A estimation model of lake depth is firstly established by using contents of the second, the third and the fourth modes. 5. The paleoenvironmental history of the eastern part of SongLiao basin is also discussed by analyzing the grain size distribution of Yushu loess-like sediments in Jilin. It was found that there is a tectonic movement before 40ka B.P. in SongLiao basin. After the movement, loess dust deposited in Yushu area as keerqin desert developed. In recent 2000 years, the climate became drier and more deserts activated in the eastern part of Song-Liao basin.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

内蒙古中部道郎和都格铜多金属矿的地质特征矿床类型及找矿模式研究

Daolangheduge copper polymetallic deposit is located on east edge of Ondor Sum-Bainaimiao metallogenic belt, which is a prospective area of porphyry copper deposit, in Xianghuangqi of central Inner Mongolia. Geotectonically, it occurred in the continental margin accretion belt along the north margin of North China Plate, south of the suture zone between North China Plate and Siberian Plate. The intrusive rocks in this area mainly consist of intermediate-acid magmatic rocks, and the quartz veins, tourmaline veins and the transitional phase are comparatively developed. According to our research, the ore-bearing rock body is mainly quartz diorite while the surrounding rock is mainly biotite granite. Besides, the wall rock alteration are mainly propylitization, pyritization and silicification, which consist of epidotization, actinolitization, chloritzation and so on. The metallic minerals are mainly chalcopyrite and pyrite. In addition, the primary ore is mainly of quartz-chalcopyrite-pyrite type. Above all, Daolangheduge copper polymetallic deposit is suggested to be categorized in the porphyry copper type. With isotopic dating and geochemical research on quartz diorite of ore-bearing rock body, the zircon LA-ICP-MS U-Pb dating of two samples yields an age of 266±2 Ma, falling into the range of late Permian Epoch. It is the first accurate age data in Xianghuangqi area, so it should play a key role in the research of deposit and magmatic rocks in this area. With the major elements and trace elements analysis of 14 samples, the quartz diorite should be among the calc-alkaline series, the geochemical characteristics show higher large-ion lithophile elements of Rb, Sr and LREE, low high-field strength elements of Nb, Ta and high transition elements of Cu, Cr . Also, the REE patterns have negative Eu anomalies. With the same analysis of 4 sample for the biotite granite, the geochemical characteristics show higher Rb, Th,, Zr, Hf and LREE, low Nb, Sm and HREE and Eu has no anomaly. It should be among the calc-alkaline series, over aluminum quality and has characteristics of Adakites. According to isotopic dating and geochemical characteristics of ore-bearing rock body, it is suggested that its materials mainly derived from upper mantle that had fractional crystallization and its magma source region may be affected by fluid metasomatism of paleo-asian ocean. It should be an extensional process of post-orogeny according to regional tectonic evolution. Consequently, because of the decrease of temperature and pressure, the ore forming fluid was raised to surface and mineralized accompanied by magmatic activity which might occur in south of the suture zone. By geological survey, further geophysical and geochemical work is needed. In this area, we have accomplished high precision magnetic prospecting, high density electrical survey, gravity prospecting, soil geochemical prospecting, X-ray fluorescence analyzer prospecting and so on. According to geophysical and geochemical abnormal and surface occurrence, 11 drills are arranged to verification. The type of ores are mainly quartz-chalcopyrite-pyrite ores within 3 drills by drill core logging. Although the grade as well as the scale of already-found Cu deposits are insufficient for industrial exploitation, the mineralization prospect in this region is supposed to be great and the potential in mineral exploration at depth is excellent.

波动方程偏移速度分析

On the subject of oil and gas exploration, migration is an efficacious technique for imagining structures underground. Wave-equation migration (WEM) dominates over other migration methods in accuracy, despite of higher computational cost. However, the advantages of WEM will emerge as the progress of computer technology. WEM is sensitive to velocity model more than others. Small velocity perturbations result in grate divergence in the image pad. Currently, Kirrchhoff method is still very popular in the exploration industry for the reason of difficult to provide precise velocity model. It is very urgent to figure out a way to migration velocity modeling. This dissertation is mainly devoted to migration velocity analysis method for WEM: 1. In this dissertation, we cataloged wave equation prestack depth migration. The concept of migration is introduced. Then, the analysis is applied to different kinds of extrapolate operator to demonstrate their accuracy and applicability. We derived the DSR and SSR migration method and apply both to 2D model. 2. The output of prestack WEM is in form of common image gathers (CIGs). Angle domain common image gathers (ADCIGs) gained by wave equation are proved to be free of artifacts. They are also the most potential candidates for migration velocity analysis. We discussed how to get ADCIGs by DSR and SSR, and obtained ADCIGs before and after imaging separately. The quality of post stack image is affected by CIGs, only the focused or flattened CIGs generate the correct image. Based on wave equation migration, image could be enhanced by special measures. In this dissertation we use both prestack depth residual migration and time shift imaging condition to improve the image quality. 3. Inaccurate velocities lead to errors of imaging depth and curvature of coherent events in CIGs. The ultimate goal of migration velocity analysis (MVA) is to focus scattered event to correct depth and flatten curving event by updating velocities. The kinematic figures are implicitly presented by focus depth aberration and kinetic figure by amplitude. The initial model of Wave-equation migration velocity analysis (WEMVA) is the output of RMO velocity analysis. For integrity of MVA, we review RMO method in this dissertation. The dissertation discusses the general ideal of RMO velocity analysis for flat and dipping events and the corresponding velocity update formula. Migration velocity analysis is a very time consuming work. Respect to computational convenience, we discus how RMO works for synthetic source record migration. In some extremely situation, RMO method fails. Especially in the areas of poorly illuminated or steep structure, it is very difficult to obtain enough angle information for RMO. WEMVA based on wave extrapolate theory, which successfully overcome the drawback of ray based methods. WEMVA inverses residual velocities with residual images. Based on migration regression, we studied the linearized scattering operator and linearized residual image. The key to WEMVA is the linearized residual image. Residual image obtained by Prestack residual migration, which based on DSR is very inefficient. In this dissertation, we proposed obtaining residual migration by time shift image condition, so that, WEMVA could be implemented by SSR. It evidently reduce the computational cost for this method.