DISTRIBUTION FUNCTION OF A DRIFTING ELECTRON-GAS FOR GENERAL ENERGY-BAND STRUCTURES

The usual application of the Lei-Ting balance equation method for treating electron transport problems makes use of a Fermi distribution function for the electron motion relative to the center of mass. It is pointed out that this presumes the existence of a moving frame of reference that is dynamically equivalent to the rest frame of reference, and this is only true for electrons with a constant effective mass. The method is thus inapplicable to problems where electrons governed by a general energy-band dispersion E(k) are important (such as in miniband conduction). It is demonstrated that this difficulty can be overcome by introducing a distribution function for a drifting electron gas by maximizing the entropy subject to a prescribed average drift velocity. The distribution function reduces directly to the usual Fermi distribution for electron motion relative to the center of mass in the special case of E(k)=($) over bar h(2)\k\(2)/2m*. This maximum entropy treatment of a drifting electron gas provides a physically more direct as well as a more general basis for the application of the balance equation method.

电离层行进式扰动的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.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

中新世秦安黄土一古土壤序列记录的古气候变化周期

Analysis of periodic oscillations of climate is very important in understanding the behavior of the climate system. Milankovitch hypothesis, which holds that the glacial-interglacial climatic cycles during the Quaternary were primarily driven by variations in orbital parameters, has been supported by substantial geological evidence. Continuous long-term and high-resolution records are crucial to detect how variations of Earth's orbital parameters affected climate before the Quaternary when the boundary conditions were significantly different. Qinan loess formed in the Miocene is nearly continuous aeolian deposit in northern China. Previous study has established a constrained chronology, which provides a basis to examine long-term climatic variations. One of important issues to untangle the mechanisms behind major climate changes is the investigation of climate cycles recorded in Qinan loess. In this paper, two climatic proxies, magnetic susceptibility and redness, are analyzed for QA-I section to evaluate climate cycles using Maximum entropy spectral analysis and Blackman-Tuckey method. Main conclusions are presented as following: Results exhibit significant peaks at periods of 100 ka, 64 ka, 41 ka, 30 ka and 23 ka, but also 1000 ka, 600 ka and 400 ka. These peaks correspond to the dominant periods of the Earth's orbit parameters, which indicates that the formation of the aeolian sediment in northern China might be primarily driven by variations in orbital parameters. Fluctuations with different cycles respectively dominated in different periods. Major shifts in the dominant cycles occurred at 20.3, 19.0, 17.9, 15.2, 12.5 and 11.3 Myr ago. The transition that happened at 17.9 Myr ago was synchronous with the uplift of the Tibetan Plateau, while others at 15.2, 12.5 and 11.3 Myr ago were in good agreement with the timing of the development of Antarctic ice sheet. Therefore we inferred that these shifts might be related to changes in global ice volume and/or the Tibetan uplift. 3. The strong period of 100 ka is observed between 17.9 and 15.2, and 12.5 and 11.3 Myr ago. Ice sheet-climate models that have been used to explain the cause of the 100 ka period since the middle Pleistocene couldn't be responsible for driving the 100 ka climate cycle in the Miocene in Northern China because of the different boundary of climatic conditions between the Quaternary and Miocene. Further investigation is needed to understand how this cycle became dominant in Qinan loess records during these two time segments.

甘肃西峰晚第三纪红粘土地层中的蜗牛化石记录及其古环境意义

