System and method for determining a frequency sweep for seismic analysis

A method for calculating a target output spectrum to be achieved by a vibrator device, the method comprising:

• receiving constraint information related to the vibrator device;
• using a mathematical relationship for achieving the target output spectrum, wherein the mathematical relationship relates the (i) the constraint information of the vibrator device, (ii) the target output spectrum, and (iii) a frequency sweep rate;
• and
• determining, based on the mathematical relationship, a frequency sweep for achieving the target output spectrum by the vibrator device in compliance with the constraint information of the vibrator device,
• wherein said receiving, said using, and said determining are performed by a processor-based device, and
• wherein the mathematical relationship is given by:
• SR(f)=β·[Fconstraint(f)/Ftarget(f)]2·[F1−F0]/SL,
• where f is an instantaneous frequency, SL is a sweep length with a start frequency F0 and an end frequency F1, Ftarget(f) is the target output spectrum, Fconstraint(f) is a force constraint related to the constraint information, SR(f) is the frequency sweep rate, and β is a constant.

A computer device for generating a target output spectrum for use in seismic exploration to be achieved by a vibrator device, the computer device comprising:

• an interface configured to receive constraint information related to the hydraulic seismic vibrator device;
• and
• a processor connected to the interface and configured to,
• use a mathematical relationship for achieving the target output spectrum, wherein the mathematical relationship relates the (i) the constraint information of the vibrator device, (ii) the target output spectrum, and (iii) a frequency sweep rate;
• and
• determine, based on the mathematical relationship, a frequency sweep for achieving the target output spectrum by the vibrator device in compliance with the constraint information of the vibrator device,
• wherein the mathematical relationship is given by,
• SR(f)=β·[Fconstraint(f)/Ftarget(f)]2·[F1−F0]/SL,
• where f is an instantaneous frequency, SL is a sweep length with a start frequency F0 and an end frequency F1, Ftarget(f) is the target output spectrum, Fconstraint(f) is a force constraint related to the constraint information, SR(f) is the frequency sweep rate, and β is a constant.