METHOD FOR EVALUATING THICKNESS AND DENSITY OF ADSORBED METHANE IN PORES CONTRIBUTED BY ORGANIC MATTER, CLAY AND OTHER MINERALS IN MUD SHALE RESERVOIR

Abstract

A method for evaluating thickness and density of adsorbed methane in pores contributed by organic matter, clay and other minerals in a mud shale reservoir, including: crushing a sample and selecting three or more subsamples with different meshes to determine TOC, kerogen, whole rock analysis, low-temperature nitrogen adsorption-desorption and methane isotherm adsorption; calculating contents of organic matter in respective subsamples from TOC and kerogen contents; normalizing contents of organic matter, clay and other minerals; evaluating the volume of pores contributed by organic matter, clay and other minerals per unit mass according to contents thereof and low-temperature nitrogen adsorption-desorption; evaluating content of adsorbed methane in organic matter, clay and other minerals per unit mass according to contents thereof and methane isotherm adsorption; and establishing a model for calculating density and thickness of adsorbed methane.

Claims

1. A method for evaluating thickness and density of adsorbed methane in pores contributed by organic matter, clay and other minerals in a mud shale reservoir, comprising: 1) crushing a mud shale reservoir sample to produce a plurality of subsamples; and selecting three or more subsamples varying in mesh for determinations of organic carbon content and kerogen content, whole rock analysis, and determinations of low temperature nitrogen adsorption-desorption and methane isotherm adsorption; wherein: mass percentages of organic carbon in respective subsamples are w.sub.TOC1.sup.0, w.sub.TOC2.sup.0, . . . and w.sub.TOCn.sup.0 (%), respectively; mass percentages of carbon in kerogen in respective subsamples are w.sub.C1, w.sub.C2, and w.sub.Cn (%), respectively; mass percentages of clay in respective subsamples are w.sub.clay1.sup.0, w.sub.clay2.sup.0, . . . and w.sub.clayn.sup.0 (%), respectively; and mass percentages of other minerals in respective subsamples are w.sub.others1.sup.0, w.sub.others2.sup.0, . . . and w.sub.othersn.sup.0 (%), respectively, pores in respective subsamples per unit mass having a size respectively of <2 nm, 2-5 nm, 5-10 nm, 10-20 nm, 20-50 nm, 50-100 nm and 100-200 nm have a volume of V.sub.ij (cm.sup.3/ g); respective subsamples per unit mass have an adsorbed methane content of Q.sub.ixy (m.sup.3/t) under a temperature of T.sub.x and a pressure of P.sub.y, wherein i is the number of respective subsamples of the mud shale reservoir, and is selected from 1, 2, 3, . . . , n; j is the number of pore sizes, selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; 2) substituting the mass percentages of organic carbon (w.sub.TOC1.sup.0, w.sub.TOC2.sup.0, . . . and w.sub.TOCn.sup.0) and the corresponding mass percentages of carbon in kerogen (w.sub.C1, w.sub.C2, . . . and w.sub.Cn) in respective subsamples into the following equation to obtain mass percentages of organic matter in respective subsamples (w.sub.TOM1.sup.0, w.sub.TOM2.sup.0, . . . and w.sub.TOMn.sup.0);
w.sub.TOMi.sup.0=w.sub.TOMi.sup.0/w.sub.Ci100%; wherein w.sub.TOMi.sup.0 (%) is an unnormalized mass percentage of organic matter in respective subsamples; w.sub.TOCi.sup.0 (%) is an experimentally measured mass percentage of organic carbon in respective subsamples; w.sub.Ci (%) is an experimentally measured mass percentage of carbon in kerogen in respective subsamples; i is the number of respective subsamples of the mud shale reservoir, and is selected from 1, 2, 3, . . . , n; and normalizing the mass percentages of organic matter, clay and other minerals in respective subsamples according to the following equations; wherein a sum of the mass percentages of organic matter, clay and other minerals in respective subsamples is 100%; the normalized mass percentages of organic matter, clay and other minerals is in respective subsamples are respectively w.sub.TOMi (%), w.sub.clayi (%) and w.sub.othersi (%);
w.sub.TOMi=w.sub.TOMi.sup.0100%
w.sub.clayi=w.sub.clayi.sup.0(100w.sub.TOMi.sup.0)/100%
w.sub.othersi=w.sub.othersi.sup.0(100w.sub.TOMi.sup.0)/100% wherein w.sub.TOMi.sup.0, w.sub.clayi.sup.0 and w.sub.othersi.sup.0 are mass percentages of organic matter, clay and other minerals in respective subsamples before normalization, respectively; i is the number of respective subsamples of the mud shale reservoir, and is selected from 1, 2, 3, . . . , n; 3) establishing a first equation set and a first target function according to the normalized mass percentages of organic matter (w.sub.TOM1, w.sub.TOM2, . . . and w.sub.TOMn), the normalized mass percentages of clay (w.sub.clay1, w.sub.clay2, . . . and w.sub.clayn), and the normalized mass percentage of other minerals (w.sub.others1, w.sub.others2, . . . and w.sub.othersn) in respective subsamples obtained in step (2) and the volume V.sub.ij (cm.sup.3/g) of pores having a size respectively of <2 nm, 2-5 nm, 5-10 nm, 10-20 nm, 20-50 nm, 50-100 nm and 100-200 nm in respective subsamples per unit mass obtained in step (1); wherein in the case of a minimum value of the first target function f(V.sub.TOMj, V.sub.clayj, V.sub.othersj), a volume of pores having a size numbered as j contributed by organic matter per unit mass is V.sub.TOMj (cm.sup.3/g), a volume of pores having a size numbered as j contributed by clay per unit mass is V.sub.clayj, (cm.sup.3/g), and a volume of pores having a size numbered as j contributed by other minerals per unit mass is V.sub.othersj (cm.sup.3/g); $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.1}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.1}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.1}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{1.Math.j}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.2}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{2.Math.j}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.3}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.3}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.3}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{13}$ $.Math.\mathrm{.Math.}.Math..Math.\mathrm{.Math.}.Math..Math.\mathrm{.Math.}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.n}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.n}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.n}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{\mathrm{nj}}$ $.Math.V_{\mathrm{TOM}.Math.\text{-}.Math.j}>0,V_{\mathrm{clay}.Math.\text{-}.Math.j}>0,V_{\mathrm{others}.Math.\text{-}.Math.j}>0$ $f\left(V_{\mathrm{TOM}.Math.\text{-}.Math.j},V_{\mathrm{clay}.Math.\text{-}.Math.j},V_{\mathrm{others}.Math.\text{-}.Math.j}\right)=\underset{i=1}{\overset{n}{.Math.}}.Math..Math.{\left(V_{\mathrm{ij}}-w_{\mathrm{TOM}.Math.\text{-}.Math.i}{V}_{\mathrm{TOM}.Math.\text{-}.Math.i}-w_{\mathrm{clay}.Math.\text{-}.Math.i}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}-w_{\mathrm{others}.Math.\text{-}.Math.i}{V}_{\mathrm{others}.Math.\text{-}.Math.j}\right)}^2;$ wherein V.sub.TOMj, V.sub.clayj and V.sub.othersj are the volumes of pores having a size numbered as j respectively contributed by organic matter, clay and other minerals per unit mass; j is a number of pore size from small to large, and is selected from 1, 2, . . . 6 and 7; i is the number of respective subsamples of the mud shale reservoir, and is selected from 1, 2, 3, . . . , n; 4) establishing a second equation set and a second target function using the normalized mass percentage of organic matter (w.sub.TOM1, w.sub.TOM2, . . . and w.sub.TOMn), the normalized mass percentage of clay (w.sub.clay1, w.sub.clay2, . . . and w.sub.clayn), and the normalized mass percentage of other minerals (w.sub.others1, w.sub.others2, . . . and w.sub.othersn) in respective subsamples obtained in step (2) and the content Q.sub.ixy of adsorbed methane existing in respective subsamples per unit mass obtained in step (1) under a temperature of T.sub.x and a pressure of P.sub.y; wherein in the case of a minimum value of the second target function f(Q.sub.TOMxy, Q.sub.clayxy, Q.sub.othersxy), a temperature of T.sub.x and a pressure of P.sub.y, a content of adsorbed methane existing in organic matter per unit mass is Q.sub.TOMxy, a content of adsorbed methane existing in clay per unit mass is Q.sub.clayxy, and a content of adsorbed methane existing in other minerals per unit mass is Q.sub.othersxy; $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.1}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.1}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.1}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{1.Math.\mathrm{xy}}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.2}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.2}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.2}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{2.Math.\mathrm{xy}}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.3}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.3}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.3}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{3.Math.\mathrm{xy}}$ $.Math.\mathrm{.Math.}.Math..Math.\mathrm{.Math.}.Math..Math.\mathrm{.Math.}$ $.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.n}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.n}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.n}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{\mathrm{nxy}}$ $.Math.Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}>0,Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}>0,Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}>0$ $f\left(Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}},Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}},Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}\right)=\underset{i=1}{\overset{n}{.Math.}}.Math..Math.{\left(Q_{\mathrm{ixy}}-w_{\mathrm{TOM}.Math.\text{-}.Math.i}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}-w_{\mathrm{clay}.Math.\text{-}.Math.i}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}-w_{\mathrm{others}.Math.\text{-}.Math.i}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}\right)}^2$ wherein, Q.sub.TOMxy (m.sup.3/t), Q.sub.clayxy (m.sup.3/t) and Q.sub.othersxy (m.sup.3/t) are the contents of adsorbed methane respectively existing in organic matter, clay and other minerals per unit mass under a temperature of T.sub.x ( C.) and a pressure P.sub.y (MPa); Q.sub.ixy (m.sup.3/t) is a content of adsorbed methane existing in subsample i per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; wherein i is the number of respective subsamples of the mud shale reservoir, and is selected from 1, 2, 3, . . . , n; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; and y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; 5) calculating V.sub.absorbed by TOMjxy according to the following equations based on step (3) by approximating the pores contributed by organic matter per unit mass to cylinders with corresponding pore size; wherein V.sub.absorbed by TOMjxy is the volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass have a size D.sub.TOMj lower than 0.38 nm, V.sub.absorbed by TOMjxy=0; when the D.sub.TOMj is not more than twice a thickness h.sub.absorbed by TOMjxy of adsorbed methane and is not less than 0.38 nm, V.sub.absorbed by TOMjxy=V.sub.TOMj; and when the D.sub.TOMj is more than twice the thickness h.sub.absorbed by TOMjxy of adsorbed methane and is not less than 0.38 nm, $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{TOM}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{TOM}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j};$ $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{cc}0.Math.& \left(D_{\mathrm{TOM}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right).Math.\\ V_{\mathrm{TOM}.Math.\text{-}.Math.j}.Math.& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{TOM}.Math.\text{-}.Math.j},D_{\mathrm{TOM}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\\ \frac{4.Math.D_{\mathrm{TOM}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{TOM}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{TOM}.Math.\text{-}.Math.j},D_{\mathrm{TOM}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array};$ wherein, V.sub.absorbed by TOMjxy (cm.sup.3/g) is the volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.TOMj (cm.sup.3/g) is the volume of pores numbered j contributed by organic matter per unit mass, D.sub.TOMj (nm) is the size of pores numbered j contributed by organic matter; h.sub.absorbed by TOMjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by organic matter; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; and y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; calculating V.sub.absorbed by clayjxy according to the following equations based on step (3) by approximating the pores contributed by clay per unit mass to cylinders with corresponding pore size; wherein V.sub.absorbed by clayjxy is the volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass have a size D.sub.clayj lower than 0.38 nm, V.sub.absorbed by clayjxy=0; when the D.sub.clayj is not more than twice a thickness h.sub.absorbed by clayjxy of adsorbed methane and is not less than 0.38 nm, V.sub.absorbed by clayjxy=V.sub.clayj; and when the D.sub.clayj is more than twice the thickness h.sub.absorbed by clayjxy of adsorbed methane and is not less than 0.38 nm. $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{clay}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{clay}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j};$ $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{cc}0.Math.& \left(D_{\mathrm{clay}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right).Math.\\ V_{\mathrm{clay}.Math.\text{-}.Math.j}.Math.& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{clay}.Math.\text{-}.Math.j},D_{\mathrm{clay}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\\ \frac{4.Math.D_{\mathrm{clay}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{clay}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{clay}.Math.\text{-}.Math.j},D_{\mathrm{clay}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array};$ wherein, V.sub.absorbed by clayjxy (cm.sup.3/g) is the volume of pores occupied by adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.clayj (cm.sup.3/g) is the volume of pores numbered j contributed by clay per unit mass; D.sub.clayj (nm) is the size of pores numbered j contributed by clay per unit mass; h.sub.absorbed by clayjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by clay; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7, x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; and y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; calculating V.sub.absorbed by othersjxy according to the following equations based on step (3) by approximating the pores contributed by other minerals per unit mass to cylinders with corresponding pore size; wherein V.sub.absorbed by othersjxy is the volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass have a size D.sub.othersj lower than 0.38 nm, V.sub.absorbed by othersjxy=0; when the D.sub.othersj is not more than twice a thickness h.sub.absorbed by othersjxy of adsorbed methane and is not less than 0.38 nm, V.sub.absorbed by othersjxy=V.sub.othersj; and when the D.sub.othersj is more than twice the thickness h.sub.absorbed by othersjxy of adsorbed methane and is not less than 0.38 nm, $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{others}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{others}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{others}.Math.\text{-}.Math.j};$ $V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{cc}0.Math.& \left(D_{\mathrm{others}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right).Math.\\ V_{\mathrm{others}.Math.\text{-}.Math.j}.Math.& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{others}.Math.\text{-}.Math.j},D_{\mathrm{others}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\\ \frac{4.Math.D_{\mathrm{others}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{others}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{others}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{others}.Math.\text{-}.Math.j},D_{\mathrm{others}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array};$ wherein , V.sub.absorbed by othersjxy (cm.sup.3/g) is the volume of pores occupied by adsorbed methane existing in pores numbered j contributed by other minerals per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.othersj (cm.sup.3/g) is the volume of pores numbered j contributed by other minerals per unit mass; D.sub.othersj (nm) is the size of pores numbered j contributed by other minerals per unit mass; h.sub.absorbed by othersjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by other minerals; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; and y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; 6) establishing a third equation set and a third target function based on steps (4) and (5) according to the facts that a density of adsorbed methane is lower than that of solid methane but greater than that of free methane; the density of adsorbed methane in pores of organic matter decreases with the increase of pore size, and the density of adsorbed methane decreases with the increase of temperature while increases with the increase of pressure; wherein in the case of a minimum value of the third target function f(.sub.absorbed by TOMjxy, h.sub.absorbed by TOMjxy), a density .sub.absorbed by TOMjxy and a thickness h.sub.absorbed by TOMjxy of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y can be obtained; $\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]=Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.7.Math.\mathrm{jxy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jmy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{my}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxz}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.7}>{}_{\mathrm{free}.Math.\text{-}.Math.x.Math..Math.1}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.7.Math.\mathrm{xy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jmy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxz}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.1}$ $f\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}},h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\right)=\underset{x=1}{\overset{m}{.Math.}}.Math..Math.\underset{y=1}{\overset{z}{.Math.}}.Math..Math.{\left(Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}-\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]\right)}^2;$ wherein, V.sub.absorbed by TOMjxy (cm.sup.3/g) is the volume of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.absorbed by TOMjxy (kg/m.sup.3) is a density of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.freexy (kg/m.sup.3) is a density of free methane under a temperature of T.sub.x and a pressure of P.sub.y; Q.sub.TOMxy (m.sup.3/t) is the content of adsorbed methane existing in organic matter per unit mass; .sub.solid (kg/m.sup.3) is a density of solid methane; h.sub.absorbed by TOMjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by organic matter; M is the molar mass of methane referring to 16.0425 g/mol; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; and y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; establishing a forth equation set and a forth target function based on steps (4) and (5) according to the facts that the density of adsorbed methane is lower than that of solid methane but greater than that of free methane, the density of adsorbed methane in the pores of clay decreases with the increase of pore size, and the density of adsorbed methane decreases with the increase of temperature while increases with the increase of pressure; wherein in the case of a minimum value of the forth target function f(.sub.absorbed by clayjxy, h.sub.absorbed by clayjxy), a density .sub.absorbed by clayjxy and a thickness h.sub.absorbed by clayjxy of adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y can be obtained; $\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]=Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.7.Math.\mathrm{jxy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jmy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{my}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxz}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.7}>{}_{\mathrm{free}.Math.\text{-}.Math.x.Math..Math.1}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.7.Math.\mathrm{xy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jmy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxz}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.1}$ $f\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}},h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\right)=\underset{x=1}{\overset{m}{.Math.}}.Math..Math.\underset{y=1}{\overset{z}{.Math.}}.Math..Math.{\left(Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}-\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]\right)}^2;$ wherein, V.sub.absorbed by clayjxy (cm.sup.3/g) is the volume of adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.absorbed by clayjxy (kg/m.sup.3) is a density of adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.freexy (kg/m.sup.3) is the density of free methane under a temperature of T.sub.x and a pressure of P.sub.y; Q.sub.clayxy (m.sup.3/t) is the content of adsorbed methane existing in clay per unit mass; .sub.solid (kg/m.sup.3) is the density of solid methane; h.sub.absorbed by clayjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by clay; M is the molar mass of methane referring to 16.0425 g/mol; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; y is the number of pressure from low to high, and is selected from 1, 2, . . . , z; establishing a fifth equation set and a fifth target function based on steps (4) and (5) according to the facts that the density of adsorbed methane is lower than that of solid methane but greater than that of free methane, the density of adsorbed methane in the pores of other minerals decreases with the increase of pore size, and the density of adsorbed methane decreases with the increase of temperature while increases with the increase of pressure; wherein in the case of a minimum value of the fifth target function f(.sub.absorbed by othersjxy, h.sub.absorbed by othersjxy), a density .sub.absorbed by othersjxy and a thickness h.sub.absorbed by othersjxy of adsorbed methane in pores numbered j contributed by other minerals per unit mass under a temperature of T.sub.x and a pressure of P.sub.y can be obtained; $\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]=Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.7.Math.\mathrm{jxy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jmy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{my}}$ ${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxz}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.7}>{}_{\mathrm{free}.Math.\text{-}.Math.x.Math..Math.1}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.7.Math.\mathrm{xy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jmy}}$ $h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxz}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.1}$ $f\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}},h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}\right)=\underset{x=1}{\overset{m}{.Math.}}.Math..Math.\underset{y=1}{\overset{z}{.Math.}}.Math..Math.{\left(Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}-\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]\right)}^2;$ wherein, V.sub.absorbed by othersjxy (cm.sup.3/g) is the volume of adsorbed methane in pores numbered j contributed by other minerals per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.absorbed by othersjxy (kg/m.sup.3) is a density of adsorbed methane in pores numbered j contributed by other minerals per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.freexy (kg/m.sup.3) is the density of free methane under a temperature of T.sub.x and a pressure of P.sub.y; Q.sub.othersxy (m.sup.3/t) is the content of adsorbed methane existing in other minerals per unit mass; .sub.solid (kg/m.sup.3) is the density of solid methane; h.sub.absorbed by othersjxy (nm) is the thickness of adsorbed methane in pores numbered j contributed by other minerals; M is the molar mass of methane referring to 16.0425 g/mol; j is the number of pore sizes from small to large, and is selected from 1, 2, . . . , 7; x is the number of temperature from low to high, and is selected from 1, 2, . . . , m; y is the number of pressure from low to high, and is selected from 1, 2, . . . , z.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The FIGURE is a flow chart showing the method of the invention for evaluating thickness and density of adsorbed methane in pores contributed by organic matter, clay and other minerals in a mud shale reservoir.

