1.一种气藏状态参数随钻数值反演分析方法，其特征在于，包括以下步骤：
1)建立控制方程，包括如下步骤：
11)建立单液相井筒压力场模型；公式如下：
式中： 为总压力梯度，Pa/m；
Zv为垂深，m；
为单液相钻井液摩阻压力梯度，Pa/m；
采用Herschel-Bulkley流变模型计算液体钻井液摩阻压力梯度；
12)建立单液相井筒换热模型，包括如下步骤：
121)建立钻柱内钻井液模型
122)建立钻柱管体模型
123)建立环空钻井液模型
124)建立地层模型
125)建立井壁与地层界面模型
13)建立气液两相流换热模型，包括如下步骤：
131)建立质量守恒方程包括：
建立钻井液质量守恒方程：
建立游离气质量守恒方程：
建立溶解气质量守恒方程：
132)建立动量守恒方程，如下式：
133)建立能量守恒方程，包括：
建立钻柱内能量守恒方程：
建立环空能量守恒方程：
建立井筒/地层交界面能量守恒方程：
式中f(t)为无因次瞬时热损失函数；
14)建立气藏渗流模型，包括以下步骤：
141)建立低压气体非达西不稳定渗流模型：
142)建立高压气体非达西不稳定渗流关系模型：
15)建立天然气PVT模型，公式如下：
其中：
ρr＝0.27pr/(Z.Tr)，pr＝p/pc,Tr＝T/Tc，Z＝pv/(RT),ρ＝1/v
其它参数如下：
A1＝0.3265,A2＝-1.0700,A3＝-0.5339,A4＝0.01569
A5＝-0.05165,A6＝0.5475,A7＝-0.7361,A8＝0.1844
A9＝0.1056,A10＝0.6134,A11＝0.7210；
步骤16)建立天然气粘度模型，公式如下：
R1:0<ρg≤0.26g/cm3,0.7≤Tr≤3.4,16≤Mw≤72.2
μ＝5.6563271×10-2ρg+4.9374602×10-3Tr+4.1949307×10-5Mw+2.93978342×10-3R2:0.26<ρg≤0.46g/cm3,0.88≤Tr≤3.4,16≤Mw≤72.2
μ＝1.3707401×10-1ρg+2.72296913×10-3Tr+7.64990184×10-6Mw+4.050623771×10-3R3:0.46<ρg≤0.595g/cm3,0.645≤Tr≤1.5,21.6≤Mw≤119
μ＝4.8834347×10-1ρg+4.437225271×10-2Tr+6.5756117×10-4Mw-2.002744453×10-1R4:0.595<ρg≤0.76g/cm3,0.5≤Tr≤1.4,35≤Mw≤119
μ＝1.500143567ρg+9.84848797×10-2Tr+1.407003797×10-3Mw-8.257944449×10-1  (16)
17)建立溶解速率模型，公式如下：
式中： 为考虑三种主要组分的饱和溶解度；
cbulb为平均溶解度；
K为质量传输系数，层流时：k＝0.026Re0.8Sc1/3D/(dout-dinner)，紊流时：
18)建立气液两相流流速模型，公式如下：
vg＝C0vm+v∞   (18)
式中：
C0——速度剖面系数；两相流中心处速度与平均流速的比值，取决于流型，流型包括：泡状流、段塞流、搅动流、环状流；
vm——气液两相流混合物平均速度；
v∞——气相的相对液相滑移速度，取决于流型；
2)求解域网格划分与节点布置，应用交错网格技术布置温度和压力速度节点位置，其中温度节点位于网格中心，压力和速度项布置于网格边界处，从而保证离散变量严格满足控制方程和求解变量的连续；
轴向坐标自井口指向井底，轴向网格编号自井口向井底逐渐增加，井口节点编号为1，井底增加一个网格轴向长度为0的虚拟节点N+1，各离散区域相同节点编号的网格轴向位置相同，过流面积改变处布置于网格边界处；
油气扩散运移模拟过程中，新时间步长内，气相运移距离不允许超过网格长度，否则会造成计算不收敛问题，即新时间步的步长应满足：
3)控制方程离散，包括如下步骤：
31)质量守恒方程离散
离散质量守恒方程7、方程8和方程9得到：
液相：
自由气相：
溶解气相：
n-1 n
t 时刻网格参数和上游(i+1)节点t时刻参数已知，故，质量守恒离散方程求解实质为tn时刻网格下游边界处液相、溶解气和自由气的质量传输项求取；环空中单元格按照内部流体相组成分为三种网格单元：
①气液两相混合网格单元；
②气相前缘网格，该网格内一部分为气液两相，一部分为纯液相；
③纯液相网格；在气相前缘所在网格上部的网格为纯液相网格；
311)气液两相网格下游边界质量传输
气液两相网格边界处流出的三种物质质量，可以根据时间步长、上一时间步长的体积含气率、上一时间步长各自的密度、和当前时间的各自速度计算得出；
液相：
自由气：
溶解气：
312)纯液相网格下游边界质量传输
纯液相网格边界处流出的只有不含溶解气的液相，可以根据时间步长和当前时间的速度计算得出；
液相：
自由气：
溶解气：
313)气相前缘网格下游边界质量传输
tn-1时刻到tn时刻下游边界处气相速度为 液相速度为 纯液相中不含溶解气；
tn-1时刻节点i中气液两相部分长为Lg，纯液相部分长度为(△Zi-Lg),气液两相部分的体积含气率为 液相密度为 气相密度为 纯液相密度为ρι；
气相前缘经过
后到达网格i下游边界处，△t时刻内纯液相速度为网格气液两相平均速度，即①Δtn≤Δt时，即
Δtn时段气相前缘未超过下游边界，这种情况下排出的只有纯液相，排出的质量为：
液相：
溶解气：
自由气：
②Δtn>Δt时
Δtn时段气相前缘超过下游边界，根据排出的物质是否含有溶解气分两种情况确定排除物质的质量；
a.
Δtn时段排出的液相体积不大于tn-1时刻网格内的纯液相体积，不排出溶解气；
液相：
溶解气：
自由气：
b.
Δtn时段排出的液相体积大于tn-1时刻网格内的纯液相体积，排出液相中包含有溶解气；
液相：
溶解气：
自由气：
式中 为i网格tn-1时刻含溶解气的液相密度，ρL为不含溶解气的液相密度；
314)网格溶解气质量求取
Δt时间段溶解气质量求取区域为网格上游边界到气相前缘，溶解速率按照n-1时刻网格上游到网格气相前缘区域内的物质构成和流动特征进行求取，忽略该时间段内地层侵入的气体产生的溶解气；
Δt时间后静止网格内的溶解气质量为：
该方程为标量非线性方程，迭代求解出 后，则得Δt时间段内网格溶解气质量增加量为：
32)动量守恒方程离散
环空钻井液向上流动，轴坐标方向自上而下增加，二者方向相反，离散时将速度项取正值，得离散格式的动量守恒方程：
含有气相前缘的网格中，对单液相区和气液两相区分别求取摩阻压力降，通过加权求和得该网格的总压力降；
33)能量守恒方程离散
应用比焓微商热力学理论，消去井筒/地层界面温度，得到用热阻表示的环空节点温度计算离散方程：
式中， 为环空钻井液到地层的总热阻；
4)气藏参数数值反演算法，包括如下步骤：
41)划分网格
根据井身结构和钻具结构划分网格，设置节点数组pNode[N]和pNodeNew[N]；pNode中存放上一时间点网格数据，pNodeNew存放新时间点数据；
42)计算网格垂向坐标，更新pNode节点数组；
根据井眼轨迹计算网格中心的垂直坐标和网格垂直长度；
43)按照单液相流动进行温度场分析，更新pNode节点数组
44)按照单液相进行网格流动参数和压力计算并；包括液相质量、流出质量、流出质量流量、液相密度、液相流速、压力梯度、压力、截面含气率赋0、气相流速赋0、气相流出质量赋
0、溶解气质量赋0、流出溶解气质量赋0，更新pNode节点数组；
45)确定气侵节点编号M
遍历网格，比对网格节点轴向位置和气层深度，确定气侵节点编号M；
46)扩散运移模拟初始条件应用
扩散运移模拟时间t＝0，气相前沿网格索引 气相前沿网格内气相长度 (点
气侵处理)；
将网格变量pNodeNew[N]同步为pNode[N]数值；
47)根据网格长度和网格气相速度确定时间步长dt；
48)时间推进，t＝t+dt，
49)气侵上游节点压力赋初值
气侵上游节点压力调整值pDiff＝0，
410)确定气侵上游节点压力，
411)根据气层参数、钻进参数和 计算气侵量
412)根据网格类型计算网格参数
自下而上遍历网格，根据网格索引号i确定网格类型，计算网格流动参数；
网格内为单液相流动
气相前沿网格，更新 和
气液两相流网格
413)比较井口节点压力与大气压，计算相对误差εp，确定底部边界压力调整值
414)若εp>10-3，转至步骤410)，重复执行，直至εp<10-3，新时间点网格参数计算结束；
415)根据新时刻网格流动参数和压力分析新温度场；
416)更新全部环空网格数据pNode为pNodeNew；
417)若时间到达分析时间，则分析结束，否则，转至步骤47)重复执行。
2.根据权利要求1所述的一种气藏状态参数随钻数值反演分析方法，其特征在于：应用气藏状态参数随钻数值反演解释方法编制气藏状态参数分析模拟器，包括功能模块如下：
ReadKickData
ReadKickData是分析模拟器的数据输入函数，完成模拟井所有数据的输入，主要包括井身结构，钻具结构，井眼轨迹，地层垂直温度分布，储层参数：温度、压力、孔隙度、渗透率、深度、供给半径，原油和天然气热物性参数，水泥、地层、钻柱等的热物性参数，钻进参数：排量、机械钻速、扭矩、钻头喷嘴，模拟参数：轴向网格尺度、模拟周期、是否考虑破碎气、是否考虑热源；
KickFlowBehavior
KickFlowBehavior函数是模拟器的总功能模块，按照程序流程图组装其它函数，完成钻井油气扩散及运移模拟分析及数据存储功能；
TemperatureAnalysis
TemperatureAnalysis函数根据网格流动参数进行钻井井筒温度场分析，将温度计算结果更新到网格节点数组pNodeNew中；
KickSimulating
KickSimulating函数根绝网格节点温度和井底与井口边界条件进行钻井油气扩散和运移模拟分析，并将流动参数计算结果更新到网格节点数组pNodeNew中；
GridGeneration
GridGeneration函数根据井身结构和钻柱结构对温度和压力场求解域进行轴向分段和网格划分，存储网格节点的轴向几何信息及径向换热对象的几何和介质信息；
NormalTPNode
NormalTPNode函数根据入口和出口都是气液两相流的网格流动参数计算方法，分析网格流动参数；
GasFrontNode
GasFrontNode函数根据气相前沿网格算法，分析网格流动参数，并更新气相前沿网格索引和气相区长度数据；
VerticalCoordinate
根据井眼轨迹计算网格中心的垂直坐标和网格垂直长度；函数功能要求调用
DirectionParaCal函数根据井深和井斜角计算垂深；
TOriginGeneration
TOriginGeneration函数根据地层的垂直温度分布数据插值产生节点垂直深度处的原始温度；
Ini
Ini函数对求解域内网格节点变量应用初始条件，赋初值；
HeatResistanceSA
计算钻主内钻井液与环空钻井液的换热热阻；实现该函数功能需要调用HPipe和HAnnu函数计算液相在钻主内和气液两相流在环空中的强迫对流换热系数；
HeatResistanceAE
计算环空气液两相流钻井液与地层之间的换热热阻；实现该函数功能需要调用HAnnu函数计算气液两相流在环空中的强迫对流换热系数和瞬时地层热损失函数计算井筒到地层的换热热阻；
HPipe
HPipe函数的功能是计算钻主内钻井液与钻柱内壁的管流强迫对流换热系数；
HAnnu
HAnnu函数的功能是计算环空气液两相流钻井液与井壁的环空流强迫对流换热系数；
PPTP
根据介质类型计算介质在给定温度和压力下的热力学性质的总集成函数，需要根据具体的介质类型原油、水、天然气；调用相应介质的热物性计算函数完成其功能；
DPDZ_TwoPhase
根据气液两相流机理模型计算气液两相流流动参数和压降梯度，气液两相流流动参数包括：气相速度、液相速度、截面含气率、气相密度、流型；
TInDrillstem
TInDrillstem函数根据单液相能量守恒方程迭代计算钻主内钻井液节点温度；
TAnnulus
TAnnulus函数根据气液两相流能量守恒方程迭代计算环空钻井液节点温度；
CpGas，计算天然气定压比热；
CpOil，计算原油定压比热；
CpWater，计算地层水定压比热；
DenGas，计算天然气密度；
DenOilSaturated，计算饱和原油密度；
DenOilUnsaturated，计算不饱和原油密度；
RsStanding，计算原油溶解气油比；
SurfaceTension，计算油水界面张力；
ThconGas，计算天然气热导率；
ThconOil，计算原油热导率；
ThconWater，计算地层水热导率；
ViscGas，计算天然气粘度；
ViscOil，计算原油粘度；
ViscWater，计算地层水粘度。
3.根据权利要求2所述的一种气藏状态参数随钻数值反演分析方法，其特征在于：气藏状态参数随钻数值反演解释方法的边界条件与初始条件设置：
边界条件
①钻柱入口处：
Qm＝C  (45-1)
T1(z＝0)＝TE  (45-2)
其中TE为钻井液入口温度，Qm为钻井液排量；
②环空出口处：
p＝patm  (46)
③井底钻柱与环空连接处，钻井液温度相同，有
初始条件
①温度场的初始条件为单液相条件下已经计算得到的井筒温度，
即：
Ti(z,0)＝Ti,st(z)  (48)
②压力场的初始条件为单液相条件下已经计算得到的井筒压力，
即：
pi(z,0)＝pi,st(z)  (49)
③气侵开始时，环空各网格内处无气相：
αi(z,0)＝0  (50)。