In this study, 260 mollusk fossil samples from a Red Clay sequence at Xifeng, Gansu province, in the northern China were analyzed quantitatively. 12 fossil species and four fossil zones have been identified. Three main ecological groups were determined based on ecological requirement of each mollusk taxon. According to fossil composition and succession of three ecological groups, the author discussed the origin and sedimentary environment of the red clay deposits, and the process of ecological environmental changes as well as the variations of the East Asia monsoons during 6.2-2.4 Ma in the Loess Plateau. A preliminary study on periodicity of paleoclimatic changes was also conducted by using spectral analysis method. The main results and conclusions are presented as follows:A continuous land mollusk fossil sequence of 6.2-2.4 Ma from Xifeng Red Clay Formation has been established, which provided a basic data for studying the environmental changes during late Miocene to Pliocene.The study of composition and preservation condition of mollusk fossils reveals a terrestrial in situ ecological population in the Red Clay Formation. All of identifiable mollusk species are composed of terrestrial taxa, which support the view that the Red Clay is an eolian origin, similar to the overlying Quaternary loess deposits.The mollusk record reveals the processes of ecological and environmental changes during 6.2-2.4 Ma in the Loess Plateau. Climatic changes experienced cold and dry from 6.2-5.4 Ma, warm and wet during 5.4-4.5 Ma, mild and moderate from 4.5-3-4 Ma, to rapid cooling and drying after 3.4 Ma. From '5.4- 2.4 Ma, climate was stepwise cooling. The cooling trend is in good agreement with a general1 0global cooling trend during this period, as documented by marine 5 0 records.4. Three remarked ecological shifts took place in mollusk assemblages from 6.2-2.4 Ma, focused on about 5.4, 4.5 and 3.4 Ma. The warming shift around 5.4 Ma was probably related to the rising of the global temperature. The cooling shifts around 4,5 and 3.4 Ma however might be closely linked to the uplift of Tibet Plateau and the development of Northern Hemisphere ice sheet.The succession in mollusk ecological groups also recorded the variability of the East Asian winter and summer monsoon. The winter monsoon dominated two periods from 6.2-5.4 Ma and from 3.4-2.4 Ma, while the summer monsoon was strong during 5.4-4.5 Ma. The variations in winter and summer monsoons were in phase during 4.5-3.4 Ma. Monsoon regimes changed with the duration about 1 Ma, which roughly corresponds to the cycle driven by tectonic activity on the time scales of ICP-IO7 years. In addition, mollusk fossils recorded the large amplitude and high frequency fluctuations overlapped on 105-107 years climate cycle.The maximum entropy spectral analysis and filter-band analysis of total mollusk individuals and three typical ecological groups suggest that the climate changes controlled mainly by solar insolation had periods about 70 ka and 40 ka on the time scales of 105 during late Miocene-Pliocene. Climatic periodicity intensified from 4.0 Ma, which reflected strengthened forcing by high latitude ice volume.

1064nm倍频波长分离膜的制备与性能研究

通过对主膜系添加匹配层并借助计算机对膜系进行优化，设计出结构规整、性能优良的1064ilm倍频波长分离膜。用电子束蒸发及光电极值监控技术在K9玻璃基底上沉积薄膜，将样品置于空气中在260℃温度下进行3h热退火处理。然后用Lambda 900分光光度计测量了样品的光谱性能；用表面热透镜(STL)技术测量了样品的弱吸收值；用调Q脉冲激光装置测试了样品的抗激光损伤阈值(LIDT)。实验结果发现，样品的实验光谱性能与理论光谱性能有很好的一致性。退火前后其光谱性能几乎没有发生温漂，说明薄膜的温度稳定性好；同时退火使

高真空室中的水蒸气对HfO2薄膜性能的影响

本文研究了在镀膜过程中真空室内水蒸气的含量对HfO2薄膜物理性能的影响。用电子束蒸发和光电极值监控的方法在BK7基底上制备HfO2薄膜。利用残余气体分析仪在线监测了真空室内的残余气氛组成。分别用Lambda 900光谱仪、X射线衍射方法、表面热透镜技术和1064nm的激光器测试了薄膜的光学性能、微结构、吸收和激光损伤阈值。实验发现，附加冷阱装置有助于我们有效控制镀膜过程中的水汽含量，且在水蒸气含量较少的真空室内镀制的薄膜具有较高的折射率，较小的晶粒尺寸，较低的弱吸收值和较高的损伤阈值。

两种后处理方法对HfO2薄膜性能的影响

利用电子束蒸发和光电极值监控技术制备了氧化铪薄膜，并分别用两种后处理方法（空气中退火和氧等离子体轰击）对样品进行了处理．然后，对样品的透过率、吸收和抗激光损伤阈值进行了测试分析．实验结果表明，两种后处理方法都能不同程度地降低了氧化铪薄膜的吸收损耗、提高了抗激光损伤阈值．实验结果还表明，氧等离子体轰击的后处理效果明显优于热退火，样品的吸收平均值在氧等离子体后处理前后分别为34．8ppm和9．0ppm，而基频（1 064nm）激光损伤阈值分别为10．0J／cm^2和21．4J／cm^2．