DETAILED DESCRIPTION OF EMBODIMENTS

Example 1

[0042] As shown in the FIGURE, the invention provided a method for evaluating thickness and density of adsorbed methane in pores contributed by organic matter, clay and other minerals in a mud shale reservoir, which was described as follows.

[0043] 1) A mud shale reservoir sample was crushed into a plurality of subsamples, of which 5 subsamples respectively of 20-40 mesh, 40-60 mesh, 60-80 mesh, 80-100 mesh and 100-120 mesh were selected for determinations of TOC content and kerogen content, whole rock analysis, and analysis of low temperature nitrogen adsorption-desorption and methane isotherm adsorption. The obtained mass percentages of TOC in respective subsamples were 1.28%, 1.10%, 2.07%, 2.22% and 2.94%, respectively; the obtained mass percentages of carbon in kerogen were 86.12%, 86.72%, 87.01%, 85.57% and 87.98%, respectively; the obtained mass percentages of clay were 41.6%, 42.2%, 23.0%, 25.7% and 30.3%, respectively; and the obtained mass percentages of other minerals were 58.4%, 57.8%, 77.0%, 74.3% and 69.7%, respectively. The obtained volume V.sub.ij (cm.sup.3/g) of pores in respective subsamples per unit mass having a size respectively of <2 nm, 2-5 nm, 5-10 nm, 10-20 nm, 20-50 nm, 50-100 nm and 100-200 nm in the low temperature nitrogen adsorption-desorption was shown in Table 1. After the methane isotherm adsorption, the obtained adsorbed methane content of Q.sub.ixy (m.sup.3/t) in respective subsamples per unit mass under 30 C. and 1 MPa, 2 MPa, 3 MPa, 4 MPa, 5 MPa, 6 MPa, 7 MPa, 8 MPa, 9 MPa and 10 MPa was shown in Table 2.

