Question: Is it right that any domain enclosed by a closed surface $\Sigma\subset\Omega$ is contained within $\Omega$ is equivalent to $H^2(\Omega)=0$. If so we have the following theorem for this domain.
Here we convention that if $\Sigma_1$ and $\Sigma_2$ have the same boundary, then the orientation of the boundary is also the same.
Let $\Omega\subset\in\mathbb R^3$ be a domain, and $P,Q,R\in C^1(G)$. If the de Rham cohomology $H^2(G)=0$. Then the following three statements are equivalence.
$(1)$ Surface integer $\displaystyle{\iint\limits_{\Sigma}P dydz+Q dzdx+R dxdy}$ only depends on the boundary of $\Omega$.
$(2)$ For every closed surface $\Sigma\in\Omega$, we have that $\displaystyle{\iint\limits_{\Sigma}P dydz+Q dzdx+R dxdy=0}$.
$(3)$ $\displaystyle{\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}=0}$.
Here $(1)\Longrightarrow (2)\Longrightarrow(3)$ is easy. We only prove $(3)\Longrightarrow(1)$. Since $H^2(G)=0$, every closed form of $G$ is exact.
Since $\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}=0$, there exists vector field $\mathbf{F}=(f_1,f_2,f_3)$ such that
$$\mathrm{rot}(f_1,f_2,f_3)=(P,Q,R),$$
that is
$$P=\frac{\partial f_3}{\partial y}-\frac{\partial f_2}{\partial z},\quad
Q=\frac{\partial f_1}{\partial z}-\frac{\partial f_3}{\partial x},\quad
R=\frac{\partial f_2}{\partial x}-\frac{\partial f_1}{\partial y}.$$
Hence for every oriented surface $\Sigma\subset G$, by Stokes's theorem we have that
\begin{eqnarray*}
% \nonumber to remove numbering (before each equation)
&& \iint\limits_{\Sigma}P dydz+Q dzdx+R dxdy\\
&=& \iint\limits_{\Sigma}\left(\frac{\partial f_3}{\partial y}-\frac{\partial f_2}{\partial z}\right) dydz+
\left(\frac{\partial f_1}{\partial z}-\frac{\partial f_3}{\partial x}\right) dzdx+\left(\frac{\partial f_2}{\partial x}-\frac{\partial f_1}{\partial y}\right) dxdy \\
&=& \oint_{\partial\Sigma}f_1dx+f_2dy+f_3dz
\end{eqnarray*}
Now for every two surfaces $\Sigma_1,\Sigma_2$, if
$$\partial\Sigma_1=\Gamma=\partial\Sigma_2,$$
then
$$\iint\limits_{\Sigma_1}P dydz+Q dzdx+R dxdy
=\oint_{\Gamma}f_1dx+f_2dy+f_3dz=\iint\limits_{\Sigma_2}P dydz+Q dzdx+R dxdy.$$