I am currently working on a similar system. But your properties (2) and (3) would not work and need change.
Instead, the following properties would work much better:
$$\int_a^c f(x) dx=\int_a^b f(x)dx+\int_b^c f(x)dx\tag{1}$$
$$\int_a^b (f(x)+g(x)) dx=\int_a^b f(x)dx+\int_a^b g(x)dx\tag{2}$$
$$\int_a^b c f(x) dx =c \int_a^b f(x) dx\tag{3}$$
$$\int_{-\infty}^0 f(x) dx=\int_0^\infty f(-x) dx\tag{4}$$
where $a,b,c,f(x)$ and $g(x)$ take values from ℝ ∪ {−∞, +∞}.
Thus, any integral
$$T=\int_a^b f(x) dx$$
represents an "extended" number.
The integrals which are regularizable by a stable method like Cesaro or Abel (as opposed to Ramanujan or Dirichlet) are taken to be equal to their regularized sum:
$$\int_0^\infty f(x)\,dx=\lim_{\epsilon\to 0}\int_0^\infty e^{-\epsilon x}f(x) \, dx\tag{6}$$
Two integrals $\int_0^\infty f(x) dx$ and $\int_0^\infty g(x) dx$ are thus equal if
$$\lim_{\epsilon\to 0}\int_0^\infty e^{-\epsilon x}(f(x)-g(x)) \, dx=0$$
We also can equate some divergent series to the integrals:
$$\sum_{k=0}^\infty f(k)=\int_{-1/2}^\infty\sum_{k=0}^\infty\operatorname{rect}(x+k)f(k)dx$$
In our notation we will consider by definition $$\sum_{k=n}^\infty f(k)=\sum_{k=0}^\infty f(k)-\sum_{k=0}^{n-1}f(k)$$
Now, we postulate that the regularized value or the integral or corresponding series represents the regular part of an extended number, while the rest is the irregular part. Among the suitable regularization methods are Cesaro, Abel, Ramanujan, Borel, Dirichlet regulaization and some others (they agree with each other when applicable). We will denote the regularized value of an extended number $w$ as $\operatorname{reg} w$
Particularly, very useful would be the Faulhaber's formula for Ramanujan's summation of analytic functions:
$$\operatorname{reg} \sum _{n=0}^{\infty} f(n)= -\sum_{n=1}^{\infty} \frac{f^{(n-1)} (0)}{n!} B_n \tag{7}$$
We will use the following symbols for the three most key integrals and series:
$$\omega_+=\sum_{k=0}^\infty 1$$
$$\omega_-=\sum_{k=1}^\infty 1=\omega_+-1$$
$$\tau=\int_0^\infty dx=\omega_+-1/2=\omega_-+1/2$$
(this can also be formally interpreted as $\tau=\pi\delta(0)$ due to Fourier transform).
By interpreting formula (7) as a Taylor series, we come to a formula which allows to generalize the analytic functions to extended numbers (at least in the sense of determining the regular part of the result):
$$\operatorname{reg} f'(\omega_-+z)= \Delta f(z)\tag{8}$$ and in particular, to the generalizations of powers of our key series:
$$\operatorname{reg}\omega_-^n=B_n\tag{9}$$
$$\operatorname{reg}\omega_+^n=B^*_n\tag{9a}$$
where $B^*$ are the second Bernoulli numbers (those which have $B^*_1=1/2$ ).
A more general formula reveals the role of the Hurwitz zeta function:
$$\operatorname{reg}(\omega_-+x)^a= B_a(x)=-a\zeta(1-a,x)$$
Based on formula (7) we even can derive an expression for a derivative of an analytic function which does not use limits:
$$f'(x)=\operatorname{reg}(f(\omega_++x)-f(\omega_-+x))=\operatorname{reg} \Delta f(\omega_-+x)$$
which works for any regular $x$.
Also, since many series expansions of trigonometric functions use Bernoulli numbers, we can interpret them as regular parts of similar series but involving extended numbers. This way, and using formula (8) we can obtain the following relations:
$$\operatorname{reg}\sin (a\omega_-+x) = \frac{a}{2} \cot \left(\frac{a}{2}\right) \sin x -\frac{a}{2} \cos x$$
$$\operatorname{reg}\cos (a\omega_-+x) = \frac{a}{2} \csc \left(\frac{a}{2}\right) \cos \left(\frac{a}{2}- x \right)$$
$$\operatorname{reg}\ln (\omega_-+z)=\psi(z)$$
$$\operatorname{reg} e^{z\omega_-}=\frac{z}{e^{z}-1}$$
Particularly, $$\operatorname{reg}\sin \omega_-=-1/2;$$ $$\operatorname{reg}\sin \omega_+=1/2;$$ $$\operatorname{reg}\ln \omega_+=-\gamma;$$ $$\operatorname{reg} e^{\omega_-}=\frac{1}{e-1};$$ $$\operatorname{reg} e^{\omega_+}=\frac{e}{e-1}.$$
Another notable thing is the possibility to express trigonometric functions via inverse trigonometric or logarithms:
$$\cot x=\operatorname{reg}\frac1{\pi }\ln \left(\frac{\omega _+-\frac{x}{\pi }}{\omega _-+\frac{x}{\pi }}\right) = \operatorname{reg}\frac2z \cos (2x\omega_\pm)$$
$$\tan x=\operatorname{reg} \frac1\pi\ln \left(\frac{\tau +\frac{x}{\pi }}{\tau -\frac{x}{\pi }}\right)$$
$$\coth x=\operatorname{reg}\frac{1}{\pi} \operatorname{arccoth}\left(\frac{\pi \omega _+}{x}\right)+\frac1x=\operatorname{reg}\frac1x \cosh (2 x\omega_\pm)$$
Following from Faulhaber's formula (for $ n\ge0 $ ),
$$ \int_0^\infty x^n dx=\frac{\left(\tau +\frac{1}{2}\right)^{n+2}-\left(\tau -\frac{1}{2}\right)^{n+2}}{(n+1)(n+2)}=\frac{\omega _+^{n+2}-\omega _-^{n+2}}{(n+1)(n+2)} $$
Interpreting Fourier transform formally, for even $ n $ we also have
$$ \int_0^\infty x^n dx=i^n\pi\delta^{(n)}(0) $$
For $ n>1 $
$$ \int_0^\infty \frac1{x^n} dx=\frac1{(n-1)!}\int_0^\infty x^{n-2} dx=\frac{\omega _+^{n}-\omega _-^{n}}{(n-1)n!} $$
Particularly,
$$\int_0^\infty 1 dx =\tau$$
$$\int_0^\infty x dx=\frac{\tau^2}2+\frac1{24}$$
$$\int_0^\infty x^2 dx=\frac{\tau^3}3+\frac{\tau}{12}$$
We introduce generalized limits in the following way:
$$
\operatorname{gen}\lim_{x\to u^+}f(x)=f(a)-\int_u^a f'(x)dx
$$
where $a>u$ and
$$\operatorname{gen}\lim_{x\to u^-}f(x)=f(a)+\int_a^u f'(x)dx$$
where $a<u$
This value can serve as a measure of the growth rate of a function.
Particularly, for $ n\ge0 $
$$ \operatorname{gen}\lim_{x\to\infty}x^n=\frac{\omega _+^{n+1}-\omega _-^{n+1}}{n+1} $$
For odd $ n $ ,
$$ \operatorname{gen}\lim_{x\to\infty}x^n=i^{n-1}\pi n\delta^{(n-1)}(0) $$
And
$$ \operatorname{gen}\lim_{x\to 0^+} \frac1{x^n}=\frac{\omega _+^{n+1}-\omega _-^{n+1}}{(n+1)!} $$
Following Urs Graf, p.36, for odd $ n $ ,
$$ \operatorname{gen}\lim_{x\to 0^+} \frac1{x^n}=\frac{i^{n-1}\pi\delta^{(n-1)}(0)}{(n-1)!} $$
For instance,
$$\int_{0^+}^\infty \frac1{x^2}=\tau$$
This, combined with the delta function form of $\tau=\pi\delta(0)$ can explain why logarithm has its imaginary part represented as a step function (which is integral of delta function).
The Cauchy principal value of an analytic function at a pole corresponds to the regular part of its generalized limits at the pole, while the order of the pole corresponds to the order of polynomial of $\tau$ of the generalized limit. Thus,
$$\operatorname{gen}\lim_{x\to0^\pm}\Gamma(0)=-\gamma\pm\tau$$
$$\operatorname{gen}\lim_{x\to{-1}^\pm}\Gamma(x)=\gamma-1\mp\tau$$
$$\operatorname{gen}\lim_{x\to{-2}^\pm}\Gamma(x)=\frac{3}{4}-\frac{\gamma }{2}\pm\frac\tau 2$$
$$\operatorname{gen}\lim_{x\to{-3}^\pm}\Gamma(x)=\frac{\gamma }{6}-\frac{11}{36}\mp\frac\tau 6$$
$$\operatorname{gen}\lim_{x\to1^\pm}\zeta(x)=\gamma\pm\tau$$
For Bernoulli polynomials we have these noticeable formulas:
$$\omega_-^n=\operatorname{gen}\lim_{x\to\infty}B_n(x)$$
$$\omega_+^n=\operatorname{gen}\lim_{x\to\infty}B_n(x+1)$$
And more generally,
$$(\omega_-+a)^n=B_n(a)+n\int_0^\infty B_{n-1}(x+a)dx$$
Moreover, generalizing the above formulas we can write down the general formula for conversion of divergent integral representation to omega-tau representation:
$$\int_0^{\infty } f(x) \, dx=\operatorname{reg}\int_0^{\infty } f(x) \, dx+\int _{\omega_-}^{\omega_+} \int _0^x f(t) dtdx$$
The integral from $\omega_-$ to $\omega_+$ should be understood as the difference of the antiderivative on those values.
Here is a table of some extended numbers and their representations in different forms:
$$
\begin{array}{cccccc}
\text{Delta form} & \text{In terms of } \tau, \omega_+,\omega_- & \text{Regular part} & \text{Integral or series form} & \text{Generalized limit form} \\
\pi \delta (0) & \tau & 0 & \int_0^{\infty } \, dx;\int_0^{\infty } \frac{1}{x^2} \, dx & \operatorname{gen}\lim_{x\to\infty}x;\operatorname{gen}\lim_{x\to0^+}\frac1x \\
\pi \delta (0)-\frac{1}{2} & \omega _-;\tau-\frac{1}{2} & -\frac{1}{2} & \sum _{k=1}^{\infty } 1 & \operatorname{gen}\lim_{x\to\infty} (x-1/2) \\
\pi \delta (0)+\frac{1}{2} & \omega _+;\tau+\frac{1}{2} & \frac{1}{2} & \sum _{k=0}^{\infty } 1 & \operatorname{gen}\lim_{x\to\infty} (x+1/2) \\
2 \pi \delta (i) & e^{\omega_+}-e^{\omega_-}-1 & 0 & \int_{-\infty }^{\infty } e^x \, dx & \operatorname{gen}\lim_{x\to\infty} e^x \\
& \frac{\tau ^2}{2}+\frac{1}{24};\frac{\omega_+^3-\omega_-^3}6 & 0 & \int_0^{\infty} x \, dx;\int_0^\infty \frac2{x^3}dx & \operatorname{gen}\lim_{x\to\infty}\frac{x^2}2;\operatorname{gen}\lim_{x\to{0^+}} \frac1{x^2}\\
& \frac{\tau ^2}{2}-\frac{1}{24} & -\frac1{12} & \sum _{k=0}^{\infty } k & \operatorname{gen}\lim_{x\to\infty} \left(\frac{x^2}2-\frac1{12}\right) \\
-\pi \delta''(0) &\frac {\tau^3}3 +\frac\tau{12};\frac{\omega_+^4-\omega_-^4}{12}& 0 & \int_0^\infty x^2dx;\int_0^\infty\frac6{x^4}dx&\operatorname{gen}\lim_{x\to\infty}\frac{x^3}3;\operatorname{gen}\lim_{x\to{0^+}} \frac2{x^3}\\
\pi^2\delta(0)^2-\pi\delta(0)+1/4&\omega_-^2&\frac16&2 \int_0^{\infty } \left(x-\frac{1}{2}\right) \, dx+\frac{1}{6}&\operatorname{gen}\lim_{x\to\infty}B_2(x)\\
\pi^2\delta(0)^2+\pi\delta(0)+1/4&\omega_+^2&\frac16&2 \int_0^{\infty } \left(x+\frac{1}{2}\right) \, dx+\frac{1}{6}&\operatorname{gen}\lim_{x\to\infty}B_2(x+1)\\
\pi^2\delta(0)^2&\tau^2&-\frac1{12}&\int_{-\infty}^{\infty } |x| \, dx-\frac{1}{12}&\operatorname{gen}\lim_{x\to\infty}B_2(x+1/2)\\
&\ln \omega_++\gamma&0&\int_1^\infty \frac{dx}x;\sum_{k=1}^\infty \frac1x -\gamma&\operatorname{gen}\lim_{x\to\infty}\ln x\\
-\pi\delta''(0)-\frac14 \pi\delta(0);\pi^3\delta(0)^3&\tau^3&0&\int_0^\infty \left(3x^2-\frac1{4}\right)dx&\operatorname{gen}\lim_{x\to\infty}B_3(x+1/2)\\
\frac{2\pi\delta(i)+1}{e-1}&e^{\omega_-}&\frac1{e-1}&\frac1{e-1}+\frac1{e-1}\int_{-\infty}^\infty e^x dx&\operatorname{gen}\lim_{x\to\infty} \frac{e^x+1}{e-1}\\
\frac{2\pi\delta(i)+1}{1-e^{-1}}&e^{\omega_+}&\frac1{1-e^{-1}}&\frac1{1-e^{-1}}+\frac1{1-e^{-1}}\int_{-\infty}^\infty e^x dx&\operatorname{gen}\lim_{x\to\infty} \frac{e^x+1}{1-e^{-1}}\\
\end{array}
$$
The system forms an integral domain which can be extended to form a field.
Still, there are many white areas, particularly, there is no straightforward way to transform a divergent integral into an expression in terms of $\omega_\pm$ and vice versa. There is no formula to construct an integral representing a product of two integrals, etc.