TABLE-US-00001 TABLE 1 Volume of Pores with Different Sizes in Respective Subsamples Per Unit Mass Having (10.sup.3 cm.sup.3/g) Subsample Size No. <2 nm 2-5 nm 5-10 nm 10-20 nm 20-50 nm 50-100 nm 100-200 nm 1 0.29 1.16 0.84 1.09 0.74 0.44 0.23 2 0.34 1.16 0.84 0.98 0.67 0.37 0.19 3 0.17 1.07 0.79 1.01 0.71 0.39 0.19 4 0.31 1.09 0.88 1.07 0.70 0.36 0.17 5 0.30 1.34 1.16 1.21 0.88 0.58 0.26 6 0.33 1.01 0.74 0.87 0.66 0.39 0.14 7 0.29 1.26 1.03 1.26 0.79 0.49 0.22 8 0.31 1.26 1.05 1.27 0.85 0.46 0.23

TABLE-US-00002 TABLE 2 Adsorbed Methane Content in Respective Subsamples Per Unit Mass under 30 C. and Different Pressures (m.sup.3/t) Subsample Pressure No. 1 MPa 2 MPa 3 MPa 4 MPa 5 MPa 6 MPa 7 MPa 8 MPa 9 MPa 10 MPa 1 1.01 1.35 1.52 1.62 1.68 1.73 1.77 1.80 1.82 1.84 2 0.97 1.30 1.46 1.56 1.62 1.67 1.70 1.73 1.75 1.77 3 1.17 1.56 1.76 1.88 1.95 2.01 2.05 2.08 2.11 2.13 4 1.21 1.61 1.82 1.94 2.02 2.07 2.12 2.15 2.18 2.20 5 1.24 1.66 1.86 1.99 2.07 2.13 2.17 2.21 2.23 2.26

