Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
In this paper, the multisensor data fusion technique of a fault tolerant integrated navigation system is discussed. A neural approach for data fusion is proposed for multisensor integrated systems. The simulation res...In this paper, the multisensor data fusion technique of a fault tolerant integrated navigation system is discussed. A neural approach for data fusion is proposed for multisensor integrated systems. The simulation results show that this neural approach for data fusion is feasible.展开更多
Remote transmission of log data is an urgent problem for service companies. Remote transmission technology of log data here refers to both the transmission solution in combination with the CifNet multi-well data manag...Remote transmission of log data is an urgent problem for service companies. Remote transmission technology of log data here refers to both the transmission solution in combination with the CifNet multi-well data management system to automate the transmission, storage, management, and retrieval of log data to reduce turn-over time. It is an applied digital signature technology to implement breakpoint transmission and error recovery and ensure the effectiveness and reliability of log data transmission.展开更多
The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid interv...The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid intervals. In this study, the cosine-modulated binomial window function (CMBWF)-based staggered-grid finite-difference method is proposed. Two new parameters, the modulated time and modulated range are used in the new window function and by adjusting these two parameters we obtain different characteristics of the main and side lobes of the amplitude response. Numerical dispersion analysis and elastic wavefield forward modeling suggests that the CMBWF method is more precise and less computationally costly than the conventional Taylor series expansion staggered-grid finite-difference method.展开更多
The constitutive relationships of Al-Mg-Si alloy deformed at various strain rates,temperatures and strains were studied.The microstructure evolution was quantitatively characterized and analyzed,including recrystalliz...The constitutive relationships of Al-Mg-Si alloy deformed at various strain rates,temperatures and strains were studied.The microstructure evolution was quantitatively characterized and analyzed,including recrystallization fraction,grain sizes,local misorientation,geometrically necessary dislocation and stored strain energy during hot deformation and subsequent heat treatment.The results show that the dislocation density and energy storage are linear with ln Z during hot deformation and subsequent heat treatment,indicating continuous recrystallization occurring in both processes.With higher ln Z,the dislocation density declines more sharply during subsequent heat treatment.When ln Z is less than 28,dislocation density becomes more stable with less reduction during subsequent heat treatment after hot deformation.As these dislocations distribute along low angle grain boundaries,the subgrain has good stability during subsequent heat treatment.The main recrystallization mechanism during hot deformation is continuous dynamic recrystallization,accompanied by geometric dynamic recrystallization at higher ln Z.展开更多
Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical mod...Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical modeling, the rotated staggered-grid method (RSM) is a modification of the standard staggered-grid method (SSM). The variable-order method is based on the method of variable-length spatial operators and wavefield propagation, and it calculates the real dispersion error by adapting different finite-difference orders to different velocities. In this study, the variable-order rotated staggered-grid method (VRSM) is developed after applying the variable-order method to RSM to solve the numerical dispersion problem of RSM in low-velocity regions and reduce the computation cost. Moreover, based on theoretical dispersion and the real dispersion error of wave propagation calculated with the wave separation method, the application of the original method is extended from acoustic to shear waves, and the calculation is modified from theoretical to time-varying values. A layered model and an overthrust model are used to demonstrate the applicability of VRSM. We also evaluate the order distribution, wave propagation, and computation time. The results suggest that the VRSM order distribution is reasonable and VRSM produces high-precision results with a minimal computation cost.展开更多
AIM To identify chromosomal copy number aberrations(CNAs) in early-stage hepatocellular carcinoma(HCC) and analyze whether they are correlated with patient prognosis.METHODS One hundred and twenty patients with early-...AIM To identify chromosomal copy number aberrations(CNAs) in early-stage hepatocellular carcinoma(HCC) and analyze whether they are correlated with patient prognosis.METHODS One hundred and twenty patients with early-stage HCC were enrolled in our study, with the collection of formalin fixed, paraffin-embedded(FFPE) specimens and clinicopathological data. Tumor areas were marked by certified pathologists on a hematoxylin and eosinstained slide, and cancer and adjacent non-cancerous tissues underwent extraction of DNA, which was analyzed with the Affymetrix Onco Scan platform to assess CNAs and loss of heterozygosity(LOH). Ten individuals with nonmalignant disease were used as the control group. Another cohort consisting of 40 patients with stage Ⅰ/Ⅱ HCC were enrolled to analyze gene expression and to correlate findings with the Onco Scan data.RESULTS Copy number amplifications occurred at chromosomes 1 q21.1-q44 and 8 q12.3-24.3 and deletions were found at 4 q13.1-q35.2, 8 p 23.2-21.1, 16 q23.3-24.3, and 17 p13.3-12, while LOH commonly occurred at 1 p32.3, 3 p21.31, 8 p23.2-21.1, 16 q22.1-24.3, and 17 p 13.3-11 in early-stage HCC. Using Cox regression analysis, we also found that a higher percentage of genome change(≥ 60%) was an independent factor for worse prognosis in early-stage HCC(P = 0.031). Among the 875 genes in the Onco Scan Gene Chip, six were independent predictors of worse disease-free survival, of which three were amplified(MYC, ELAC2, and SYK) and three were deleted(GAK, MECOM, and WRN). Further, patients with HCC who exhibited ≥ 3 CNAs involving these six genes have worse outcomes compared to those who had < 3 CNAs(P < 0.001). Similarly, Asian patients with stage I HCC from The Cancer Genome Atlas harboring CNAs with these genes were also predicted to have poorer outcomes.CONCLUSION Patients with early-stage HCC and increased genome change or CNAs involving MYC, ELAC2, SYK, GAK, MECOM, or WRN are at risk for poorer outcome after resection.展开更多
In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
A novel adaptive fault-tolerant control scheme in the differential algebraic framework was proposed for attitude control of a heavy lift launch vehicle (HLLV). By using purely mathematical transformations, the decou...A novel adaptive fault-tolerant control scheme in the differential algebraic framework was proposed for attitude control of a heavy lift launch vehicle (HLLV). By using purely mathematical transformations, the decoupled input-output representations of HLLV were derived, rendering three decoupled second-order systems, i.e., pitch, yaw and roll channels. Based on a new type of numerical differentiator, a differential algebraic observer (DAO) was proposed for estimating the system states and the generalized disturbances, including various disturbances and additive fault torques. Driven by DAOs, three improved proportional-integral- differential (PID) controllers with disturbance compensation were designed for pitch, yaw and roll control. All signals in the closed-loop system were guaranteed to be ultimately uniformly bounded by utilization of Lyapunov's indirect method. The convincing numerical simulations indicate that the proposed control scheme is successful in achieving high performance in the presence of parametric perturbations, external disturbances, noisy corruptions, and actuator faults.展开更多
Mixed integer linear programming (MILP) approach for simultaneous gross error detection and data reconciliation has been proved as an efficient way to adjust process data with material, energy, and other balance con...Mixed integer linear programming (MILP) approach for simultaneous gross error detection and data reconciliation has been proved as an efficient way to adjust process data with material, energy, and other balance constrains. But the efficiency will decrease significantly when this method is applled in a large-scale problem because there are too many binary variables involved. In this article, an improved method is proposed in order to gen- erate gross error candidates with reliability factors before data rectification. Candidates are used in the MILP objec- tive function to improve the efficiency and accuracy by reducing the number of binary variables and giving accurate weights for suspected gross errors candidates. Performance of this improved method is compared and discussed by applying the algorithm in a widely used industrial example.展开更多
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
Filtration and cross-flow ultrafiltration techniques were used to separate culture media of Prorocentrurn donghaiense at the exponential growth, stationary and decline stages into 〈0.45 μm filtrate, 100 kDa-0.45 μm...Filtration and cross-flow ultrafiltration techniques were used to separate culture media of Prorocentrurn donghaiense at the exponential growth, stationary and decline stages into 〈0.45 μm filtrate, 100 kDa-0.45 μm, 1%100 kDa and 1-10 kDa retentate and 〈1 kDa ultrafiltrate fractions. The fluorescence properties of different molecular weights of dissolved organic matter (DOM) were measured by excitation-emission matrix spectra. Protein-like and humic-like fluorophores were observed in the DOM produced by P. donghaiense. The central positions of protein-like fluorophores showed a red shift with prolonged growth duration, shifting from tyrosine-like properties at the exponential growth stage to tryptophan-like properties at the stationary and decline stages. The excitation wavelengths of protein-like fluorophores exhibited some change in the exponential growth and stationary stages with increased molecular size, but showed little change in the decline stage. However, the emission wavelengths in the decline stage exhibited a blue shift. Very distinct C type and A type peaks in humic-like fluorophores were observed. With a prolonged culture time, the intensities of both of the peaks became strong and the excitation wavelengths of peak A showed a red shift, while the A:C ratios fell. More than 94% of fluorescent DOM was in the lower than 1 kDa molecular weight fraction.展开更多
Multi-source multi-class classification methods based on multi-class Support Vector Machines and data fusion strategies are proposed in this paper. The centralized and distributed fusion schemes are applied to combine...Multi-source multi-class classification methods based on multi-class Support Vector Machines and data fusion strategies are proposed in this paper. The centralized and distributed fusion schemes are applied to combine information from several data sources. In the centralized scheme, all information from several data sources is centralized to construct an input space. Then a multi-class Support Vector Machine classifier is trained. In the distributed schemes, the individual data sources are proc-essed separately and modelled by using the multi-class Support Vector Machine. Then new data fusion strategies are proposed to combine the information from the individual multi-class Support Vector Machine models. Our proposed fusion strategies take into account that an Support Vector Machine (SVM) classifier achieves classification by finding the optimal classification hyperplane with maximal margin. The proposed methods are applied for fault diagnosis of a diesel engine. The experimental results showed that almost all the proposed approaches can largely improve the diagnostic accuracy. The robustness of diagnosis is also improved because of the implementation of data fusion strategies. The proposed methods can also be applied in other fields.展开更多
In data post-processing for quantum key distribution, it is essential to have a highly efficient error reconciliation protocol. Based on the key redistribution scheme, we analyze a one-way error reconciliation protoco...In data post-processing for quantum key distribution, it is essential to have a highly efficient error reconciliation protocol. Based on the key redistribution scheme, we analyze a one-way error reconciliation protocol by data simulation. The relationship between the error correction capability and the key generation efficiency of three kinds of Hamming code are demonstrated. The simulation results indicate that when the initial error rates are (0,1.5%], (1.5,4%], and (4,11%], using the Hamming (31,26), (15,11), and (7,4) codes to correct the error, respectively, the key generation rate will be maximized. Based on this, we propose a modified one-way error reconciliation protocol which employs a mixed Hamming code concatenation scheme. The error correction capability and key generation rate are verified through data simulation. Using the parameters of the posterior distribution based on the tested data, a simple method for estimating the bit error rate (BER) with a given confidence interval is estimated. The simulation results show that when the initial bit error rate is 10.00%, after 7 rounds of error correction, the error bits are eliminated completely, and the key generation rate is 10.36%; the BER expectation is 2.96×10^-10, and when the confidence is 95% the corresponding BER upper limit is 2.17×10^-9. By comparison, for the single (7,4) Hamming code error reconciliation scheme at a confidence of 95%,the key generation rate is only 6.09%, while the BER expectation is 5.92x 10"9, with a BER upper limit of 4.34×10^-8. Hence, our improved protocol is much better than the original one.展开更多
Background error covariance plays an important role in any variational data assimilation system, because it determines how information from observations is spread in model space and between different model variables. ...Background error covariance plays an important role in any variational data assimilation system, because it determines how information from observations is spread in model space and between different model variables. In this paper, the use of orthogonal wavelets in representation of background error covariance over a limited area is studied. Based on the WRF model and its 3D-VAR system, an algorithm using orthogonal wavelets to model background error covariance is developed. Because each wavelet function contains information on both position and scale, using a diagonal correlation matrix in wavelet space gives the possibility to represent some anisotropic and inhomogeneous characteristics of background error covariance. The experiments show that local correlation functions are better modeled than spectral methods. The formulation of wavelet background error covariance is tested with the typhoon Kaemi (2006). The results of experiments indicate that the subsequent forecasts of typhoon Kaemi’s track and intensity are significantly improved by the new method.展开更多
The purpose of this article is to provide the inversion relationships between the reciprocal sum S(1, 2,…, m) and the alternating sum T(1, 2,…, m) for generalized Lucas numbers which generalizes the Melham's re...The purpose of this article is to provide the inversion relationships between the reciprocal sum S(1, 2,…, m) and the alternating sum T(1, 2,…, m) for generalized Lucas numbers which generalizes the Melham's results.展开更多
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
文摘In this paper, the multisensor data fusion technique of a fault tolerant integrated navigation system is discussed. A neural approach for data fusion is proposed for multisensor integrated systems. The simulation results show that this neural approach for data fusion is feasible.
文摘Remote transmission of log data is an urgent problem for service companies. Remote transmission technology of log data here refers to both the transmission solution in combination with the CifNet multi-well data management system to automate the transmission, storage, management, and retrieval of log data to reduce turn-over time. It is an applied digital signature technology to implement breakpoint transmission and error recovery and ensure the effectiveness and reliability of log data transmission.
基金supported by the National Major Research Equipment Development Projects(No.ZDYZ2012-1-02-04)the National Natural Science Foundation of China(No.41474106)
文摘The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid intervals. In this study, the cosine-modulated binomial window function (CMBWF)-based staggered-grid finite-difference method is proposed. Two new parameters, the modulated time and modulated range are used in the new window function and by adjusting these two parameters we obtain different characteristics of the main and side lobes of the amplitude response. Numerical dispersion analysis and elastic wavefield forward modeling suggests that the CMBWF method is more precise and less computationally costly than the conventional Taylor series expansion staggered-grid finite-difference method.
基金Project(2016YFB0300901)supported by the National Key R&D Program of ChinaProject(TC190H3ZV/2)supported by the National Building Project of Application Demonstration Platform on New Materials Products,ChinaProject(15308469779)supported by Key Laboratory of National Science and Technology for Equipment Pre-research,China。
文摘The constitutive relationships of Al-Mg-Si alloy deformed at various strain rates,temperatures and strains were studied.The microstructure evolution was quantitatively characterized and analyzed,including recrystallization fraction,grain sizes,local misorientation,geometrically necessary dislocation and stored strain energy during hot deformation and subsequent heat treatment.The results show that the dislocation density and energy storage are linear with ln Z during hot deformation and subsequent heat treatment,indicating continuous recrystallization occurring in both processes.With higher ln Z,the dislocation density declines more sharply during subsequent heat treatment.When ln Z is less than 28,dislocation density becomes more stable with less reduction during subsequent heat treatment after hot deformation.As these dislocations distribute along low angle grain boundaries,the subgrain has good stability during subsequent heat treatment.The main recrystallization mechanism during hot deformation is continuous dynamic recrystallization,accompanied by geometric dynamic recrystallization at higher ln Z.
基金supported by the National Science and Technology Major Project of China(No.2011ZX05004-003)the National Basic Research Program of China(No.2013CB228602)the National High Tech Research Program of China(No.2013AA064202)
文摘Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical modeling, the rotated staggered-grid method (RSM) is a modification of the standard staggered-grid method (SSM). The variable-order method is based on the method of variable-length spatial operators and wavefield propagation, and it calculates the real dispersion error by adapting different finite-difference orders to different velocities. In this study, the variable-order rotated staggered-grid method (VRSM) is developed after applying the variable-order method to RSM to solve the numerical dispersion problem of RSM in low-velocity regions and reduce the computation cost. Moreover, based on theoretical dispersion and the real dispersion error of wave propagation calculated with the wave separation method, the application of the original method is extended from acoustic to shear waves, and the calculation is modified from theoretical to time-varying values. A layered model and an overthrust model are used to demonstrate the applicability of VRSM. We also evaluate the order distribution, wave propagation, and computation time. The results suggest that the VRSM order distribution is reasonable and VRSM produces high-precision results with a minimal computation cost.
基金Supported by the Chang Gung Memorial Hospital in Taiwan,No.CMRPG 3C0951-3 and No.CMRPG 3A0671 to Yu MC,and No.CMRPD3F0011 to Tsai CN
文摘AIM To identify chromosomal copy number aberrations(CNAs) in early-stage hepatocellular carcinoma(HCC) and analyze whether they are correlated with patient prognosis.METHODS One hundred and twenty patients with early-stage HCC were enrolled in our study, with the collection of formalin fixed, paraffin-embedded(FFPE) specimens and clinicopathological data. Tumor areas were marked by certified pathologists on a hematoxylin and eosinstained slide, and cancer and adjacent non-cancerous tissues underwent extraction of DNA, which was analyzed with the Affymetrix Onco Scan platform to assess CNAs and loss of heterozygosity(LOH). Ten individuals with nonmalignant disease were used as the control group. Another cohort consisting of 40 patients with stage Ⅰ/Ⅱ HCC were enrolled to analyze gene expression and to correlate findings with the Onco Scan data.RESULTS Copy number amplifications occurred at chromosomes 1 q21.1-q44 and 8 q12.3-24.3 and deletions were found at 4 q13.1-q35.2, 8 p 23.2-21.1, 16 q23.3-24.3, and 17 p13.3-12, while LOH commonly occurred at 1 p32.3, 3 p21.31, 8 p23.2-21.1, 16 q22.1-24.3, and 17 p 13.3-11 in early-stage HCC. Using Cox regression analysis, we also found that a higher percentage of genome change(≥ 60%) was an independent factor for worse prognosis in early-stage HCC(P = 0.031). Among the 875 genes in the Onco Scan Gene Chip, six were independent predictors of worse disease-free survival, of which three were amplified(MYC, ELAC2, and SYK) and three were deleted(GAK, MECOM, and WRN). Further, patients with HCC who exhibited ≥ 3 CNAs involving these six genes have worse outcomes compared to those who had < 3 CNAs(P < 0.001). Similarly, Asian patients with stage I HCC from The Cancer Genome Atlas harboring CNAs with these genes were also predicted to have poorer outcomes.CONCLUSION Patients with early-stage HCC and increased genome change or CNAs involving MYC, ELAC2, SYK, GAK, MECOM, or WRN are at risk for poorer outcome after resection.
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金Foundation item: Project(2012M521538) supported by China Postdoctoral Science Foundation Project suppolted by Postdoctoral Science Foundation of Central South University
文摘A novel adaptive fault-tolerant control scheme in the differential algebraic framework was proposed for attitude control of a heavy lift launch vehicle (HLLV). By using purely mathematical transformations, the decoupled input-output representations of HLLV were derived, rendering three decoupled second-order systems, i.e., pitch, yaw and roll channels. Based on a new type of numerical differentiator, a differential algebraic observer (DAO) was proposed for estimating the system states and the generalized disturbances, including various disturbances and additive fault torques. Driven by DAOs, three improved proportional-integral- differential (PID) controllers with disturbance compensation were designed for pitch, yaw and roll control. All signals in the closed-loop system were guaranteed to be ultimately uniformly bounded by utilization of Lyapunov's indirect method. The convincing numerical simulations indicate that the proposed control scheme is successful in achieving high performance in the presence of parametric perturbations, external disturbances, noisy corruptions, and actuator faults.
基金Supported by the National High Technology Research and Development Program of China (2007AA40702 and 2007AA04Z191)
文摘Mixed integer linear programming (MILP) approach for simultaneous gross error detection and data reconciliation has been proved as an efficient way to adjust process data with material, energy, and other balance constrains. But the efficiency will decrease significantly when this method is applled in a large-scale problem because there are too many binary variables involved. In this article, an improved method is proposed in order to gen- erate gross error candidates with reliability factors before data rectification. Candidates are used in the MILP objec- tive function to improve the efficiency and accuracy by reducing the number of binary variables and giving accurate weights for suspected gross errors candidates. Performance of this improved method is compared and discussed by applying the algorithm in a widely used industrial example.
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
基金Supported by the High Technology Research and Development Program of China (863 Program) (Nos 2006AA09Z180 and 2004AA639790)the National Natural Science Foundation of China (No 40106013)the National Basic Research Program of China (973 Program) (No 2001CB409703)
文摘Filtration and cross-flow ultrafiltration techniques were used to separate culture media of Prorocentrurn donghaiense at the exponential growth, stationary and decline stages into 〈0.45 μm filtrate, 100 kDa-0.45 μm, 1%100 kDa and 1-10 kDa retentate and 〈1 kDa ultrafiltrate fractions. The fluorescence properties of different molecular weights of dissolved organic matter (DOM) were measured by excitation-emission matrix spectra. Protein-like and humic-like fluorophores were observed in the DOM produced by P. donghaiense. The central positions of protein-like fluorophores showed a red shift with prolonged growth duration, shifting from tyrosine-like properties at the exponential growth stage to tryptophan-like properties at the stationary and decline stages. The excitation wavelengths of protein-like fluorophores exhibited some change in the exponential growth and stationary stages with increased molecular size, but showed little change in the decline stage. However, the emission wavelengths in the decline stage exhibited a blue shift. Very distinct C type and A type peaks in humic-like fluorophores were observed. With a prolonged culture time, the intensities of both of the peaks became strong and the excitation wavelengths of peak A showed a red shift, while the A:C ratios fell. More than 94% of fluorescent DOM was in the lower than 1 kDa molecular weight fraction.
文摘Multi-source multi-class classification methods based on multi-class Support Vector Machines and data fusion strategies are proposed in this paper. The centralized and distributed fusion schemes are applied to combine information from several data sources. In the centralized scheme, all information from several data sources is centralized to construct an input space. Then a multi-class Support Vector Machine classifier is trained. In the distributed schemes, the individual data sources are proc-essed separately and modelled by using the multi-class Support Vector Machine. Then new data fusion strategies are proposed to combine the information from the individual multi-class Support Vector Machine models. Our proposed fusion strategies take into account that an Support Vector Machine (SVM) classifier achieves classification by finding the optimal classification hyperplane with maximal margin. The proposed methods are applied for fault diagnosis of a diesel engine. The experimental results showed that almost all the proposed approaches can largely improve the diagnostic accuracy. The robustness of diagnosis is also improved because of the implementation of data fusion strategies. The proposed methods can also be applied in other fields.
基金supported in part by the Foundation for Key Program of Chinese Ministry of Education under Grant No.212177Scientific Research Foundation of the Education Department of Shaanxi Province under Grant No.12JK0973
文摘In data post-processing for quantum key distribution, it is essential to have a highly efficient error reconciliation protocol. Based on the key redistribution scheme, we analyze a one-way error reconciliation protocol by data simulation. The relationship between the error correction capability and the key generation efficiency of three kinds of Hamming code are demonstrated. The simulation results indicate that when the initial error rates are (0,1.5%], (1.5,4%], and (4,11%], using the Hamming (31,26), (15,11), and (7,4) codes to correct the error, respectively, the key generation rate will be maximized. Based on this, we propose a modified one-way error reconciliation protocol which employs a mixed Hamming code concatenation scheme. The error correction capability and key generation rate are verified through data simulation. Using the parameters of the posterior distribution based on the tested data, a simple method for estimating the bit error rate (BER) with a given confidence interval is estimated. The simulation results show that when the initial bit error rate is 10.00%, after 7 rounds of error correction, the error bits are eliminated completely, and the key generation rate is 10.36%; the BER expectation is 2.96×10^-10, and when the confidence is 95% the corresponding BER upper limit is 2.17×10^-9. By comparison, for the single (7,4) Hamming code error reconciliation scheme at a confidence of 95%,the key generation rate is only 6.09%, while the BER expectation is 5.92x 10"9, with a BER upper limit of 4.34×10^-8. Hence, our improved protocol is much better than the original one.
基金National Natural Science Foundation of China (40775064)
文摘Background error covariance plays an important role in any variational data assimilation system, because it determines how information from observations is spread in model space and between different model variables. In this paper, the use of orthogonal wavelets in representation of background error covariance over a limited area is studied. Based on the WRF model and its 3D-VAR system, an algorithm using orthogonal wavelets to model background error covariance is developed. Because each wavelet function contains information on both position and scale, using a diagonal correlation matrix in wavelet space gives the possibility to represent some anisotropic and inhomogeneous characteristics of background error covariance. The experiments show that local correlation functions are better modeled than spectral methods. The formulation of wavelet background error covariance is tested with the typhoon Kaemi (2006). The results of experiments indicate that the subsequent forecasts of typhoon Kaemi’s track and intensity are significantly improved by the new method.
基金Supported by the Natural Science Foundation of the Education Department of Henan Province(2003110009)
文摘The purpose of this article is to provide the inversion relationships between the reciprocal sum S(1, 2,…, m) and the alternating sum T(1, 2,…, m) for generalized Lucas numbers which generalizes the Melham's results.