1.一种模拟气井动态积液和产量变化的方法，其特征在于，所述方法包括以下步骤：S01根据生产历史或试采数据计算及拟合得到产能方程；S02根据油管中液柱和流体的压力分布曲线诊断井是否积液，所述两条压力分布曲线若无交点则无积液，若有交点则有积液，交点对应的深度值为井筒积液液面深度，所述流体压力分布曲线为气井井筒沿程压力分布曲线；S03对比当前气井产量，采用试算法逐步增加或减小产液量，直至当前气井产量与计算得到的标况下的临界携液气量一致，对应的试算法拟定的产液量为当前条件下气井的临界携液量；S04根据气井的临界携液量和从地层流入井底的液量计算当前积在井筒中的液量，根据当前积在井筒中的液量和油管内径计算得到下一时间步长井筒积液高度增量，进而计算下一时间步长井筒积液高度，根据当前时间步的井底流压计算出下一时间步长的井底流压增量，进而计算下一时间步长的井底流压；S05根据物质平衡方程、实际油藏建模结果或历史拟合经验公式计算，得到下一时间步长地层压力；S06根据地层压力和井底流压，通过产能方程计算当前步下的地层产气量；S07重复S02～S05，获取气井积液高度、产气量、地层压力和井底流压随时间变化的数值，实现动态预测。
2.根据权利要求1所述的方法，其特征在于，所述数据包括储层数据、流体物性参数、气井生产数据和井身结构数据。
3.根据权利要求1所述的方法，其特征在于，所述油管中液柱的压力分布以管鞋处压力为初值向上采用静液柱压降计算模型计算得出。
4.根据权利要求3所述的方法，其特征在于，根据所述产能方程计算结果，沿油套环空向下根据静气柱压降计算模型计算得出所述管鞋处压力。
5.根据权利要求1所述的方法，其特征在于，所述油管中流体的压力分布以井口油压为初值采用气液两相流动方程向下计算得出。
6.根据权利要求1所述的方法，其特征在于，所述气井临界携液量基于液膜模型计算获得，所述液膜模型是根据环状流界面剪切应力无量纲表达式计算得到垂直管道界面剪切应力与无量纲液膜厚度的关系，随所述无量纲液膜厚度变化的界面剪切应力存在最小值，为积液临界剪切应力值。
7.根据权利要求1所述的方法，其特征在于，所述气井临界携液量的计算还需要计算气液界面摩擦因子、临界携液流速和临界携液气量。
8.根据权利要求1所述的方法，其特征在于，所述下一时间步长气井积液高度的计算根据气量和气液比计算出当前从地层流入井底的液量。
9.根据权利要求8所述的方法，其特征在于，所述下一时间步气井积液高度为井筒积液初始高度加上下一时间步长井筒积液高度增量。
10.根据权利要求1所述的方法，其特征在于，所述下一时间步长井筒积液高度增量计算中当前积在井筒中的液量为当前条件下气井的临界携液量减去当前从地层流入井底的液量。
11.根据权利要求1所述的方法，其特征在于，所述下一时间步长的井底流压增量，为井筒初始液面高度时的压力值与下一时间步长井筒积液高度时的压力值之差。
12.根据权利要求1所述的方法，其特征在于，所述下一时间步长的井底流压为上一时间步长的井底流压与下一时间步长的井底流压增量的和。
13.根据权利要求1所述的方法，其特征在于，所述物质平衡方程根据累计产气量、水侵量、累计产水量、天然气体积系数、原始状态下天然气体积系统和水的体积系数来计算动态储量。
14.根据权利要求1所述的方法，其特征在于，所述地层压力通过历年的累产气量和地层压力建立物质平衡方程，结合气井目前的累产气量，推算出目前的地层压力。
15.根据权利要求1所述的方法，其特征在于，所述产能方程根据拟合系数、地层压力和井底流压计算地层产气量。
