In order to assure structural integrity for operating welded structures, it is necessary to evaluate crack growth rate and crack propagation direction for each observed crack non-destructively. Here, three dimensional (3D) welding residual stresses must be evaluated to predict crack propagation. Today, X-ray diffraction is used and the ultrasonic method has been proposed as non-destructive method to measure residual stresses. However, it is impossible to determine residual stress distributions in the thickness direction. Although residual stresses through a depth of several tens of millimeters can be evaluated non-destructively by neutron diffraction, it cannot be used as an on-site measurement technique. This is because neutron diffraction is only available in special irradiation facilities. Author pays attention to the bead flush method based on the eigenstrain methodology. In this method, 3D welding residual stresses are calculated by an elastic Finite Element Method (FEM) analysis from eigenstrains which are evaluated by an inverse analysis from released strains by strain gauges in the removal of the reinforcement of the weld. Here, the removal of the excess metal can be regarded as non-destructive treatment because toe of weld which may become crack starters can be eliminated. The effectiveness of the method has been proven for welded plates and pipes even with relatively lower bead height. In actual measurements, stress evaluation accuracy becomes poorer because measured values of strain gauges are affected by processing strains on the machined surface. In the previous studies, the author has developed the bead flush method that is free from the influence of the affecting strains by using residual strains on surface by X-ray diffraction. However, stress evaluation accuracy is not good enough because of relatively poor measurement accuracy of X-ray diffraction. In this study, a method to improve the estimation accuracy of residual stresses in this method is formulated, and it is shown numerically that inner welding residual stresses can be estimated accurately from the residual strains measured by X-ray diffraction.