Episodic evolution of prolactin gene in primates

In the present study, we obtained exon 2-5 of prolactin (PRL) gene from four primate species by PCR and sequencing. Adding other genes available in GenBank, we calculate amino acid substitution rates for prolactin gene in primate. Comparison of nonsynonymous substitution rate to synonymous substitution rate ratios shows no evidence of positive selection for any lineage of primate prolactin gene. According to this and the facts that (i) no sites under positive selection are inferred by using maximum-likelihood method; (ii) among 32 amino acid replacement that occurred along the rapid evolutionary phase, only two are included in the 40 functionally important residues, indicating that amino acid replacement tends to occur in those functionally unimportant residues; (iii) partial of prolactin function is replaced by placental lactogen in primate at the rapid evolutionary phase of prolactin gene, we thus deem that it is relaxation of purifying selection to some extent rather than positive selection that enforces the rapid evolution of primate prolactin gene.

Phylogenetic analyses of Fringillidae (Aves : Passeriformes) using mitochondrial tRNA gene sequences and secondary structure

Fringillidae is a large and diverse family of Passeriformes. So far, however, Fringillidae relationships deduced from morphological features and by a number of molecular approaches have remained unproven. Recently, much attention has been attracted to mitochondrial tRNA genes, whose sequence and secondary structural characteristics have shown to be useful for Acrodont Lizards and deep-branch phylogenetic studies. In order to identify useful phylogenetic markers and test Fringillidae relationships, we have sequenced three major clusters of mitochondrial tRNA genes from 15 Fringillidae, taxa. A coincident tree, with coturnix as outgroup, was obtained through Maximum-likelihood method using combined dataset of 11 mitochondrial tRNA gene sequences. The result was similar to that through Neighbor-joining but different from Maximum-parsimony methods. Phylogenetic trees constructed with stem-region sequences of 11 genes had many different topologies and lower confidence than with total sequences. On the other hand, some secondary structural characteristics may provide phylogenetic information on relatively short internal branches at under-genus level. In summary, our data indicate that mitochondrial tRNA genes can achieve high confidence on Fringillidae phylogeny at subfamily level, and stem-region sequences may be suitable only at above-family level. Secondary structural characteristics may also be useful to resolve phylogenetic relationship between different genera of Fringillidae with good performance.

中国松鼠科 Sciuridae (Non-flying squirrels) 动物的分类和分布及系统发育研究

松鼠科隶属于哺乳纲、啮齿目，是一类世界广布性动物，全世界共计36 属 234 种。本文通过形态分析及分子生物学方法，研究了我国松鼠科动物的分类及 系统发育关系，内容包括以下两个方面： 一、中国松鼠科动物的分类 本研究查看了中科院昆明动物研究所、中科院动物研究所、云南大学、四川 省林业科学研究院、四川大学、陕西师范大学、陕西动物研究所、西北大学、上 海自然历史博物馆、广东昆虫所馆藏该科各类群标本。同时，结合相关文献，对 各类群的外形特征进行了详细的描述和系统分类整理，我国共计10 属34 种。测 取可量性状数据共计20800 个（完整头骨共计1300 号），而且还运用多变量及单 变量分析方法进行数据分析，对部分类群的种下分类结果进行分析：1.探讨了分 布于中国境内的明纹花鼠（Tamiops mcclellandii）3 亚种的有效性。2.研究了 隐纹花鼠（Tamiops swinhoei）在中国的分化并描记了分布于四川马尔康地区的 隐纹花鼠为一新亚种（Tamiops swinhoei markamensis）。3.描记了分布于云南 省东北部昭通地区的赤腹松鼠为一新亚种：昭通亚种（Callosciurus erythraeus zhaotongensis）。4. 提出了分布于云南省无量山区的泊氏长吻松鼠为一新亚种： 无量山亚种（Dremomys pernyi wuliangshanensis）。 二、中国松鼠科动物的系统发育研究 2004－2005 年相继到滇西北高黎贡山、香格里拉，滇西南临沧地区，滇东 南马关、麻栗坡、富宁，滇东北昭通地区进行标本采集，基本收集到我国南方松 鼠科属、种样品进行分子生物学研究。 1.中国松鼠科动物各属的系统发育关系 在分析线粒体细胞色数b 基因（1040bp）序列的基础上，运用贝叶斯法 （Bayesian phylogenetic reconstruction approaches）、最简约法（Maximum parsimony method）、邻接法（Neighbor Joining method）构建中国松鼠科动物 所有属间系统发育树。结果表明：所分析的10 属分别聚集为3 群，巨松鼠属Ratufa单独成为一群，丽松鼠属Callosciurus、长吻松鼠属Dremomys、花松鼠属 Tamiops 以及线松鼠属Menetes 成为一群，松鼠属Sciurus、岩松鼠属 Sciurotamias、花鼠属Tamias、旱獭属Marmota 以及黄鼠属Spermophilus 成 为一群。在此基础上我们进一步探讨了中国松鼠科动物的分化与环境变化的关 系。 2.长吻松鼠属的系统发育关系 在分析线粒体细胞色数b 基因（378bp）序列的基础上，我们重建了长吻松 鼠属内5 个类群间的系统发育关系，结合对5 类群外部形态特征的详细比较，我 们探讨并论证了红腿长吻松鼠、橙喉长吻松鼠作为种的地位的有效性，同时论证 了该属内泊氏长吻松鼠最早分化，其次为橙腹长吻松鼠，第三是红腿长吻松鼠， 红颊长吻松鼠和橙喉长吻松鼠是最后分化的类群。