[0044] 2) The mass percentages of organic matter without normalization in 5 subsamples (1.49%, 1.27%, 2.38%, 2.59% and 3.34%) were obtained by substituting the mass percentage of TOC in 5 subsamples (1.28%, 1.10%, 2.07%, 2.22% and 2.94%), the mass percentage of carbon in kerogen (86.12%, 86.72%, 87.01%, 85.57% and 87.98%) into the following equation.

w.sub.TOMi.sup.0=w.sub.TOCi.sup.0/w.sub.Ci100%;

[0045] where w.sub.TOMi.sup.0 (%) was an unnormalized mass percentage of organic matter in respective subsamples; w.sub.TOCi.sup.0 (%) was an experimentally measured mass percentage of organic carbon in respective subsamples; w.sub.Ci (%) was an experimentally measured mass percentage of carbon in kerogen in respective subsamples; i was the number of respective subsamples of the mud shale reservoir, and was selected from 1, 2, 3, . . . , n.

[0046] Then the mass percentages of organic matter, clay and other minerals in respective subsamples were normalized according to the following equations, where a sum of the mass percentages of organic matter, clay and other minerals in respective subsamples was 100%. The obtained normalized mass percentages of organic matter were respectively 1.49%, 1.27%, 2.38%, 2.59% and 3.34%; the obtained normalized mass percentages of clay were respectively 40.98%, 41.67%, 22.47%, 25.05% and 29.32%; and the obtained normalized mass percentages of other minerals were respectively 57.53%, 57.08%, 75.21%, 72.42% and 67.45% in 5 subsamples.

