地震资料反褶积、波阻抗反演是提高分辨率和储层预测识别的关键技术。反褶积的基本原理，是利用采集到的地震数据和已知的子波经过反褶积得到反射系数;波阻抗反演则是由反射系数计算出反应地层物性变化的物理参数波阻抗。　　 为了提高分辨率、计算速度，针对波阻抗反演的稀疏优化问题，根据凸优化理论，本论文提出了基于原始-对偶内点法的波阻抗稀疏反演最优化方法。将波阻抗稀疏反演的问题表示为不等式约束的最优化问题，可以用原始-对偶内点法寻找最优解。内点法分为内外两层循环，外循环使用逐步逼近(SequentialUnconstrainedMinimizationTechnique，SUMT)的最优化方法，以对偶间隙(DualityGap)作为停机准则;内循环为牛顿法解无约束最小化问题，求解当前对偶间隙下的最优解。利用该方法得出的结果提高了分辨率，同时与原始数据的拟合程度没有显著的降低，保持了较高的信噪比。另外，本论文采用了预条件后的最小二乘QR分解(LeastSquaresQR-factorization，LSQR)方法求解线性方程组，精度高，收敛快，占用内存较其它方法小，提高了计算速度。