You need $\lceil d\log_2(10)\rceil$ qubits to store a $d$ decimal digits number.
So you can generate a superposition state of binary representation of a 16 decimal digits number using the following circuit:
from qiskit import QuantumRegister, QuantumCircuit, Aer
import numpy as np
num_of_decimal_digits = 16
num_of_qubits = int(np.ceil(num_of_decimal_digits * np.log2(10)))
qr = QuantumRegister(num_of_qubits, 'q')
circ = QuantumCircuit(num_of_qubits)
circ.h(qr)
circ.measure_all()
And to get the number:
simulator = Aer.get_backend('qasm_simulator')
result = simulator.run(circ, shots = 1).result()
counts = result.get_counts()
bitstring = next(iter(counts))
decimal = int(bitstring, 2)
print('Binary:', bitstring)
print('Decimal:', f'{decimal:016d}')
Note: You need all these qubits to have a superposition. If you just want to generate a random 16-digits number, you can use a single qubit and run the circuit repeatedly (by setting shots = number of bits) then append the result bits.