w.sub.TOMi=w.sub.TOMi.sup.0100%

w.sub.clayi=w.sub.clayi.sup.0(100w.sub.TOMi.sup.0)/100%

w.sub.othersi=w.sub.othersi.sup.0(100w.sub.TOMi.sup.0)/100%

[0047] where w.sub.TOMi (%), w.sub.clayi (%) and w.sub.othersi (%) were normalized mass percentages of organic matter, clay and other minerals in respective subsamples, respectively; w.sub.TOMi.sup.0, w.sub.clayi.sup.0 and w.sub.othersi.sup.0 were mass percentages of organic matter, clay and other minerals in respective subsamples before normalization, respectively; i was the number of respective subsamples of the mud shale reservoir, and was selected from 1, 2, 3, . . . , n.

[0048] 3) A first equation set and a first target function were established according to the normalized mass percentages of organic carbon (1.49%, 1.27%, 2.38%, 2.59% and 3.34%), the normalized mass percentages of clay (40.98%, 41.67%, 22.47%, 25.05% and 29.32%), and the normalized mass percentages of other minerals (57.53%, 57.08%, 75.21%, 72.42% and 67.45%) in respective subsamples obtained in step (2) and the volume V.sub.ij (referring to Table 1) of pores having a size respectively of <2 nm, 2-5 nm, 5-10 nm, 10-20 nm, 20-50 nm, 50-100 nm and 100-200 nm in respective sub samples per unit mass obtained in step (1),

[0049] where in the case of a minimum value of the first target function f(V.sub.TOMj, V.sub.clayj, V.sub.othersj), a volume of pores having a size numbered as j contributed by organic matter per unit mass was V.sub.TOMj (cm.sup.3/g), a volume of pores having a size numbered as j contributed by clay per unit mass was V.sub.clayj (cm.sup.3/g), and a volume of pores having a size numbered as j contributed by other minerals per unit mass was V.sub.othersj (cm.sup.3/g). The results were shown in Table 3.

