Conventional full-waveform inversion (FWI) based on the least-squares misfit function faces problems in converging to the global minimum when using gradient methods because of the cycle-skipping phenomena. An initial model producing data that are at most a half-cycle away from the observed data is needed for convergence to the global minimum. Low frequencies are helpful in updating low-wavenumber components of the velocity model to avoid cycle skipping. However, low enough frequencies are usually unavailable in field cases. The multiplication of wavefields of slightly different frequencies adds artificial low-frequency components in the data, which can be used for FWI to generate a convergent result and avoid cycle skipping. We generalize this process by multiplying the wavefield with itself and then applying a smoothing operator to the multiplied wavefield or its square to derive the nonlinearly smoothed wavefield, which is rich in low frequencies. The global correlation-norm-based objective function can mitigate the dependence on the amplitude information of the nonlinearly smoothed wavefield. Therefore, we have evaluated the use of this objective function when using the nonlinearly smoothed wavefield. The proposed objective function has much larger convexity than the conventional objective functions. We calculate the gradient of the objective function using the adjoint-state technique, which is similar to that of the conventional FWI except for the adjoint source. We progressively reduce the smoothing width applied to the nonlinear wavefield to naturally adopt the multiscale strategy. Using examples on the Marmousi 2 model, we determine that the proposed FWI helps to generate convergent results without the need for low-frequency information.