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.展开更多
We derived an equation for saturation in carbonate reservoirs based on the electrical efficiency model in the case of lacking core data. Owing to the complex pore structure and strong heterogeneity in carbonate reserv...We derived an equation for saturation in carbonate reservoirs based on the electrical efficiency model in the case of lacking core data. Owing to the complex pore structure and strong heterogeneity in carbonate reservoirs, the relation between electrical efficiency and water porosity is either complex or linear. We proposed an electrical efficiency equation that accounts for the relation between electrical efficiency and water porosity. We also proposed a power-law relation between electrical efficiency and deep-formation resistivity and analyzed the factors controlling the error in the water saturation computations. We concluded that the calculation accuracy of the electrical efficiency is critical to the application of the saturation equation. The saturation equation was applied to the carbonate reservoirs of three wells in Iraq and Indonesia. For relative rock electrical efficiency error below 0.1, the water saturation absolute error is also below 0.1. Therefore, we infer that the proposed saturation equation generally satisfies the evaluation criteria for carbonate reservoirs.展开更多
A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan ...A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan variance method and each error coefficient is identified. Furthermore, a random drift error model for IFOG is built by the method of time series analysis. The conclusion provides supports for improving IFOG design and compensating for errors of IFOG in practice.展开更多
To analyze the errors of processing data, the testing principle for jet elements is introduced and the property of testing system is theoretically and experimentally studied. On the basis of the above, the method of p...To analyze the errors of processing data, the testing principle for jet elements is introduced and the property of testing system is theoretically and experimentally studied. On the basis of the above, the method of processing data is presented and the error formulae, which are the functions of the testing system property, are derived. Finally, the methods of reducing the errors are provided. The measured results are in correspondence with the theoretical conclusion.展开更多
In this paper, we define a new class of biased linear estimators of the vector of unknown parameters in the deficient_rank linear model based on the spectral decomposition expression of the best linear minimun bias es...In this paper, we define a new class of biased linear estimators of the vector of unknown parameters in the deficient_rank linear model based on the spectral decomposition expression of the best linear minimun bias estimator. Some important properties are discussed. By appropriate choices of bias parameters, we construct many interested and useful biased linear estimators, which are the extension of ordinary biased linear estimators in the full_rank linear model to the deficient_rank linear model. At last, we give a numerical example in geodetic adjustment.展开更多
A condition number is an amplification coefficient due to errors in computing. Thus the theory of condition numbers plays an important role in error analysis. In this paper, following the approach of Rice, condition n...A condition number is an amplification coefficient due to errors in computing. Thus the theory of condition numbers plays an important role in error analysis. In this paper, following the approach of Rice, condition numbers are defined for factors of some matrix factorizations such as the Cholesky factorization of a symmetric positive definite matrix and QR factorization of a general matrix. The condition numbers are derived by a technique of analytic expansion of the factor dependent on one parameter and matrix-vector equation. Condition numbers of the Cholesky and QR factors are different from the ones previously introduced by other authors, but similar to Chang's results. In Cholesky factorization, corresponding with the condition number of the factor matrix L , K _L is a low bound of Stewart's condition number K .展开更多
Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network mode...Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.展开更多
The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was de...The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was devised. The technique regarded the least square sphere center as the initial center of the concentric spheres containing all measurement points, and then the center was moved gradually to reduce the radial separation till the minimum radial separation center was got where the constructed concentric spheres conformed to the minimum zone condition. The method was modeled firstly, then the geometric approximation process was analyzed, and finally,the software for data processing was programmed. As evaluation example, five steel balls were measured and the measurement data were processed with the developed program. The average iteration times of the approximation technique is 4.2, and on average the obtained sphericity error is 0. 529μm smaller than the least square solution,with accuracy increased by 7. 696%.展开更多
The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
Using lAP AGCM simulation results for the period 1961-2005, summer hot days in China were calculated and then compared with observations. Generally, the spatial pattern of hot days is reasonably reproduced, with more ...Using lAP AGCM simulation results for the period 1961-2005, summer hot days in China were calculated and then compared with observations. Generally, the spatial pattern of hot days is reasonably reproduced, with more hot days found in northern China, the Yangtze and Huaihe River basin, the Chuan-Yu region, and southern Xinjiang. However, the model tends to overestimate the number of hot days in the above-mentioned regions, particularly in the Yangtze and Huaihe River basin where the simulated summer-mean hot days is 13 days more than observed when averaged over the whole region, and the maximum overestimation of hot days can reach 23 days in the region. Analysis of the probability distribution of daily maximum temperature (Trnax) suggests that the warm bias in the model-simulated Tmax contributes largely to the overestimation of hot days in the model. Furthermore, the discrepancy in the simulated variance of the Tmax distribution also plays a non- negligible role in the overestimation of hot days. Indeed, the latter can even account for 22% of the total bias of simulated hot days in August in the Yangtze and Huaihe River basin. The quantification of model bias from the mean value and variability can provide more information for further model improvement.展开更多
Array calibration is important in engineering practice. In this paper, fast calibration methods for a ULA's gain and phase errors both in far and near fields are proposed. In the far field, using a single sound so...Array calibration is important in engineering practice. In this paper, fast calibration methods for a ULA's gain and phase errors both in far and near fields are proposed. In the far field, using a single sound source without exact orientation, this method horizontally rotates the array exactly once, performs eigen value decomposition for the covariance matrix of received data, then computes the gain and phase error according to the formulas. In the near field, using the same single sound source, it is necessary to rotate the array horizontally at most three times, build equations according to geometric relations, then solve them. Using the formula proposed in this paper, spherical waves are modified into plane waves. Then eigen values decomposition is performed. These two calibration methods were shown to be valid by simulation and are fast, accurate and easy to use. Finally, an analysis of factors influencing estimation precision is given.展开更多
A new multi-sensor data fusion algorithm based on EMD-MMSE was proposed.Empirical mode decomposition(EMD)is used to extract the noise of every time series for estimating the variance of the noise.Then minimum mean squ...A new multi-sensor data fusion algorithm based on EMD-MMSE was proposed.Empirical mode decomposition(EMD)is used to extract the noise of every time series for estimating the variance of the noise.Then minimum mean square error(MMSE)estimator is used to calculate the weights of the corresponding series.Finally,the fused signal is the weighted addition of all these series.The experiments in lab testified the efficiency of this method.In addition,the comparison in fusion time and fusion results with existing fusion method based on wavelet and average technique shows the advantage of this method greatly.展开更多
This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound cor...This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.展开更多
Branched continued fractions are one of the multidimensional generalization of the continued fractions. Branched continued fractions with not equivalent variables are an analog of the regular C-fractions for multiple ...Branched continued fractions are one of the multidimensional generalization of the continued fractions. Branched continued fractions with not equivalent variables are an analog of the regular C-fractions for multiple power series. We consider 1-periodic branched continued fraction of the special form which is an analog fraction with not equivalent variables if the values of that variables are fixed. We establish an analog of the parabola theorem for that fraction and estimate truncation error bounds for that fractions at some restrictions. We also propose to use weight coefficients for obtaining different parabolic regions for the same fraction without any additional restriction for first element.展开更多
基金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 Science and Technology Major Project(2011ZX05030)
文摘We derived an equation for saturation in carbonate reservoirs based on the electrical efficiency model in the case of lacking core data. Owing to the complex pore structure and strong heterogeneity in carbonate reservoirs, the relation between electrical efficiency and water porosity is either complex or linear. We proposed an electrical efficiency equation that accounts for the relation between electrical efficiency and water porosity. We also proposed a power-law relation between electrical efficiency and deep-formation resistivity and analyzed the factors controlling the error in the water saturation computations. We concluded that the calculation accuracy of the electrical efficiency is critical to the application of the saturation equation. The saturation equation was applied to the carbonate reservoirs of three wells in Iraq and Indonesia. For relative rock electrical efficiency error below 0.1, the water saturation absolute error is also below 0.1. Therefore, we infer that the proposed saturation equation generally satisfies the evaluation criteria for carbonate reservoirs.
文摘A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan variance method and each error coefficient is identified. Furthermore, a random drift error model for IFOG is built by the method of time series analysis. The conclusion provides supports for improving IFOG design and compensating for errors of IFOG in practice.
文摘To analyze the errors of processing data, the testing principle for jet elements is introduced and the property of testing system is theoretically and experimentally studied. On the basis of the above, the method of processing data is presented and the error formulae, which are the functions of the testing system property, are derived. Finally, the methods of reducing the errors are provided. The measured results are in correspondence with the theoretical conclusion.
文摘In this paper, we define a new class of biased linear estimators of the vector of unknown parameters in the deficient_rank linear model based on the spectral decomposition expression of the best linear minimun bias estimator. Some important properties are discussed. By appropriate choices of bias parameters, we construct many interested and useful biased linear estimators, which are the extension of ordinary biased linear estimators in the full_rank linear model to the deficient_rank linear model. At last, we give a numerical example in geodetic adjustment.
文摘A condition number is an amplification coefficient due to errors in computing. Thus the theory of condition numbers plays an important role in error analysis. In this paper, following the approach of Rice, condition numbers are defined for factors of some matrix factorizations such as the Cholesky factorization of a symmetric positive definite matrix and QR factorization of a general matrix. The condition numbers are derived by a technique of analytic expansion of the factor dependent on one parameter and matrix-vector equation. Condition numbers of the Cholesky and QR factors are different from the ones previously introduced by other authors, but similar to Chang's results. In Cholesky factorization, corresponding with the condition number of the factor matrix L , K _L is a low bound of Stewart's condition number K .
基金Project(51321065)supported by the Innovative Research Groups of the National Natural Science Foundation of ChinaProject(2013CB035904)supported by the National Basic Research Program of China(973 Program)Project(51439005)supported by the National Natural Science Foundation of China
文摘Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.
基金Supported by National Natural Science Foundation of China(No.50175081) andTianjin Municipal Science and Technology Commission (No.0431835116).
文摘The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was devised. The technique regarded the least square sphere center as the initial center of the concentric spheres containing all measurement points, and then the center was moved gradually to reduce the radial separation till the minimum radial separation center was got where the constructed concentric spheres conformed to the minimum zone condition. The method was modeled firstly, then the geometric approximation process was analyzed, and finally,the software for data processing was programmed. As evaluation example, five steel balls were measured and the measurement data were processed with the developed program. The average iteration times of the approximation technique is 4.2, and on average the obtained sphericity error is 0. 529μm smaller than the least square solution,with accuracy increased by 7. 696%.
文摘The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
基金supported by the Special Scientific Research Fund of the Meteorological Public Welfare Profession of China[grant number GYHY01406021]National Key Research and Development Program[grant number 2016YFC0402702]the National Natural Science Foundation of China[grant numbers 41575095,41175073]
文摘Using lAP AGCM simulation results for the period 1961-2005, summer hot days in China were calculated and then compared with observations. Generally, the spatial pattern of hot days is reasonably reproduced, with more hot days found in northern China, the Yangtze and Huaihe River basin, the Chuan-Yu region, and southern Xinjiang. However, the model tends to overestimate the number of hot days in the above-mentioned regions, particularly in the Yangtze and Huaihe River basin where the simulated summer-mean hot days is 13 days more than observed when averaged over the whole region, and the maximum overestimation of hot days can reach 23 days in the region. Analysis of the probability distribution of daily maximum temperature (Trnax) suggests that the warm bias in the model-simulated Tmax contributes largely to the overestimation of hot days in the model. Furthermore, the discrepancy in the simulated variance of the Tmax distribution also plays a non- negligible role in the overestimation of hot days. Indeed, the latter can even account for 22% of the total bias of simulated hot days in August in the Yangtze and Huaihe River basin. The quantification of model bias from the mean value and variability can provide more information for further model improvement.
文摘Array calibration is important in engineering practice. In this paper, fast calibration methods for a ULA's gain and phase errors both in far and near fields are proposed. In the far field, using a single sound source without exact orientation, this method horizontally rotates the array exactly once, performs eigen value decomposition for the covariance matrix of received data, then computes the gain and phase error according to the formulas. In the near field, using the same single sound source, it is necessary to rotate the array horizontally at most three times, build equations according to geometric relations, then solve them. Using the formula proposed in this paper, spherical waves are modified into plane waves. Then eigen values decomposition is performed. These two calibration methods were shown to be valid by simulation and are fast, accurate and easy to use. Finally, an analysis of factors influencing estimation precision is given.
基金The National High Technology Research and Development Program of China(863Program)(No.2001AA602021)
文摘A new multi-sensor data fusion algorithm based on EMD-MMSE was proposed.Empirical mode decomposition(EMD)is used to extract the noise of every time series for estimating the variance of the noise.Then minimum mean square error(MMSE)estimator is used to calculate the weights of the corresponding series.Finally,the fused signal is the weighted addition of all these series.The experiments in lab testified the efficiency of this method.In addition,the comparison in fusion time and fusion results with existing fusion method based on wavelet and average technique shows the advantage of this method greatly.
文摘This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.
文摘Branched continued fractions are one of the multidimensional generalization of the continued fractions. Branched continued fractions with not equivalent variables are an analog of the regular C-fractions for multiple power series. We consider 1-periodic branched continued fraction of the special form which is an analog fraction with not equivalent variables if the values of that variables are fixed. We establish an analog of the parabola theorem for that fraction and estimate truncation error bounds for that fractions at some restrictions. We also propose to use weight coefficients for obtaining different parabolic regions for the same fraction without any additional restriction for first element.