[00009]$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.1}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.1}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.1}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{1.Math.j}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.2}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{2.Math.j}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.3}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.3}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.3}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{3.Math.j}$$.Math.\begin{array}{ccc}\mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\end{array}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.n}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}+w_{\mathrm{clay}.Math.\text{-}.Math.n}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}+w_{\mathrm{others}.Math.\text{-}.Math.n}{V}_{\mathrm{others}.Math.\text{-}.Math.j}=V_{\mathrm{nj}}$$.Math.V_{\mathrm{TOM}.Math.\text{-}.Math.j}>0,V_{\mathrm{clay}.Math.\text{-}.Math.j}>0,V_{\mathrm{others}.Math.\text{-}.Math.j}>0$$f\left(V_{\mathrm{TOM}.Math.\text{-}.Math.j},V_{\mathrm{clay}.Math.\text{-}.Math.j},V_{\mathrm{others}.Math.\text{-}.Math.j}\right)={.Math.}_{i=1}^n.Math.{\left(V_{\mathrm{ij}}-w_{\mathrm{TOM}.Math.\text{-}.Math.i}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}-w_{\mathrm{clay}.Math.\text{-}.Math.i}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}-w_{\mathrm{others}.Math.\text{-}.Math.i}{V}_{\mathrm{others}.Math.\text{-}.Math.j}\right)}^2$

[0050] where V.sub.TOMj, V.sub.clayj, V.sub.othersj were the volumes of pores having a size numbered as j respectively contributed by organic matter, clay and other minerals per unit mass; j was a number of pore size from small to large, and was selected from 1, 2, . . . , 7; i was the number of respective subsamples of the mud shale reservoir, and was selected from 1, 2, 3, . . . , 5.

TABLE-US-00003 TABLE 3 Volume of Pores Having Different Sizes (10.sup.3 cm.sup.3/g) Size Component <2 nm 2-5 nm 5-10 nm 10-20 nm 20-50 nm 50-100 nm 100-200 nm Organic Matter 2.2406 20.5508 19.9608 25.1899 15.9194 10.1544 4.7596 Clay 0.7188 2.0910 1.3509 1.4891 1.0590 0.5302 0.2689 Other Minerals 0.0111 0.0349 0.0201 0.0498 0.0446 0.0301 0.0263

[0051] 4) A second equation set and a second target function were established using the normalized mass percentages of organic matter (1.49%, 1.27%, 2.38%, 2.59% and 3.34%), the normalized mass percentages of clay (40.98%, 41.67%, 22.47%, 25.05% and 29.32%), and the normalized mass percentages of other minerals (57.53%, 57.08%, 75.21%, 72.42% and 67.45%) in respective subsamples in obtained step (2) and the content Q.sub.ixy of adsorbed methane existing in respective subsamples per unit mass (referring to Table 2) obtained in step (1) under a temperature of T.sub.x and a pressure of P.sub.y,

[0052] where in the case of a minimum value of the second target function f(Q.sub.TOMxy, Q.sub.clayxy, Q.sub.othersxy), a temperature of 30 C. and a pressure respectively of 1 MPa, 2 MPa, 3 MPa, 4 MPa, 5 MPa, 6 MPa, 7 MPa, 8 MPa, 9 MPa and 10 MPa, a content of adsorbed methane existing in organic matter per unit mass was Q.sub.TOMxy, a content of adsorbed methane existing in clay per unit mass was Q.sub.clayxy, and a content of adsorbed methane existing in other minerals per unit mass was Q.sub.othersxy. The results were shown in Table 4.

[00010]$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.1}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.1}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.1}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{1.Math.\mathrm{xy}}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.2}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.2}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.2}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{2.Math.\mathrm{xy}}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.3}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.3}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.3}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{3.Math.\mathrm{xy}}$$.Math.\begin{array}{ccc}\mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\end{array}$$.Math.w_{\mathrm{TOM}.Math.\text{-}.Math.n}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{clay}.Math.\text{-}.Math.n}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}+w_{\mathrm{others}.Math.\text{-}.Math.n}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}=Q_{\mathrm{nxy}}$$.Math.Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}>0,Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}>0,Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}>0$$f\left(Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}},Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}},Q_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}\right)={.Math.}_{i=1}^n.Math.{\left(Q_{\mathrm{ixy}}-w_{\mathrm{TOM}.Math.\text{-}.Math.i}{Q}_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}-w_{\mathrm{clay}.Math.\text{-}.Math.i}{Q}_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}-w_{\mathrm{others}.Math.\text{-}.Math.i}{Q}_{\mathrm{others}.Math.\text{-}.Math.\mathrm{xy}}\right)}^2$

[0053] where, Q.sub.TOMxy (m.sup.3/t), Q.sub.clayxy, (m.sup.3/t) and Q.sub.othersxy (m.sup.3/t) were the contents of adsorbed methane respectively existing in organic matter, clay and other minerals per unit mass under a temperature of T.sub.x ( C.) and a pressure P.sub.y (MPa);

[0054] Q.sub.ixy (m.sup.3/t) was a content of adsorbed methane existing in subsample i per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; i was the number of respective subsamples of the mud shale reservoir, and was selected from 1, 2, 3, . . . , 5; x was the number of temperature, and was 1; y was the number of pressure from low to high, and was selected from 1, 2, . . . , 10.

