In this paper we collect the isoseismal data of 180 earthquakes from 1990 to 2021 in eastern China,adopt the long- and short-axis ellipse model,and apply the Levenberg-Marquardt method and the least squares method to re-fit the regional seismic intensity attenuation relationship in the study region. The results show that the two regression models maintain a good consistency,and the seismic intensity attenuation relationships in the two regional statistical units(Northeast China and North China and Central China and South China) are significantly different,so different partitions should be treated. Moreover,the comparison with the previous results of earthquake intensity attenuation relationship shows that the results of the two methods in this paper better reflect the regional characteristics of intensity attenuation. Taking long,short axis and area calculations of the Ningjiang M5.7 earthquake in Jilin,the Badong M5.1 earthquake in Hubei and the M7.8 earthquake in Hebei Tangshan,we found that the calculation results of the intensity attenuation relationship in this paper are slightly more accurate than that of previous studies. So that the results of this paper can reflect the characteristics of small earthquakes and large damages in the eastern region,and are suitable for engineering earthquake research and applications in this region. The applicable magnitude range of the regional earthquake intensity attenuation relationship in eastern China given in this paper is M3.0～8.0,with the distance of 0～660km.