专利名称：一种地层回波信号提取方法及装置的制作方法
技术领域：
本发明涉及ー种井下核磁共振测井仪的地层回波信号的提取技木，尤其涉及ー种地层回波信号提取方法及装置。
背景技术：
井下核磁共振仪器测井过程中，由于地层回波信号非常微弱，使得回波信号的提取十分困难。现有核磁共振测井仪中普遍使用数字相敏检波器(DPSD)的检波算法或者小波提取算法，对输入信号进行幅度与相位的提取。但是目前的提取方法误差较大，精度较低，且极易受到外界噪声的干扰影响回波信号的提取质量。

发明内容
本发明所要解决的技术问题是克服目前提取核磁共振的地层回波信号误差较大、精度较低的缺陷。 为了解决上述技术问题，本发明提供了一种地层回波信号提取方法，包括如下步骤采集核磁共振的回波原始信号；将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量；对所述正交分量进行小波分解和小波重构，获得正交输出信号；对所述同相分量进行小波分解和小波重构，获得同相输出信号；采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。优选地，对所述正交分量进行小波分解和小波重构的步骤，包括对所述正交分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构；对所述同相分量进行小波分解和小波重构的步骤，包括对所述同相分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。优选地，对所述正交分量或同相分量以db5小波基进行七层小波分解的步骤中，每ー层小波分解均采用如下表达式进行d) = d\' + a(s\1 + )；s) ニ S1' + j3yd) +；= d) + /(S'l + (I)；= S11 + SlCi^z + §2d;n ；St=i.S;-;以及dt=d;/^；
其中，第一层小波分解中，<和ぜ分别为对所述正交分量或同相分量进行奇偶分裂操作得到的奇地址数据和偶地址数据；第二至第七层小波分解中，<和^*分别为对上ー层小波分解得到的dt进行奇偶分裂操作得到的奇地址数据和偶地址数据；<+1表示<向左移动ー个时间单位后得到的数据；表示式向右平移一个时间単位后得到的数据；表示ゼ向右平移一个时间単位后得到的数据；尤2表示ゴム向左延迟ー个时间周期后得到的数据；ゴム表示<向左平移一个时间周期后得到的数据；
a , ^ ^ 2> y , 6 6 2 以及 2,为系数。优选地，a= -I. 5861343420693648 ； ^ : = I. 0796367753628087 ； ^ 2=-0. 0529801185718856 ； y = -0. 8829110755411875 ^1 = O. 4435068520511142 ； S2 =
I.576237461483639 ;以及 I = -I. 1496043988602418。本发明提供了一种地层回波信号提取装置，包括采集模块，用于采集核磁共振的回波原始信号；运算模块，用于将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量；分解重构模块，用于对所述正交分量进行小波分解和小波重构，获得正交输出信号；用于对所述同相分量进行小波分解和小波重构，获得同相输出信号；合成模块，用于采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。优选地，所述分解重构模块用于对所述正交分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构；所述分解重构模块用于对所述同相分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。优选地，所述分解重构模块用于采用如下表达式进行七层小波分解中的每ー层小波分解d) = 4° + び(< + <+1)；S1t = sf + pyd\ + fi2d)_y ；=d) + y{s]+s]^)\sf = sj + + S2d;+l ；ち=#.42;以及dt=df ;其中，第一层小波分解中，<和sf分别为对所述正交分量或同相分量进行奇偶分裂操作得到的奇地址数据和偶地址数据；第二至第七层小波分解中，<和 < 分别为对上ー层小波分解得到的dt进行奇偶分裂操作得到的奇地址数据和偶地址数据；<+1表示^>向左移动ー个时间单位后得到的数据；
表示式向右平移一个时间単位后得到的数据；表示べ向右平移一个时间単位后得到的数据；<+2表示向左延迟一个时间周期后得到的数据；表示<向左平移一个时间周期后得到的数据；a , ^ ^ 2> y , 6 6 2 以及(为系数。优选地，a= -I. 5861343420693648 ； ^ : = I. 0796367753628087 ； ^ 2=-0. 0529801185718856 ； y = -0. 8829110755411875 J1 = O. 4435068520511142 ； S2 =I. 576237461483639 ;以及 I = -I. 1496043988602418。与现有技术相比，本发明的实施例将提升小波算法与DPSD相敏检波算法两者结合在一起应用于核磁共振仪器微弱信号的提取设计，兼容了两种算法各自的优点。通过将 两种算法结合设计来实现对核磁微弱信号的提取，即具有小波算法的优点又同时兼备DPSD算法的优点，可以高效地滤除信号中的噪声干扰，比传统核磁共振测井仪所使用的信号提取算法效率更高，提取有用信号的信噪比更大，大大提高了核磁共振测井中回波信号的提取质量。本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。


附图用来提供对本发明技术方案的进ー步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明的技术方案，并不构成对本发明技术方案的限制。图I为本发明实施例的地层回波信号提取方法的流程示意图。图2a为本发明实施例中对回波信号的正交分量和同相分量进行小波分解和小波重构的流程示意图。图2b为本发明实施例中对回波信号进行小波分解和小波重构的原理示意图。图3为本发明实施例中对回波信号的正交分量和同相分量进行小波分解的分解模块的原理示意图。图4为本发明实施例中进行小波重构的重构模块的原理示意图。图5为本发明实施例的地层回波信号提取装置的结构示意图。
具体实施例方式以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。本发明实施例以及实施例中的各个特征在不相冲突前提下的相互结合，均在本发明的保护范围之内。如图I所示，本发明实施例的地层回波信号提取方法主要包括如下步骤。步骤S110，采集获得核磁共振的回波原始信号x(n)。步骤S120，将回波原始信号x(n)与同频率的正弦信号Sin(CociIi)和余弦信号cos(co0n)分别进行正交乘法运算，获得回波原始信号的正交分量Q(n)和同相分量I(n)。
步骤S130，对正交分量Q(n)进行小波分解和小波重构，获得正交输出信号Q_out(n);对同相分量I(n)进行小波分解和小波重构，获得同相输出信号I_out(n)。步骤S140，采用数字相敏检波技术对正交输出信号Q_out(n)和同相输出信号1_out (n)进行相位合成和幅值合成的处理，获得核磁共振回波信号的初始相位デ(《)和包络信号幅值y(n)。本发明的实施例中，步骤S130中，对回波信号的正交分量Q(n)和同相分量I (n)进行小波分解和小波重构的过程可以采用如图2a所示的流程来进行。步骤S131，将回波信号的正交分量Q(n)和同相分量I (n)分别以db5小波基进行七层小波分解，各自获得七组细节系数和ー组最底层近似系数(第七层分解后的近似系数)。步骤S132，将七组细节系数全部置为零。
步骤S133，利用该最底层近似系数以及置为零的七组细节系数进行小波重构，得到两路原始信号しout (n)与Q_out(n)以完成对其的噪声滤波处理。如图2b所示，将输入信号按照图2b所示进行抽取后，得到的奇偶两组序列与行列式/"7(一 f进行乘法运算，最后得出的值即为分解后的奇地址分解数据dt和偶地址分解数据
StO图2b 中
r n …、「 I Op (A+ル)]「I 0]「1 +0],、./(J=4 ハ,^n - 式
_a(l + z ) I 0 I _ ブ(1 + z) I 0I0 i
(A)P(z~')r是对表达式A中行列式P(Z)进行转置求倒后得出的，保证POT1)7'* P (め=I。其中，a、运” @2、Y、S2以及(等都为系数。z表示时域信号作频域Z变换的基本变量。图2b中的户(ff与P(Z)的意义表示输入的时域信号进行频域z变换后与或者p(z)作行列式频域乘法运算，得出的值在转换到相对应的时域上，即得到图3与图4所示信号分解与重构的时域操作流程图。小波分解过程即是完成将输入信号的奇偶两组序列与户(厂1/进行行列式相乘得出两组分解后的系数。进ー步将其输入信号的奇偶两组序列与POT1,行列式相乘的运算进行数学上的分步拆分化简，最后得出了如图3所示小波分解的步骤与相关计算公式，如图4所示的重构过程也类似与分解过程。将分解后所得的最底层近似系数与置为零的细节系数乘以表达式P(Z)进行重构，将该步行列式乘法运算进行数学的分解化简后，得出如图4所示结构与相关计算公式。上述步骤S131中以db5小波基对正交分量Q (n)和同相分量I (n)进行七层小波分解中，每ー层小波分解都可以采用如图3所示的分解模块来进行。如图3所示，分解模块主要包括信号分离单元310，第一、第二、...、第八延时单元321、322、 ,328,以及第一、第二、 、第六预测更新单元 331、332、 、336。在第一层小波分解中，信号分离单兀(SPLIT) 310用于对输入的正交分量Q(n)或同相分量I (n)进行奇偶分裂操作；按照对正交分量Q (n)或同相分景T (n)的采样信号x(t)的地址编号的奇偶性，将采样信号分为偶地址数据 < 和奇地址数据J
其中=x(/)',式(I)4。=冲):,ト1 式⑵表达式(I)和表达式(2)中的下脚标t为大于等于0的整数。x(t)2t表示偶地址编号的采样信号x(t)的集合，x(t)2t+1表示奇地址编号的采样信号x(t)的集合。偶地址数据<输入到第一延时单元321进行第一延时处理，经过第一延时处理的偶地址数据^)输入到第三延时单元323进行第三延时处理；奇地址数据c/ 入到第二延时単元322进行第二延时处理。经过第一延时处理的偶地址数据^>、经过第三延时处理的偶地址数据4以及经过第二延时处理的奇地址数据ゴf一起输入到第一预测更新単元331。第一预测更新単元331，分别完成预测、更新及归ー化操作。通过如下表达式进行计算，得到第一奇地址预测数据劣。 rf; =<+び(<+A)式⑶表达式(3)中，系数a在本发明的实施例中取值为-I. 5861343420693648 C1表示偶地址数据<向左移动ー个时间单位(超前预测)后得到的数据，比如^序列的当前值是< 序列下ー个时刻值。第一预测更新単元331的输出信号分别输入到第四延时单元324和第二预测更新单元332。第二预测更新单元332接入第三延时单元323的输出、第一预测更新单元331的输出以及第四延时单元324的输出，根据如下表达式进行计算，得到第一偶地址预测数据= S0t + 賦 + ル/ム式(4)表达式(4)中，參数P1在本发明的实施例中取值为1.0796367753628087 ;參数 在本发明的实施例中取值为-0. 0529801185718856 表示劣向右平移ー个时间单位
后得到的数据，比如 < 序列当前值是<4序列下一个时刻的值。第二预测更新単元332的输出信号分别输入到第五延时单元325和第三预测更新単元333。第五延时单元325的输出信号输入到第六延时单元326进行第六延时处理。第三预测更新単元333接入第四延时单元324的输出、第二预测更新単元332的输出以及第五延时单元325的输出，根据如下表达式进行计算，得到第二奇地址预测数据42。42 = 41 + YisI + 4—1 试(5)表达式(5)中，參数Y在本发明的实施例中取值为-0. 8829110755411875 表示4向右平移一个时间単位后得到的数据，比如く i序列对应的值是べ序列对应上一个时间时刻的值。第三预测更新単元333的输出信号分别输入到第七延时单元327和第四预测更新単元334。第七延时单元327的输出信号输入到第八延时单元328进行第八延时处理。第四预测更新単元334接入第六延时单元326的输出、第三预测更新単元333的输出以及第七延时单元327的输出，根据如下表达式进行计算，得到第二偶地址预测数据e2
°sf = s) + S'dl2 + S2d~t ェ式(6)表达式(6)中，參数S1在本发明的实施例中取值为o. 4435068520511142 ;參数 S2在本发明的实施例中取值为I. 576237461483639。ゴム表示Jw向左延迟一个时间周期后得到的数据，当前4序列对应的值是Ai1序列对应下ー时刻的值；式^表示或2向左平移ー个时间周期后得到的数据，当前<+1序列对应的值是 < 序列对应下ー时刻的值。第四预测更新単元334将第二偶地址预测数据ろ2输入到第五预测更新単元335，第五预测更新単元335根据如下表达式(7)进行计算，输出一次小波分解后的偶地址分解数据st。第八延时单元328的输出信号输入到第六预测更新単元336，第六预测更新単元336根据如下表达式(8)进行计算，输出一次小波分解后的奇地址分解数据dt。ろ=彳式(7)4=42 K式⑶表达式(7)和表达式⑶中，參数I在本发明的实施例中取值为-I. 1496043988602418。将上述的奇地址分解数据dt再次进行信号采样形成采样信号X (t)，输入到图3所示的分解模块继续进行分解操作，后续每ー层小波分解的输入信号，均是上ー层小波分解得到的奇地址分解数据。直到对正交分量Q(n)和同相分量I(n)进行七层小波分解，最后得到最底层近似系数。本发明的实施例中，各个预测更新单元采用位移的方式来实现乘法运算，如此可以高效地降低资源量，提高运行速度，方便硬件实现。步骤S133的小波重构过程是步骤S131的小波分解过程的逆过程，可以完成对信号的重构操作。上述步骤S133中，采用如图4所示的重构模块来进行小波重构，得到两路原始信号 I_out (n)与 Q_out(n)。经过七层分解后，得到的奇地址分解数据dt和偶地址分解数据St分别输入到如图4所示的第一预测更新単元411和第二预测更新単元412进行小波重构运算。此处的小波重构运算是前述小波分解运算的逆运算。在这个逆运算的处理过程中，使用前述小波分解所得到的最底层近似系数，以及置为零的各细节系数。第一预测更新単元411的输出信号<输入到第一延时单元421以及第三预测更新単元413，第一延时单元421输出的信号也输入到第三预测更新単元413。第二预测更新单元412的输出信号#输入到第三预测更新単元413。第一延时单元421的输出信号输入到第三延时单元423，第三延时单元423的输出信号输入到第四预测更新単元414。第三预测更新単元413产生数据^分别发送给第四延时单元424和第四预测更新単元414，第四延时单元424的输出信号也输入到第四预测更新单元414。第四预测更新単元414产生数据J分别发送给第五延时单元425和第五预测更新单兀415。第四延时单兀424的输出信号以及第五延时单兀425的输出信号也输入到第五预测更新単元415。第五延时单元425的输出信号还输入到第六延时单元426。
第五预测更新単元415产生数据<分别发送给第七延时单元427和第六预测更新単元416，第四延时单元424的输出信号也输入到第六预测更新単元416，第七延时单元427的输出信号在输入到第八延时单元428的同时也输入到第六预测更新単元416。第六预测更新単元416产生数据<以及第八延时单元428输出的数据^一起输入到信号合成単元430中产生合成信号y (t)，据此可以获得对应的原始信号I_out(n)与Q_
out (nノ。本发明的实施例通过上述小波分解与重构计算处理后的正交输出信号Q_out(n)和同相输出信号I_out(n)，从回波原始信号中滤除了混杂在其中的噪声信号，大大提高了输出信号的信噪比。如图5所示，本发明实施例的地层回波信号提取装置，主要包括采集模块510、运算模块520、分解重构模块530以及合成模块540，其中采集模块510，用于采集核磁共振的回波原始信号；运算模块520，与采集模块510相连，用于将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量；分解重构模块530，与运算模块520相连，用于对所述正交分量进行小波分解和小波重构，获得正交输出信号；用于对所述同相分量进行小波分解和小波重构，获得同相输出信号;合成模块540，与分解重构模块530相连，用于采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。所述分解重构模块530用于对所述正交分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。所述分解重构模块530用于对所述同相分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。所述分解重构模块530用于采用如下表达式进行七层小波分解中的每ー层小波分解d)=《+ a{ sf + C1)；s) = sf + p'dt + JS2Jll ；d; = d) + y(s] + )；sマ=sj + SlCi;: + S2d~+l ;St = K ;以及dt=(i； I其中，第一层小波分解中，<和^>分别为对所述正交分量或同相分量进行奇偶分裂操作得到的奇地址数据和偶地址数据；第二至第七层小波分解中，<和<分别为对上ー层小波分解得到的dt进行奇偶分裂操作得到的奇地址数据和偶地址数据；sf+1表示<向左移动ー个时间单位后得到的数据；足1表示式向右平移一个时间单位后得到的数据；
表示べ向右平移一个时间単位后得到的数据；表示向左延迟ー个时间周期后得到的数据；式2+1表示<向左平移一个时间周期后得到的数据；a、P1. P 2、Y、S2以及る为系数，本发明实施例中的取值请參考前述本发明实施例地层回波信号提取的描述。本发明的实施例结合提升式小波算法与DPSD相敏检波算法的优点，高效地滤除了回波中含有的噪声信号，有效地提取其中的回波信息，更好地反映地层流体信息，将噪声的影响减小到最低程度。本发明的实施例能够大大提高回波信号的提取效率与精确度，为后续更好地进行地层流体分析奠定了良好的基础。本领域的技术人员应该明白，上述的本发明实施例所提供的装置和/或系统的各组成部分，以及方法中的各步骤，可以用通用的计算装置来实现，它们可以集中在单个的计 算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
权利要求
1.一种地层回波信号提取方法，包括如下步骤 采集核磁共振的回波原始信号； 将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量； 对所述正交分量进行小波分解和小波重构，获得正交输出信号； 对所述同相分量进行小波分解和小波重构，获得同相输出信号； 采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。
2.根据权利要求I所述的方法，其中 对所述正交分量进行小波分解和小波重构的步骤，包括对所述正交分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构； 对所述同相分量进行小波分解和小波重构的步骤，包括对所述同相分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。
3.根据权利要求2所述的方法，其中 对所述正交分量或同相分量以db5小波基进行七层小波分解的步骤中，每ー层小波分解均采用如下表达式进行d)=《+ a(s]' + C1);^1 = 5/0+M1+M1-I ； =d)+y{s]+s)_,)； s； = sj + SjdfI2 + 6zdf+l ； A ろ2;以及 4=42,#; 其中，第一层小波分解中，< 和4分别为对所述正交分量或同相分量进行奇偶分裂操作得到的奇地址数据和偶地址数据；第二至第七层小波分解中，4°和ぜ分别为对上ー层小波分解得到的dt进行奇偶分裂操作得到的奇地址数据和偶地址数据； 表示<向左移动ー个时间单位后得到的数据； 表示41向右平移一个时间単位后得到的数据； 表示ろ1向右平移一个时间単位后得到的数据； 式I2表示或2+1向左延迟一个时间周期后得到的数据； 表示42向左平移一个时间周期后得到的数据； a > ^ ^ 2> y > 8 6 2 以及 2,为系数。
4.根据权利要求3所述的方法，其中 a = -I. 5861343420693648 ；^1 = L 0796367753628087 ；^2 = -O. 0529801185718856 ；Y = -0. 8829110755411875 ；S1 = O. 4435068520511142 ； S2=L 576237461483639 ；以及 I = -I. 1496043988602418。
5.一种地层回波信号提取装置，包括 采集模块，用于采集核磁共振的回波原始信号； 运算模块，用于将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量； 分解重构模块，用于对所述正交分量进行小波分解和小波重构，获得正交输出信号；用于对所述同相分量进行小波分解和小波重构，获得同相输出信号； 合成模块，用于采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。
6.根据权利要求5所述的装置，其中 所述分解重构模块用于对所述正交分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构； 所述分解重构模块用于对所述同相分量以db5小波基进行七层小波分解，获得七组细节系数和ー组最底层近似系数；将七组细节系数置为零，利用所述最底层近似系数以及置为零的七组细节系数进行小波重构。
7.根据权利要求6所述的装置，其中 所述分解重构模块用于采用如下表达式进行七层小波分解中的每ー层小波分解
8.根据权利要求7所述的装置，其中 a = -I. 5861343420693648 ；^1 = L 0796367753628087 ；β 2 = -0. 0529801185718856 ； γ = -0.8829110755411875 ；δ1: = 0.4435068520511142 ；δ2 = 1. 576237461483639 ；以及 ξ = -1. 1496043988602418。
全文摘要
本发明公开了一种地层回波信号提取方法及装置，克服目前提取核磁共振的地层回波信号误差较大、精度较低的缺陷，该方法包括采集核磁共振的回波原始信号；将所述回波原始信号与同频率的正弦信号和余弦信号分别进行正交乘法运算，获得所述回波原始信号的正交分量和同相分量；对所述正交分量进行小波分解和小波重构，获得正交输出信号；对所述同相分量进行小波分解和小波重构，获得同相输出信号；采用数字相敏检波技术对所述正交输出信号和同相输出信号进行合成，获得核磁共振回波信号的初始相位和包络信号幅值。本发明可以高效地滤除信号中的噪声干扰。
文档编号G01V3/38GK102830435SQ20121029234
公开日2012年12月19日 申请日期2012年8月16日 优先权日2012年8月16日
发明者张嘉伟, 宋公仆, 薛志波 申请人:中国海洋石油总公司, 中海油田服务股份有限公司