TABLE-US-00004 TABLE 4 Content of Adsorbed Methane Existing in Respective Components Per Unit Mass under Different Pressures (m.sup.3/t) Pressure Component 1 MPa 2 MPa 3 MPa 4 MPa 5 MPa 6 MPa 7 MPa 8 MPa 9 MPa 10 MPa Organic Matter 27.551 40.678 47.540 51.241 53.681 55.490 56.800 57.860 58.751 59.369 Clay 1.371 2.009 2.383 2.651 2.851 2.991 3.111 3.201 3.271 3.339 Other Minerals 0.019 0.026 0.031 0.035 0.037 0.039 0.040 0.041 0.042 0.042

[0055] 5) V.sub.absorbed by TOMjxy was calculated according to the following equations based on step (3) by approximating the pores contributed by organic matter per unit mass to cylinders with corresponding pore size,

[0056] where V.sub.absorbed by TOMjxy was a volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass had a size D.sub.TOMj lower than 0.38 nm, V.sub.absorbed by TOMjxy=0; when the D.sub.TOMj was not more than twice a thickness h.sub.absorbed by TOMjxy of adsorbed methane and was not less than 0.38 nm, V.sub.absorbed by TOMjxy=V.sub.TOMj; and when the D.sub.TOMj was more than twice the thickness h.sub.absorbed by TOMjxy of adsorbed methane and was not less than 0.38 nm,

[00011]$V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{TOM}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{TOM}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}..Math.V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{c}\begin{array}{cc}0& \left(D_{\mathrm{TOM}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \begin{array}{cc}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{TOM}.Math.\text{-}.Math.j},D_{\mathrm{TOM}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \frac{4.Math.D_{\mathrm{TOM}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{TOM}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{TOM}.Math.\text{-}.Math.j}\\ \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{TOM}.Math.\text{-}.Math.j},D_{\mathrm{TOM}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}$

[0057] where, V.sub.absorbed by TOMjxy (cm.sup.3/g) was the volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.TOMj (cm.sup.3/g) was the volume of pores numbered j contributed by organic matter per unit mass; D.sub.TOMj (nm) was the size of pores numbered j contributed by organic matter; h.sub.absorbed by TOMjxy (nm) was the thickness of adsorbed methane in pores numbered j contributed by organic matter; j was the number of pore sizes from small to large, and was selected from 1, 2, . . . , 7; x was the number of temperature, and was 1; y was the number of pressure from low to high, and was selected from 1, 2, . . . , 10.

[0058] V.sub.absorbed by clayjxy was calculated according to the following equations based on step (3) by approximating the pores contributed by clay per unit mass to cylinders with corresponding pore size,

[0059] where V.sub.absorbed by clayjxy was a volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass had a size D.sub.clayj lower than 0.38 nm, V.sub.absorbed by clayjxy=0; when the D.sub.clayj was not more than twice a thickness h.sub.absorbed by clayjxy of adsorbed methane and was not less than 0.38 nm, V.sub.absorbed by clayjxy=V.sub.clayj; and when the D.sub.clayj was more than twice the thickness h.sub.absorbed by clayjxy of adsorbed methane and was not less than 0.38 nm,

[00012]$V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{clay}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{clay}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}..Math.V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{c}\begin{array}{cc}0& \left(D_{\mathrm{clay}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \begin{array}{cc}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{clay}.Math.\text{-}.Math.j},D_{\mathrm{clay}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \frac{4.Math.D_{\mathrm{clay}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{clay}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{clay}.Math.\text{-}.Math.j}\\ \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{clay}.Math.\text{-}.Math.j},D_{\mathrm{clay}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}$

[0060] where, V.sub.absorbed by clayjxy (cm.sup.3/g) was the volume of pores occupied by adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.clayj (cm.sup.3/g) was the volume of pores numbered j contributed by clay per unit mass; D.sub.clayj (nm) was the size of pores numbered j contributed by clay per unit mass; h.sub.absorbed by clayjxy (nm) was the thickness of adsorbed methane in pores numbered j contributed by clay; j was the number of pore sizes from small to large, and was selected from 1, 2, . . . , 7; x was the number of temperature, and was 1; y was the number of pressure from low to high, and was selected from 1, 2, . . . , 10.

[0061] calculating V.sub.absorbed by othersjxy according to the following equations based on step (3) by approximating the pores contributed by other minerals per unit mass to cylinders with corresponding pore size;

[0062] where V.sub.absorbed by othersjxy was a volume of pores occupied by adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; when the pores j contributed by organic matter per unit mass had a size D.sub.othersj lower than 0.38 nm, V.sub.absorbed by othersjxy=0; when the D.sub.othersj was not more than twice a thickness h.sub.absorbed by othersjxy of adsorbed methane and was not less than 0.38 nm, V.sub.absorbed by othersjxy=V.sub.othersj; and when the D.sub.othersj was more than twice the thickness h.sub.absorbed by othersjxy of adsorbed methane and was not less than 0.38 nm.

[00013]$V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}=\frac{4.Math.D_{\mathrm{others}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{others}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{others}.Math.\text{-}.Math.j}..Math.V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}=\{\begin{array}{c}\begin{array}{cc}0& \left(D_{\mathrm{others}.Math.\text{-}.Math.j}<0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \begin{array}{cc}{V}_{\mathrm{others}.Math.\text{-}.Math.j}& \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}{D}_{\mathrm{others}.Math.\text{-}.Math.j},D_{\mathrm{others}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}\\ \frac{4.Math.D_{\mathrm{others}.Math.\text{-}.Math.j}{h}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}-4.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}^2}{D_{\mathrm{others}.Math.\text{-}.Math.j}^2}{V}_{\mathrm{others}.Math.\text{-}.Math.j}\\ \left(2.Math.h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{others}.Math.\text{-}.Math.\mathrm{jxy}}<D_{\mathrm{others}.Math.\text{-}.Math.j},D_{\mathrm{others}.Math.\text{-}.Math.j}0.38.Math..Math.\mathrm{nm}\right)\end{array}$