Cloning and sequence analysis of the partial sequence of the rbcL from Bryopsis hypnoides

The partial sequence of the rbcL from Bryopsis hypnoides, including the sequences of the upstream, extron and partial intron, was amplified by PCR and their sequences were determined. With Spinacia oleracea as the outgroup, neighbor-joining method and maximum parsimony method were used respectively to build phylogenetic trees according to the rbcL exon sequence among 13 species that were the typical species of six phyla. Two kinds of trees showed clearly that there were two groups among those species, the green lineage and the non-green lineage. And the relationships of algae in the green lineage were similar in the two trees but those in the non-green lineage were not consistent. Analysis of codon preference indicated that the codon preference of the rbcL exon of Bryopsis hypnoides distinctly differed from that of the relevant sequence of photosynthetic bacteria.

层状反倾向结构岩质边坡变形机理及稳定性研究

Toppling is a major failure model in anti-dip layered rock slopes. Because of the limited by testing means and experimental apparatus, present research on the deformation mechanism and stability analysis are mainly focus on the 2-Dimensional deformation, and the research really based on 3-Dimension is still limited. Therefore, based on the present research station, the article rely on the important hydroelectric project of typical anti-dip layered rock slopes -- The left bank slope of Long-tan hydropower-station in Guang Xi, China, and focused on the influencing factors, deformation mechanism and stability analysis of anti-dip layered rock slopes, three problems as follows are researched in this paper. (1) Deformation influencing factor analysis on ant-dip layered rock slopes Three influencing factors are included: geological factor, engineering factor and environmental factor. It is concluded that the toppling deformation of anti-dip layered rock slopes are more sensitive to geological and engineering factors, but less sensitive to environmental factor. In addition, the sensitivity of various factors to the rock toppling deformation is also arranged sequentially as follows: construction, gravitation, rainfall (underground water) and rock structure intensity, etc. (2) 3D deformation study on the anti-dip layered toppling rock slopes Used 3D Distinct Element Method (3DEC) analyzed the 3D deformation characteristic of anti-dip layered rock slops. It can be seen that the toppling characteristics are obvious when the inter-angle between slope direction and layer striking direction is under 20o, when the inter-angle is over 20o and equal or less than 40o，the toppling deformation characteristics decrease sharply with increase of inter-angle, when the inter-angle is over 40o , the slope deformation is not controlled by joints but influenced by other failure mode. Therefore, in order to quantify the toppling characteristics, a differential value of displacement vector angle between layered rock slope and block rock slope is proposed as a key index to distinguish failure model for anti-dip layered rock slopes, and it was applied to study the toppling of the rock slopes at Guangxi Long-tan hydropower station, China. The results indicated that the index was effective and instructive for analyzing the anti-dip layered rock slopes. (3) Stability analysis methods Because of the imperfection of some present slope analysis methods, based on slope failure mode and those three influencing factors, “slope stability entropy” method is defined in this paper, which makes good use of the sensitivity of relational matrix to influencing factors on slope stability and the qualification characteristics for information entropy to the irregularity of slope deformation. By this method, not only the randomness of geologic body on the base of dynamic analysis of slope failure mode is fully concerned, but also it makes the analysis time-saving and simple. Finally, the research findings were used to the engineering example successfully, and rational conclusion has been obtained.

