Curvature estimates for 4-dimensional complete gradient expanding Ricci solitons