1.一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：(a)基础曲线预处理；

(b)地震合成记录，建立时深关系曲线；

(c)曲线频谱分析及滤波处理；

(d)多频率曲线融合；

(e)基于频率融合曲线的地震反演；

其中，所述步骤(a)包括如下内容：

(a1)对声波时差曲线开展标准化处理，选取全区稳定发育的泥岩段作为标准层，计算标准井标准层段的测井曲线平均值 同时计算每口井目的层段测井曲线平均值求出校正差，从测井曲线中减去校正量，即完成校正，校正公式为：其中，AC标为标准化后测井曲线值；AC为标准化前测井曲线值； 为目的层段测井曲线平均值； 为标准井标准层段的测井曲线平均值；

(a2)对自然电位曲线开展基线校正，选取全区稳定发育、测井响应具有一致性的稳定的泥岩段作为标准层，读取标准层自然电位曲线平均值作为自然电位曲线基线漂移的基值，选取目的层内各个泥岩层深度及其对自然电位幅值，按照自然电位曲线采样率进行垂向上线性插值，生成校正基线曲线SP基，对自然电位曲线进行基线校正，利用公式：其中， 为第k个采样点处基线校正后的自然电位幅度值；Spk为第k个采样点处基线校正前的自然电位幅度值； 为第k个采样点处校正基准曲线幅度值；SP标准为标准泥岩段自然电位平均幅值；

(a3)将步骤(a2)校正结果SPNEW曲线归一化到波阻抗曲线量纲范围内，根据读取波阻抗曲线开展频率直方图统计，得出波阻抗曲线的最大幅值 最小幅值 由于自然电位幅值极性与波阻抗曲线幅值极性相反，对步骤(a2)校正结果SPNEW曲线开展反极性求取SP-NEW＝SPNEW，并对SP-NEW曲线进行频率直方图统计，统计SP-NEW曲线最大值 和最小值将SPNEW曲线量纲到波阻抗曲线的范围之内，均一化公式如下：其中， 为第k个采样点处均一化后的自然电位幅度值； 为第k个采样点处均一化前的自然电位幅度值； 和 为波阻抗曲线的最大幅值和最小幅值； 和为井待校正自然电位曲线的最大幅值和最小幅值；

(a4)引用步骤(a3)所读取的波阻抗量曲线最大幅值 和最小幅值 对步骤(a2)校正结果SPNEW曲线进行频率直方图统计，得出电阻率曲线最大值 和最小值将电阻率曲线量纲到波阻抗曲线的范围之内，均一化公式如下：其中， 为第k个采样点处均一化后的电阻率幅度值； 为第k个采样点处均一化前的电阻率幅度值； 和 为波阻抗曲线的最大幅值和最小幅值； 和为井待校正电阻率曲线的最大幅值和最小幅值；

通过步骤(a3)和步骤(a4)将完成自然电位曲线和电阻率曲线的归一化处理，将自然电位曲线和电阻率曲线归一到波阻抗曲线量纲范围，以此达到各井曲线的一致性和可对比性；

所述步骤(b)包括如下内容：

(b1)制作地震合成记录，建立时深关系曲线，利用步骤(a1)标准化的声波时差曲线和密度曲线，根据公式：k k k

其中，R 为第k个采样点处换算反射系数值；ρ为第k个采样点处密度曲线幅度值；v 为第k个采样点处速度曲线幅值，其大小为换算出反射系数后，并选用合适频率的雷克子波δ，将反射系数曲线与子波褶积得到合成记录：其中，Xi为第i个采样点合成记录值；δk为第k个采样点处子波值；Ri为第i个采样点处反射系数；M1、M2分别为子波的起始样点和末尾样点序号；

以合成记录为基础，建立地震数据的时间域与测井深度域之间的对应关系；

(b2)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到SPlog曲线转换为地震反射时间与自然电位幅值相对应的SPtime曲线；

(b3)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到RLLDlog曲线转换为地震反射时间与电阻率幅值相对应的RLLDtime曲线；

所述步骤(c)包括如下内容：

(c1)对步骤(b1)中的得到SPtime曲线进行快速傅里叶变换(FFT)，得到SPtime曲线频率谱H(f)SP，利用低通滤波器公式：H(f)SPlow＝C(f)×H(f)SP

其中，C(f)为低通滤波函数；f为频率；fL为低通滤波器的截止频率；H(f)SP为自然电位曲线频谱函数；H(f)SPlow为自然电位曲线低频分量频谱函数；

通过低通滤波处理得到自然电位曲线低频分量频谱函数H(f)SPlow；

(c2)对步骤(b2)中的得到RLLDlog曲线进行快速傅里叶变换(FFT)得到RLLDtime曲线频谱H(f)RLLD，根据高通滤波器公式：H(f)RLLDhigh＝C(f)×H(f)RLLD

其中，C(f)为高通滤波函数；f为频率；fH为高通滤波器的截止频率；H(f)RLLD为自然电位曲线频谱函数；H(f)RLLDhigh为电阻率曲线高频分量频谱函数；

通过高通滤波处理得到电阻率曲线高频分量频谱函数H(f)RLLDhigh；

截止频率fL和fH的选取为地震反演所应用到的地震数据的有效频带的最高值；

所述步骤(d)包括如下内容：

通过融合技术将步骤(c1)得到的自然电位曲线低频分量频谱函数H(f)SPlow和步骤(c2)得到的电阻率曲线高频分量频谱函数H(f)RLLDhigh进行融合生成频谱曲线H(f)Ip；针对频谱曲线H(f)Ip开展反傅里叶变换，即得到基于频率融合拟波阻抗曲线Ip拟合；

所述步骤(e)包括：将步骤(d)所得频率融合拟波阻抗曲线Ip拟合应用于地震反演运算过程。