超参数控制下的波阻抗约束反演

The seismic survey is the most effective geophysical method during exploration and development of oil/gas. As a main means in processing and interpreting seismic data, impedance inversion takes up a special position in seismic survey. This is because the impedance parameter is a ligament which connects seismic data with well-logging and geological information, while it is also essential in predicting reservoir properties and sand-body. In fact, the result of traditional impedance inversion is not ideal. This is because the mathematical inverse problem of impedance is poor-pose so that the inverse result has instability and multi-result, so it is necessary to introduce regularization. Most simple regularizations are presented in existent literature, there is a premise that the image(or model) is globally smooth. In fact, as an actual geological model, it not only has made of smooth region but also be separated by the obvious edge, the edge is very important attribute of geological model. It's difficult to preserve these characteristics of the model and to avoid an edge too smooth to clear. Thereby, in this paper, we propose a impedance inverse method controlled by hyperparameters with edge-preserving regularization, the inverse convergence speed and result would be improved. In order to preserve the edge, the potential function of regularization should satisfy nine conditions such as basic assumptions edge preservation and convergence assumptions etc. Eventually, a model with clear background and edge-abnormity can be acquired. The several potential functions and the corresponding weight functions are presented in this paper. The potential functionφLφHL andφGM can meet the need of inverse precision by calculating the models. For the local constant planar and quadric models, we respectively present the neighborhood system of Markov random field corresponding to the regularization term. We linearity nonlinear regularization by using half-quadratic regularization, it not only preserve the edge, and but also simplify the inversion, and can use some linear methods. We introduced two regularization parameters (or hyperparameters) λ2 and δ in the regularization term. λ2 is used to balance the influence between the data term and the transcendental term; δ is a calibrating parameter used to adjust the gradient value at the discontinuous position(or formation interface). Meanwhile, in the inverse procedure, it is important to select the initial value of hyperparameters and to change hyperparameters, these will then have influence on convergence speed and inverse effect. In this paper, we roughly give the initial value of hyperparameters by using a trend- curve of φ-(λ2, δ) and by a method of calculating the upper limit value of hyperparameters. At one time, we change hyperparameters by using a certain coefficient or Maximum Likelihood method, this can be simultaneously fulfilled with the inverse procedure. Actually, we used the Fast Simulated Annealing algorithm in the inverse procedure. This method overcame restrictions from the local extremum without depending on the initial value, and got a global optimal result. Meanwhile, we expound in detail the convergence condition of FSA, the metropolis receiving probability form Metropolis-Hasting, the thermal procession based on the Gibbs sample and other methods integrated with FSA. These content can help us to understand and improve FSA. Through calculating in the theoretic model and applying it to the field data, it is proved that the impedance inverse method in this paper has the advantage of high precision practicability and obvious effect.

Burst event detection in wall turbulence by WVITA method

Wavelet Variable Interval Time Average (WVITA) is introduced as a method incorporating burst event detection in wall turbulence. Wavelet transform is performed to unfold the longitudinal fluctuating velocity time series measured in the near wall region of a turbulent boundary layer using hot-film anemometer. This unfolding is both in time and in space simultaneously. The splitted kinetic of the longitudinal fluctuating velocity time series among different scales is obtained by integrating the square of wavelet coefficient modulus over temporal space. The time scale that related to burst events in wall turbulence passing through the fixed probe is ascertained by maximum criterion of the kinetic energy evolution across scales. Wavelet transformed localized variance of the fluctuating velocity time series at the maximum kinetic scale is put forward instead of localized short time average variance in Variable Interval Time Average (VITA) scheme. The burst event detection result shows that WVITA scheme can avoid erroneous judgement and solve the grouping problem more effectively which is caused by VITA scheme itself and can not be avoided by adjusting the threshold level or changing the short time average interval.