[0063] where, V.sub.absorbed by othersjxy (cm.sup.3/g) was the volume of pores occupied by adsorbed methane existing in pores numbered j contributed by other minerals per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; V.sub.othersj (cm.sup.3/g) was the volume of pores numbered j contributed by other minerals per unit mass; D.sub.othersj (nm) was the size of pores numbered j contributed by other minerals per unit mass; h.sub.absorbed by othersjxy (nm) was the thickness of adsorbed methane in pores numbered j contributed by other minerals; j was the number of pore sizes from small to large. and was selected from 1, 2, . . . , 7; x was the number of temperature, and was 1; y was the number of pressure from low to high, and was selected from 1, 2, . . . , 10.

[0064] 6) A third equation set and a third target function were established based on steps (4) and (5) according to the facts that a density of adsorbed methane was lower than that of solid methane but greater than that of free methane; the density of adsorbed methane in pores of organic matter decreased with the increase of pore size; and the density of adsorbed methane decreased with the increase of temperature while increased with the increase of pressure, where in the case of a minimum value of the third target function f(.sub.absorbed by TOMjxy, h.sub.absorbed by TOMjxy), a density .sub.absorbed by TOMjxy and a thickness h.sub.absorbed by TOMjxy of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y were obtained. The results were shown in Tables 5 and 6.

[00014]$\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]=Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.7.Math.\mathrm{jxy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jmy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{my}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxz}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.7}>{}_{\mathrm{free}.Math.\text{-}.Math.x.Math..Math.1}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.7.Math.\mathrm{xy}}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jmy}}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxz}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.1}$$f\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}},h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\right)=\underset{x=1}{\overset{m}{.Math.}}.Math..Math.\underset{y=1}{\overset{z}{.Math.}}.Math..Math.(Q_{\mathrm{TOM}.Math.\text{-}.Math.\mathrm{xy}}-{\left(\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{TOM}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]\right)}^2$

[0065] where, V.sub.absorbed by TOMjxy (cm.sup.3/g) was the volume of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.absorbed by TOMjxy (kg/m.sup.3) was a density of adsorbed methane in pores numbered j contributed by organic matter per unit mass under a temperature of T.sub.x and a pressure of P.sub.y; .sub.freexy (kg/m.sup.3) was a density of free methane under a temperature of T.sub.x and a pressure of P.sub.y; Q.sub.TOMxy (m.sup.3/t) was the content of adsorbed methane existing in organic matter per unit mass; .sub.solid (kg/m.sup.3) was a density of solid methane; h.sub.absorbed by TOMjxy (nm) was the thickness of adsorbed methane in pores numbered j contributed by organic matter; M was the molar mass of methane referring to 16.0425 g/mol, j was the number of pore sizes from small to large, and was selected from 1, 2, . . . , 7; x was the number of temperature, and was 1; y was the number of pressure from low to high, and was selected from 1, 2, . . . , 10.

[0066] A forth equation set and a forth target function were established based on steps (4) and (5) according to the facts that a density of adsorbed methane was lower than that of solid methane but greater than that of free methane; the density of adsorbed methane in pores of organic matter decreased with the increase of pore size; and the density of adsorbed methane decreased with the increase of temperature while increased with the increase of pressure,

[0067] where in the case of a minimum value of the forth target function f(.sub.absorbed by clayjxy, h.sub.absorbed by clayjxy), a density .sub.absorbed by clayjxy and a thickness h.sub.absorbed by clayjxy of adsorbed methane in pores numbered j contributed by clay per unit mass under a temperature of T.sub.x and a pressure of P.sub.y were obtained. The results were shown in Tables 5 and 6.

[00015]$\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]=Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.7.Math.\mathrm{jxy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jmy}}>{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{my}}$${}_{\mathrm{solid}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxz}}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>{}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.7}>{}_{\mathrm{free}.Math.\text{-}.Math.x.Math..Math.1}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.1.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.2.Math.\mathrm{xy}}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.6.Math.\mathrm{xy}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.7.Math.\mathrm{xy}}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.1.Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j.Math..Math.2.Math.y}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.j\left(m-1\right).Math.y}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jmy}}$$h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxz}}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}\left(z-1\right)}>\mathrm{.Math.}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.2}>h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jx}.Math..Math.1}$$f\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}},h_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\right)=\underset{x=1}{\overset{m}{.Math.}}.Math..Math.\underset{y=1}{\overset{z}{.Math.}}.Math..Math.(Q_{\mathrm{clay}.Math.\text{-}.Math.\mathrm{xy}}-{\left(\frac{22.4}{M}.Math.\underset{j=1}{\overset{7}{.Math.}}.Math..Math.\left[V_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}\left({}_{\mathrm{absorbed}.Math..Math.\mathrm{by}.Math..Math.\mathrm{clay}.Math.\text{-}.Math.\mathrm{jxy}}-{}_{\mathrm{free}.Math.\text{-}.Math.\mathrm{xy}}\right)\right]\right)}^2$