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.展开更多
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 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, 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 .展开更多
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.
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.展开更多
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.展开更多
Butterfly spring-relief valve, a crucial safety attachment of pressure vessel, is used to prevent pressuresystem from exceeding allowable limit value. Safe, expeditious and accurate calibration of safety valves is con...Butterfly spring-relief valve, a crucial safety attachment of pressure vessel, is used to prevent pressuresystem from exceeding allowable limit value. Safe, expeditious and accurate calibration of safety valves is consequentlyof vital importance to safe and economic operation of generating units. NSH CALIBRATOR could complete, not only theon-line performance and parameter tests of safety valves within two to five seconds with opening pressure of safetyvalves and steam medium pressure automatically recorded, but also could complete the on-line adjustment of safetyvalves verified unqualified. It saves energy consumption, decreases noise pollution and improves accuracy and efficiencyof safety valve calibration.[展开更多
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.展开更多
To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat s...To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.展开更多
Based on the property that high degree variable nodes within an irregular LowDensity Parity-Check (LDPC) code have more powerful error-correcting capability than that of low degree variable nodes, a group of irregular...Based on the property that high degree variable nodes within an irregular LowDensity Parity-Check (LDPC) code have more powerful error-correcting capability than that of low degree variable nodes, a group of irregular LDPC codes with Unequal Error Protection (UEP)property is designed in this letter. Simulation results show that the transmission quality of the image may be effectively improved with this class of irregular LDPC code.展开更多
The theoretical model of axial ultrasonic vibration grinding force is built on the basis of a mathematical model of cutting deforming force deduced from the assumptions of thickness of the undeformed debris under Rayl...The theoretical model of axial ultrasonic vibration grinding force is built on the basis of a mathematical model of cutting deforming force deduced from the assumptions of thickness of the undeformed debris under Rayleigh distribution and a mathematical model of friction based on the theoretical analysis of relative sliding velocity of abrasive and workpiece. Then, the coefficients of the ultrasonic vibration grinding force model are calculated through analysis of nonlinear regression of the theoretical model by using MATLAB, and the law of influence of grinding depth, workpiece speed, frequency and amplitude of the mill on the grinding force is summarized after applying the model to analyze the ultrasonic grinding force. The result of the above-mentioned law shows that the grinding force decreases as frequency and amplitude increase, while increases as grinding depth and workpiece speed increase; the maximum relative error of prediction and experimental values of the normal grinding force is 11.47% and its average relative error is 5.41%; the maximum relative error of the tangential grinding force is 10.14% and its average relative error is 4.29%. The result of employing regression equation to predict ultrasonic grinding force approximates to the experimental data, therefore the accuracy and reliability of the model is verified.展开更多
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.展开更多
基金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.
文摘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.
基金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.
文摘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 .
基金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.
基金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.
文摘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.
文摘Butterfly spring-relief valve, a crucial safety attachment of pressure vessel, is used to prevent pressuresystem from exceeding allowable limit value. Safe, expeditious and accurate calibration of safety valves is consequentlyof vital importance to safe and economic operation of generating units. NSH CALIBRATOR could complete, not only theon-line performance and parameter tests of safety valves within two to five seconds with opening pressure of safetyvalves and steam medium pressure automatically recorded, but also could complete the on-line adjustment of safetyvalves verified unqualified. It saves energy consumption, decreases noise pollution and improves accuracy and efficiencyof safety valve calibration.[
文摘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.
基金Supported by China Postdoctoral Science Foundation under Grant No.20090460102
文摘To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.
文摘Based on the property that high degree variable nodes within an irregular LowDensity Parity-Check (LDPC) code have more powerful error-correcting capability than that of low degree variable nodes, a group of irregular LDPC codes with Unequal Error Protection (UEP)property is designed in this letter. Simulation results show that the transmission quality of the image may be effectively improved with this class of irregular LDPC code.
基金Project(51275530)supported by the National Natural Science Foundation of China
文摘The theoretical model of axial ultrasonic vibration grinding force is built on the basis of a mathematical model of cutting deforming force deduced from the assumptions of thickness of the undeformed debris under Rayleigh distribution and a mathematical model of friction based on the theoretical analysis of relative sliding velocity of abrasive and workpiece. Then, the coefficients of the ultrasonic vibration grinding force model are calculated through analysis of nonlinear regression of the theoretical model by using MATLAB, and the law of influence of grinding depth, workpiece speed, frequency and amplitude of the mill on the grinding force is summarized after applying the model to analyze the ultrasonic grinding force. The result of the above-mentioned law shows that the grinding force decreases as frequency and amplitude increase, while increases as grinding depth and workpiece speed increase; the maximum relative error of prediction and experimental values of the normal grinding force is 11.47% and its average relative error is 5.41%; the maximum relative error of the tangential grinding force is 10.14% and its average relative error is 4.29%. The result of employing regression equation to predict ultrasonic grinding force approximates to the experimental data, therefore the accuracy and reliability of the model is verified.
基金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.